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Chapter 1 
Introduction 



To grasp the importance of understanding elastic materials one just has to 
take a look around. Many of the things one will see then are "elastic" . These 
include all kinds of fluids (fluid dynamics is a special case of elasticity), as 
well as solid materials. From very large to very small scales elasticity theory 
applies. Phenomena like earth-quakes can be understood as the propagation 
of waves in an elastic medium (see pQ). The very same theory applies to 
complete different field like human medicine: tumor growth (see |50j ) or the 
elastic properties of arteria (see can be considered. 

These fields are covered by classical non-relativistic elasticity theory. 
There the material is identified with a certain region in space and its mo- 
tion is modeled by time- dependent mappings between this region and the 
actual position of the material. Usually one assumes the existence of a flat 
(Riemannian) metric on space. 

See or for an introductory textbook into non-relativistic theory. 
A more mathematical approach is given in jH], |H] or [Tl|. For a (very) 
short but nevertheless enlightening compendium including information on 
global properties we refer to 

Attempts to formulate a relativistic theory of elasticity can be traced back 
as early as 1911 for special relativity (|2E]) and 1916 for general relativity 
Q43J). The promptness of these results (being published only a short time 
after the advent of special (general) relativity) indicates the importance of 
elasticity to the scientific community back then. 

Since these early attempts there have been a number of approaches to 
obtain a relativistic theory of elastic materials. We ought to mention the 
work of Carter and Quintana [5| together with its review by Ehlers ^7j. We 
also refer to the works of Maugin [12], Kijowski and Magli [34^ and Tahvildar- 
Zadeh 

In this thesis we follow the formulation of relativistic theory as given 
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by Beig and Schmidt [5] or Christodoulou [T2] . This means that we treat 
relativistic elasto-dynamics as a Lagrangian field theory. We derive the cor- 
responding equations of motion and prove local well-posedness for a variety 
of different setups. 

The following are the basic ingredients: 

• a space-time manifold equipped with a smooth Lorentzian metric 

• a three-dimensional "material manifold" equipped with a smooth vol- 
ume form represents the material 

• a stored energy function, which plays the role of an "equation of state" 

The motion of the elastic medium is described by mappings from space-time 
onto the material manifold (space-time description) or by time-dependent 
mappings from the material manifold onto space-like hyper- surfaces of space- 
time (material description). The equations of motion are derived from an ac- 
tion principle. In this sense we restrict ourselves (in terms of non-relativistic 
theory) to hyper-elastic materials. In difference to standard non-relativistic 
theory we do not assume the material manifold to be equipped with a flat 
reference metric, since the volume form turns out to be sufficient for all kine- 
matical purposes. An overview on the formulation of relativistic elasticity 
as used in this thesis can be found in For a review emphasizing the 
hyperbolic structure of the equations of motion we recommend j^j. 

While in this thesis we consider elasto-dynamics (i.e. time-dependent 
motions), we note that there exist related work within the same framework 
by Beig and Schmidt on stationary scenarios. They investigate the boundary 
value problem for time-independent formations of elastic matter (like a solid 
star in rigid rotation). For the non-relativistic case we refer to [E], while 
the relativistic problem is treated in [Jj. An application of the formalism to 
static elastic shells can be found in |37j . 

This thesis is organized as follows: 

In chapter 121 we review basic definitions from the theory of hyperbolic 
quasi-linear second order systems of partial differential equations and collect 
local well-posedness results for both bounded and unbounded domains. 

Basic kinematics for relativistic elastic materials in developed in chapter 
|3J For both space-time ("Euler") description and material ("Lagrange") 
description we identify the basic unknowns and derive a covariant formulation 
of relativistic elasticity in terms of these unknowns. For both descriptions the 
Lagrangian is derived in terms of material quantities. The energy-momentum 
tensor is given and a discussion of the reference state (a "realization" of the 
(abstract) material in the actual world) is included. 



9 



Chapter 0] deals with the dynamic properties of elastic materials. For both 
descriptions we derive the equations of motion and prove local well-posedness 
results under various assumptions. We include a study of conserved quanti- 
ties in the presence of space-time symmetries and compute an explicit formula 
for the solution of the initial value problem for the linearized equations of mo- 
tion. Then the self-interaction due to gravity is included and a well-posedness 
result for this setup is proved. Finally we consider the initial-boundary value 
problem and obtain well-posedness for natural boundary conditions. 

A summary of the results of this thesis is given in chapter 0D We finish 
by stating various open problems and possible ways to generalize the results 
obtained in the course of the present work. 

The appendix finally collects remarks on energy conditions satisfied by 
the energy-momentum tensor, non-relativistic limit for the elastic equations 
of motion and 3+1 decomposition for the Chistoffel symbols. It also provides 
us with an easy 1+1 dimensional example to get used to the basic definitions 
as introduced in chapter El 



CHAPTER 1. INTRODUCTION 



Chapter 2 

Hyperbolic Second Order 
Systems 

As soon as one tries to describe a phenomenon in objective terms, mathe- 
matics comes necessarily into play. If the phenomenon is dynamical it is very 
likely that one will not be able to proceed without the tool of partial differen- 
tial equations. If these equations governing this dynamical process arise from 
an action principle (which usually is the case in mathematical physics) these 
will naturally form a quasi-linear (second order) system of partial differen- 
tial equations. We will find various examples for such systems (arising from 
action principles) in the course of this thesis. It is thus favorable to have a 
well-oiled machinery, when it comes to questions concerning well-posedness. 
To provide the necessary tools to prove such results for the systems arising 
in the context of elasticity in chapter 0] is the aim of this chapter: 

We will essentially state two results. The first is concerned with the pure 
initial value problem and is due to Kato et al [301 - T ne second deals with 
the initial boundary value problem and is due to Koch [3^] • We will not use 
these results directly but rather give corollaries, which will turn out to be 
more appropriate for later application. 

2.1 Basics 

We start with a some general definitions concerning quasi-linear second order 
systems of partial differential equations: the latter are of the form 

A%,dpd v f B +B A = (2.1) 

where the unknown f A appears linearly in the highest order derivatives, i.e. 
dfj,d u f B '. The coefficients A A V B and Ba are allowed to depend on lower order 
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derivatives of the unknown and may also depend on x M , the coordinates on 
the base space. One finds that most information about the system can be 
read off from the highest order terms. We thus define 

Definition. The principal part of a system of PDE is the part involving the 
highest order derivatives. 

For the case of equation (|2.1|) the principal part is given by 

45,W B ( 2 - 2 ) 

Definition. The principal symbol corresponding to the principal part \2.2]) 
(and thus to the system A2.1\) ) is given by 

A AB (k) :=A%k,k v (2.3) 

Finally the characteristic polynomial is given by the determinant of the prin- 
cipal symbol 

A(k) := det A AB {k) (2.4) 

Let upper case Latin indices range from one to N. Then the characteristic 
polynomial is homogenous of order 2N in the co-vector k. Note that both the 
principal symbol and the characteristic polynomial depend on the unknown 
/ and its first order derivatives. In this sense the unknown determines the 
character of the equation. This fact has to be kept in mind for the following 
definitions, where the dependence on the unknown will be suppressed. 

In the following we restrict ourselves to a certain type of systems: 

Definition. The system \2.1\) is called weakly hyperbolic w.r.t. a co-vector k, 
if A{k) 7^ and the polynomial A(k + sk) admits only real zeros (si, . . . , S2n) 
for any choice of k, i.e. if the characteristic polynomial is a hyperbolic poly- 
nomial. Such a k is then called sub- characteristic. 

This conditions turns out to be sufficient to obtain well-posedness for lin- 
ear systems (see 02]), but it is too weak to establish well-posedness for the 
quasi-linear equations (|2.1j) . There are various stronger concepts of hyper- 
bolicity, which are not equivalent. Without the attempt to give an overview 
(which can f.e. be found in |3] or rj2]) we state the notion used in this thesis. 

Definition. We call the system \2.1\) regular hyperbolic if it is regular hyper- 
bolic w.r.t. a pair {X^^uJ), i.e. if there exist a pair {X^^uJ) G TM x T*M 
satisfying X^uj^ > and constants (ei, 62) > (0, 0) such that for arbitrary ip A 
we have 

-A% B u^ A *b B > ei\tjj\ 2 (2.5) 
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while 



(2.6) 



for all k^ G T*M satisfying k^X^ = 0. 

Remark. This is an open condition (it involves real inequalities). Suppose 
everything is sufficiently smooth and the system is regular hyperbolic for the 
unknown f . Then the system is also regular hyperbolic in some neighborhood 
of this f (in a suitable topology). 

The second condition in the definition above is of some importance on its 
own. For later need we define: 

Definition. An objects C % l B is called strongly elliptic if the following inequal- 



for a positive real number e. This condition is also called Legendre-Hadamard 
condition or rank one positivity. 

The two notions of hyperbolicity given so far are connected via the fol- 
lowing lemma: 

Lemma 1. Regular hyperbolicity implies weak hyperbolicity. Every w.r.t. 
a pair (X,u) regular hyperbolic operator is also weakly hyperbolic w.r.t the 
co-vector uj. 

A proof of this statement can be found f.e. in [3J. In short the argument 
works like this: Hyperbolicity is defined as a condition on the zeros of a 
determinant. This is equivalent to an eigenvalue problem, which due to the 
symmetry of the involved objects can be shown to admit only real eigenvalues. 

There exists another even weaker notion of hyperbolicity: 

Definition. A co-vector k is called non- characteristic, if A(k) ^ 0. In turn 
a co-vector k, for which A(k) = is called characteristic. 

Clearly a sub-characteristic co-vector is non-characteristic but not the 
contrary has not necessarily to be true. This can be seen directly from the 
definitions. 



ity holds 



Ci B ^\ A \ B > eM 2 |A| 2 




Remark. Note that hyperbolicity in the forms given here is a point-wise (i.e. 
local in particular) condition. 
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To get used to the notation we discuss the scalar wave equation on 
Minkowski space as an easy example: 



The principal part is simply given by the d'Alembert operator acting on 
the unknown: rf v d . Thus the principal symbol is rf ,v k ll k v . Because 
we are only dealing with the scalar case in this example, the characteristic 
polynomial coincides with the principal symbol. We find that light-like co- 
vectors are characteristic and that all other co- vectors are non-characteristic. 
Choosing k = dt we get that the zeros of 



are all real. Since we can always find coordinates such that k = dt at a 
certain point for time-like k, we find that the wave equation is weakly hy- 
perbolic w.r.t. every time-like co-vector. (Thus every time-like co-vector is 
sub-characteristic.) It is easy to see that it is not weakly hyperbolic w.r.t. 
any space-like or light-like vector. Finally we note that the wave equation 
is regular hyperbolic w.r.t. any pair k v ) where and k v are both 
time-like. 

We have mentioned, that the notion of weak hyperbolicity is not strong 
enough to guarantee well-posedness of the Cauchy problem. In the next 
sections we will see that regular hyperbolicity is a sufficient assumption for 
that matter. 

2.2 Local existence with no boundaries present 

The ultimate goal of this section is to give a local existence result based 
purely on the concept of regular hyperbolicity. To this purpose we first state 
a local existence result which is a simplification of a theorem due to Kato et 
al. (SO] and then modify it. The final corollary will in some sense be weaker 
than the original theorem, but will be more natural from a space-time point 
of view, since no preferred foliation is included. For this purpose we make use 
of an domain of dependence result due to Christodoulou [T2]. The original 
theorem by Kato et al. reads as follows: 
Given a quasi-linear second order system 



<" v dM = l.o. 




{k + s) 2 + k 2 



(2.9) 



A%d„dJ B + B A = 



(2.10) 



for an unknown / which is a mapping from [0, T] x M 3 to R 3 . Greek indices 
run from to 3, where x° = t e [0,T]. The coefficients A^ B and Ba depend 
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smoothly on the unknown, its first order derivatives as well as on time t and 
on the point x 1 G M 3 , i.e. 



A% G C b °°([0, T]xR 3 x O) and B A G C fe °°([0, r]xi 3 x tt) (2.11) 



Here Q C M 3 x M 12 is an open set. We assume that (/, df) take values in Q. 

In addition to these differentiability properties we require the principal 
part to show the following symmetry: 



We also assume that the principal part is regular hyperbolic w.r.t. (d t ,dt) 
for some reference state in Q. This means that there exists a real e > such 
that 



for all (t, x, g, dg) G [0, T] x M 3 x Q as well as for all f G T*M 3 and all ^ G TM 3 . 
Under these assumptions the Cauchy problem with initial conditions 



is considered, where we assume that (/o,/i,^i/o) ^- Then the following 
theorem holds: 

Theorem 1. Assume s > 3 and let the preceding assumptions hold. Then 
the Cauchy problem for A2.1U\) is well-posed in the following sense: Given 
(/o,/i) e # S+1 (M 3 ,M 3 ) x # S (M 3 ,M 3 ) suc/i tfitrf (/ ,/i,9i/ ) e 0, ^ere zs a 
neighborhood V of (/o, /i) m x if s and a positive number T* < T such 
that, for any initial conditions in V , \2. 1\) has a unique solution f(t, ■) for 
t G [0,T*] wit/j G fl. Moreover / G C' r ([0, T*], iiP +1_r ) /orO<r<s 

and the map (/o,/i) ^ (f,d t f) is continuous in the topology of H s+1 x IP 
uniformly in t G [0,T*]. 

A proof of a more general version can be found in [30] • 
We will alter this result in two directions. On the one hand we would 
like to generalize it to arbitrary regular hyperbolic systems. This on the 
other hand brings up a new problem, namely the need for a version of the 
theorem, which is local in space. This need arises since regular hyperbolicity 
in general does not lead to a situation covered directly by the theorem. This 
will be seen in the next lemma. Then we will show in a final lemma how 




(2.12) 



(2.13) 
(2.14) 



/(t = 0) = / G^ +1 (R 3 ,M 3 ) 
d t f(t = 0) = f 1 eH s (R 3 ,R 3 ) 



(2.15) 
(2.16) 
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to deal with this problem by giving a way how to continue initial data on 
an open set U C I 3 to initial data on all of 1R 3 in a proper way so that 
the theorem can still be applied. But first we show how the requirement of 
regular hyperbolicity locally leads to a problem of the above type. 

Lemma 2. Assume that the system \2. 10\) obeys the symmetry \2. 12\) and 
is regular hyperbolic w.r.t. a pair (X m ,clv). If the coefficients are regular 
(smooth for example) there exists a coordinate system, every point of space- 
time has a neighborhood where we can find coordinates such that the require- 
ments of the theorem are met. 

In the proof we will not deal with the differentiability requirements but 
will only care about the geometry. First note that the symmetry ()2.12j) is 
unchanged under coordinate transformations. The goal is that of finding a 
foliation, such that to is the co-normal and X is tangent to the time-flow, 
i.e. a coordinate system in which (X,uj) oc (dt,dt). Having this picture in 
mind, we first assume that u> is closed, i.e. that dio = 0. This implies that 
uj = da for a scalar potential a. By choosing this a as new time coordinate, 
we obtain u = dt. For the other coordinates we have to solve 



This means the new coordinate has to be constant along the flow of X. Let 
<fi v s (x^) denote this flow. Then (since X is time-like) the equation <fi s (t,x) = 
(0, z) uniquely determines the parameter s(t). Now we define 



Defined this way, y is constant along the integral curves of X and thus obeys 
the above equation. Hence (a, y l ) constitute a coordinate system satisfying 
the claims of the lemma. 

If u is not closed we have to proceed along a slightly different route. 
Fix an arbirtary point on the manifold we can find a hyper-surface with co- 
normal Co such that at the chosen point Co coincides with ui. Since we assumed 
the system to be regular hyperbolic w.r.t. (X,u), by continuity there exists 
a neighborhood of the chosen point, where the system is regular hyperbolic 
w.r.t. the pair (X, Co). Since Co is hyper-surface-orthogonal by definition, we 
can use the first part of the proof to conclude the result of the lemma. 

Remark. For a regular hyperbolic system we can find data, such that the 
requirements of the theorem are locally met. To apply the theorem, we have 
to extend the initial data to all ofM. 3 . This is the gaol of the following lemma. 



X v d v y l = 



(2.17) 




(2.18) 
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Lemma 3. Let U be a non-empty open subset of M. 3 and let (g , gx) G 
H S+1 (U, M. 3 ) x H S (U, M. 3 ) be small initial data on U, i.e. a pair of functions 
satisfying the requirements for the theorem restricted to the open set U. Let 
U be a non-empty open subset of U. We also assume the existence of a 
reference solution, which induces initial data in the sense of the theorem on 
all ofM. 3 . Then there exists a pair of functions (/o,/i) £ H S+1 {R 3 , R 3 ) x 
H S (R 3 , M. 3 ), which is an initial data set in the precise sense of the theorem 
with the property that (go, g{) = (fo, f\) in U. The development of the initial 
data in the domain of dependence of U is independent of the continuation 
and depends solely on the choice for (go,gi). 

Remark. This lemma allows us to use the theorem also for the case when 
data are only given in a subregion of M. 3 . For the notion of "domain of 
dependence" we refer to jiljj/ . For the given context we might sloppily translate 
the last part of the lemma as: "there exists a unique local solution". 

We denote the initial data on ¥L 3 coming from the reference solution by 
(d ,di). Then on U we define 

(Sgo, Sgi) = (g , gx) - (d , d x ) (2.19) 

One could think of multiplying this with the characteristic function for U 
and adding this then to (do,di). These data would satisfy the algebraic 
conditions, but would not meet the differentiability requirements in general. 
We thus have to smoothen these data. Let x be a smooth nonnegative 
function on all of R 3 such that 

/ 1 on U 

X ~ { OonR m+1 \U ( ° } 

Then x($9o, Sgi) is defined on all of M. 3 and has the required differentiability. 
Thus 

(fo,fi):=(.d ,d 1 ) + x(.8g ,5g 1 ) (2.21) 

has the required differentiability too. We have to check that the hyperbolicity 
conditions are not violated. By assumption they obviously hold in U and 
M 3 \f7. They also hold in the region U\U, since the contribution of x(<^7o, <5<?i) 
can be made arbitrarily small by choosing (5 go, Sgi) small enough (in the 
proper norms, f.e. the supremum norm). This works by the smoothness of 
X- In the region U (fo, fi) and (go,gi) clearly coincide. This proves the first 
part of the lemma. The second part of the lemma (concerning domain of 
dependence) is actually proven in theorem 5.10 of ^2] which we can directly 
apply to the given setup. 
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Corollary 1. Given a regular hyperbolic system of the form \2.1U\j . Then 
the Cauchy problem is locally well-posed in the sense, that for local (in space) 
initial data there exists a unique local (in space and time) solution for the 
initial value problem, which depends continuously on the data. 

The involved function spaces are the same as in the theorem. 

Remark. Note that strong ellipticity of the spatial part (i.e. the Legendre- 
Hadamard condition) was sufficient to guarantee well-posedness. In the pres- 
ence of boundaries the need for more restrictive conditions will arise. 



2.3 Local existence in the presence of bound- 
aries 

In analogy to the previous section we are going to state a local existence result 
on the initial-boundary value problem for quasi-linear second order systems. 
The main theorem is due to Koch [SB] . A similar result by Chrzestzczyk can 
be found in , although the Koch version is much closer to the Lagrangian 
approach used in this work as we will see. Although the theorem in its gen- 
uine form allows for mixed boundary conditions we will only state a version 
involving just one kind of boundary conditions namely such of Neumann- 
type. As before we will first state the main theorem and then discuss some 
of its properties and point out the main differences to the unbounded case. 
Let Q be a bounded domain in IR 3 with smooth boundary dQ. Let '■= 
[0,T)xfl and dVL T : = [0, T) x dtt for < T < oo. Let denote the outer 
co-normal to Qt- Greek indices run from zero to three, while Latin indices 
range from one to three. The independent variables are = (t,x k ). 
Suppose we are given a system of the form 

J-r A = B A (2.22, 

in VLt subjected to the boundary condition 

9> M U T =0 (2.23) 

at 8Qt- Here both 5" and B are functions of x\ the unknown u A and its 
first order derivatives d^u A . Thus the given system is quasi-linear. We seek 
solution for the following initial conditions at {t = 0} x Q: 



(2.24) 
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We assume that (u ,Ui) G H s+1 x H s for s > 3. The principal part is 
obtained from 5F by variation w.r.t. the gradient of the unknown: 

A * B := (2 ' 25) 

Let U be an open neighborhood of graph(wo, u%, diUo) in (QUdQ) xl 3 x R 12 . 
We assume that J and B are smooth functions on U . We additionally assume 
the following symmetry of the principal part: 

Kb = A% A (2.26) 

Note that this relation is automatic, if the system is derived from an action 
principle. Now comes the crucial input, namely the hyperbolicity require- 
ments. From the previous section we expect point-wise positivity conditions 
on certain portions of the principal part. This view will only partly persist. 
For the time- components we require that for any v G U 

-A™ B {v)4> A 4> B >e\<t>\ 2 (2.27) 

for some positive constant e. This is the same condition as we had before in 
the unbounded case. The difference comes with the spatial part: Instead of 
a point-wise condition we now have the following coerciveness condition: for 
any v G C 1 (fi), v\ G C(Q) with graph(f , V\, diVo) G U there exist constants 
e > and K > 0, such that for all smooth cf) 

A%{x, v^, 8^)84%^ > el^H^ -K\\<j>\\l, (2.28) 

This inequality is usually called Garding inequality (see f.e. |Hlj). 

Remark. We will see later in this section, that the Garding inequality implies 
the Legendre-Hadamard condition \2. 7| ). The converse statement however is 
in general not true. 

If k = 0, then this coercivity condition tells us that the H l norm can 
be estimated from above by some bilinear form. If A %3 AB is smooth, then it 
is bounded too and the left had side in ()2.28)1 can itself be estimates from 
above by the H 1 norm. Thus the bilinear form generates an equivalent norm. 
Along this way, one can try to prove existence results. In the general case 
however k^O and this does not work immediately. 

If one is only interested in Dirichlet boundary conditions, the unknown 
can be assumed to vanish at the boundary. In that case the coercivity condi- 
tion is equivalent to the Legendre-Hadamard condition. This will be shown 
in a lemma following the main theorem. 
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Since boundaries are present, a smallness assumption on the initial data 
will not be sufficient anymore. We will have to guarantee that the initial 
data go along with the boundary conditions. This requirement is reflected 
in the so called compatibility conditions, which represent restrictions on the 
initial data. A discussion of the actual meaning of these conditions will be 
given following the theorem. 

Theorem 2. Suppose the compatibility conditions are satisfied up to order 
s. Then under the above assumptions there exists a unique < < T and 
a unique classical solution u G C 2 (Qt, U dQr») of the system \2.22)) - \2.24\j 
with cFuit) G L 2 (Q) ifO < a < s + 1 (d denotes all derivatives). In addition 
the solution depends continuously on the initial data. 

Remark. Note that the assumptions of the theorem involve mainly open 
conditions. Thus if we can show them to hold for some initial data, we 
immediately obtain a neighborhood of the original data as well as admissible 
data. Consequently the existence of a reference solution implies the existence 
of solutions, which are small perturbations. Then the continuity part of the 
theorem tells us, that these solutions remain close to the reference state at 
least for short times. 

As was mentioned in the introduction, a proof of (an extended version 
of) this theorem by Koch can be found [36 . 

We still need to discuss the compatibility conditions. Roughly speaking 
they guarantee, that the initial data satisfy the boundary conditions up to a 
certain order in time. Higher time derivatives of the initial data are expressed 
via the field equations. F.e. we say that the compatibility conditions are 
satisfied to zeroth order, if the boundary conditions hold, where (u, d t u, diU) 
is substituted by (u , u%, diU ). It is quite lengthy to actually write down 
these conditions for general order, so at this stage we only give the recipe 
how they are obtained. 

We introduce (suppressing the x*-dependence) 



n,dm =: r^ST'V, u A n , d t u A , d t u A H , . . . , d^u\ , d? +l u A ) (2.29) 



To express the time derivatives in terms of the initial data we use the dif- 
ferential equations 1)2.22)1 and their time derivatives, d^u for example is 
recovered as a function of the lower time derivatives by solving 



for the highest order in time. This is possible, since the coefficient is given by 
A°a B which is non-degenerate by the assumptions of the theorem. The solu- 
tion is called u^ A . By induction we can eliminate all higher time derivatives 




(2.30) 



2.3. LOCAL EXISTENCE IN THE PRESENCE OF BOUNDARIES 21 



and are left with functions, which are completely determined by the initial 
data and their spatial derivatives. Abusing notation we denote these again 
by u A ( m \ Then the compatibility condition reads 

Definition. We say, the compatibility condition is satisfied up to order s, 
provided 

n,3f^u A , ut 4,i , 4,b u^ A , B , u^ A )\ dn = 
is satisfied for all < k < s . 

Remark. This intuitive definition only works if the ingredients are suffi- 
ciently smooth. 

In general one will have to deal with non-continuous functions, i.e. ele- 
ments of certain Sobolev spaces. In this setup the restriction of a function to 
the boundary is not defined in a strict sense. This problem is overcome by 
introducing the trace operator (see f.e. jTH] or [EI]), which assigns boundary 
values in a proper way. One can show (see ^Hj P- 259) that the trace of a 
function / G H 1 does vanish on the boundary if and only if / G Hq. Using 
this we give an alternative definition for the compatibility condition: 

Definition. We say, the compatibility condition is satisfied up to order s, 
provided 

n^(u A , ut 4,i , utn u^ A , B , u^ A ) G H s - k (Q) x H^Q) 
for all < k < s . 

In general it is not clear whether the compatibility conditions can be 
satisfied. In the actual case we will show, that under certain assumptions a 
large class of initial data can be given satisfying the compatibility condition. 
Assume we are given a time-independent reference solution for the system 
()2.22|) - ()2.23p . We call this solution /. It induces initial data of the form 
(f,d t f)\t=o = (fo,0). Further assume that the principal part satisfies the 
Legendre-Hadamard condition w.r.t. the co-normal n^. 

A(x, f, df r» B n^ A <f > C|0| 2 (2.31) 

for arbitrary and a positive constant C. Then we have the following lemma 

Lemma 4. Assume that A°2 B (x, f, df) = 0. (There are no mixed components 
in the reference deformation.) Then under the two assumptions made above, 
we can derive initial conditions close to those for the reference deformation, 
which obey the compatibility conditions. 
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First of all note that the zero-component of the co-normal vanishes since 
the co-normal is annihilated by dt by virtue of the setup: 

n = (2.32) 

Let us first prove the lemma for the zeroth order. Let d n := n M <9 M denote a 
vector field transversal to the boundary dQ. (If a metric g^ u is given we can 
choose f.e. = g^ u n y ) Then the zeroth order condition reads |an = 0. 

Clearly 

The right hand side is positive definite near the reference solution by the 
Legendre-Hadamard condition. Thus the finite dimensional inverse function 
theorem (see f.e. [JT]) allows us to solve n^ A \an = for d n u . This means 
that the zeroth order compatibility condition is satisfied by a suitable choice 
of a transversal derivative of the initial data at the boundary. 
For the higher orders we first note that 

= n^ B d k t u B + l.o.t. (2.34) 

Again the Legendre-Hadamard condition allows us to use the implicit func- 
tion theorem to solve this for d\d n u. This can be rewritten as d n u^ k \ For 
k = 1 this is nothing but the normal derivative of the initial datum d n Ui. 
For higher orders k > 2 we have to show that d n u^ can be solved for higher 
normal derivatives of the initial data. For this we use the differentiated field 
equation. For k > 2 we can write ()2.3()j) as 

A^ B d^d v d k t - 2 u B = l.o.t. (2.35) 

Using that no = 0, this can be rewritten as 

-A A ° B d n u^ B = n,n v A^ B dlu^ B + l.o.t. (2.36) 

where the terms denoted by l.o.t. may be of the same order but involve A° AB 
and thus cause no trouble as we will see. Note that both A AB and n^n v A^ B 
are non-degenerated in some neighborhood of the reference solution by the 
hyperbolicity conditions. Thus we can iteratively use the above equation to 
determine u^ B in terms of normal derivatives. We have that 

d n u^ B = M B )A d k+l u A + l.o.t. (2.37) 

for even k where M B k ^ A denotes a non- degenerated matrix formed from A° AB 
and n^riyA^ (non-degenerated at least close to the reference solution). Fi- 
nally for odd k we obtain 

8 n u {k)B = Mf k _ 1)A d k u^ A + l.o.t. (2.38) 
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Now we can use the implicit function theorem on both these expressions 
because: For the reference solution this equation is identically satisfied by 
assumption. The l.o.t.-terms are either lower order in the normal derivative 
or do vanish due to the A° AB = assumption. On the other hand Mn,s A 
involves only the initial data and first order normal derivatives. Thus all the 
requirements are met to apply the implicit function theorem. 

We have found that the compatibility condition can in fact be interpreted 
as a condition on the odd-numbered normal derivatives of the initial data at 
the boundary. Thus by a suitable choice of these derivatives, the compati- 
bility conditions can be satisfied. Besides these restrictions the data can be 
given freely, as long as they stay in a certain neighborhood of the reference 
solution. In particular this means that the field, its first order time derivative 
and all of their even normal derivatives are free data. 

Remark. Note that by the above constructions we can choose a large class 
of initial data close to those generated by the reference solution, which satisfy 
the compatibility condition up to arbitrary order. 

In the next lemma we establish the connection between the Legendre- 
Hadamard condition ()2.7j) . which was sufficient for local well-posedness in 
the unbounded case and the Garding inequality ()2.28j) : 

Lemma 5. Let Q be a bounded domain in R 3 . Then the Garding inequal- 
ity for u G Hq(Q) and the Legendre-Hadamard condition \2. 7| ) are 
equivalent. 

First we show, how the Legendre-Hadamard condition implies the Gard- 
ing inequality: For this we will follow Giaqunita [22J. The argument consists 
of three parts: first we assume that A l AB {x) = A %d AB is constant, then we 
derive the inequality for a finite region and in the last step we cover Q by 
such regions. The basic tool used for the first step is Fourier transforma- 
tion. Because Hq(Q) is the closure of C^°(f2) it suffices to prove the result 
for u G C^°(Q). By the Plancherel theorem (see f.e. ^Hj) the L 2 norm of a 
square-integrable function is unchanged under Fourier-transform, thus 

\\d t u A \\ 2 L2 = \\k lU A \\ 2 L , (2.39) 

where an upper hat denotes the Fourier transform. Thus 

A(u, u) = 4vr 2 i^ f kikju^dk^c [ \k\ 2 \ii\ 2 dk = c [ \du\ 2 dx (2.40) 
Jn Jn Jn 

Next suppose that u G C^(fl) such that the support is a subset of U t (xo). 
On U t (x ) we can add and subtract a term involving A 1 ab (xq). Using the 
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inequality above this leads to 

f A i i B {x)a i u A d j u B dx>c [ \du\ 2 dx-u [ \du\ 2 dx (2.41) 
Jn Jn Jn 

The first term comes from the Garding inequality for the constant operator 
A^ B (xo), while the second is obtained by uniform continuity (on U e (x ), the 
difference of A l AB {x) from A 1 ab (xq) must not exceed a;). By choosing e small 
enough we can obtain c — to > 0. 

Since Q is bounded subset of M. 3 it is compact and can hence be covered 
by a finite number of balls B e (xj) C £1 for j = 1, . . . , N. Then we choose 

N 

G C^°(B e (xj)) such that ^ 0| = 1. Then A(u,u) can be written as 




A' 



a^ A a^ B = V / ^^(0^)^(0^) + . . . (2.42) 

The dots stand for terms which are at most linear in derivatives of u. Those 
terms can be handled easily by the triangle inequality using the fact that 
both A %3 AB and (f>i are bounded functions. Thus we obtain bounds from below 
by the L 2 norms of u and du respectively, we obtain 

N 

A(u,u)> c y2 \d((f)iu)\ 2 dx - d\\u\\ L 2\\du\\ L 2 - e\\u\\ L 2 (2.43) 
i=i J n 

As indicated above the constants (especially d and e) depend on (j) and the 
operator A l3 AB . We now may pull (pi out of the bracket in the first term which 
can lead to a change in e and d. We obtain 

A(u,u)>c \du\ 2 - d\\u\\ L 2\\du\\ L 2 - e\\u\\ L 2 (2.44) 
Jn 

Note that e is non-negative in this procedure. We now can apply Young's 
inequality (see f.e. (THj) to obtain the desired result. 

The other direction is much easier to prove. We only have to find an 
appropriate function to insert into the Garding inequality. A good choice is 

<f) A = \ A sm(a{i i x i )xu (2.45) 

where xu denotes an arbitrary smooth function with support in an open 
subset U C ft. Inserting this into the inequality leads to 



a 2 



/ A l AB X A X B jj^i fijxlr cos 2 (a /j,ix l ) > ca 2 / |/i| 2 |A| 2 X[/ cos 2 (a//; a;') + l.o. 
Ju Ju 

(2.46) 



2.3. LOCAL EXISTENCE IN THE PRESENCE OF BOUNDARIES 25 



The terms denoted by l.o. are those of lower order in a. Note that the 
coefficients of the a are all bounded, since we integrate smooth functions 
over bounded regions. We get rid of these lower order terms by a high 
frequency limit. Since all terms are bounded, there exists an ao such that 
for all a > ao the above inequality holds without the lower order terms. We 
have that 

/ cos 2 {a^x l )xl {A i i B X A X B ^ j - c\fi\ 2 \X\ 2 ) > (2.47) 
Ju 

Recall that U was taken to be an arbitrary open subset of Q. This implies 
together with the smoothness of the integrand that the above estimate holds 
point- wise. This in turn allows us to drop the cosine and the function xu- 
The result is the desired Legendre-Hadamard condition. 

Although this section is devoted to boundary value problems, we add the 
following remark: the Lemma can be used to prove the analogous result for 
Q = M 3 : 

Corollary 2. The lemma still holds with the bounded domain fl being re- 
placed by M. 3 . 

To show that the Legendre-Hadamard condition follows from the Garding 
inequality in this more general context we can use the very same argument 
as above without any change at all. The only difficulty arises when we want 
to prove the other direction. The techniques used in the proof of the lemma 
do not carry over, since there we made use of the compactness of Q, which 
now is not given anymore. We have to work along a different route. 

Let therefore <fi G Cq°. Then there exists a compact subset Q such that 
the support of <j) is contained within Q. Hence by lemma El the Garding 
inequality holds on Q. But since (ft has compact support all integrals can be 
taken over all of M 3 . The Garding inequality on all of M 3 is obtained. The 
analogous procedure can be carried out for all 4> G Co°- The desired result 
follows. 

Remark. This lemma tells us that the Legendre-Hadamard condition is still 
sufficient to treat the Dirichlet problem. 

We conclude this section by stating two other positivity notions beside 
the Legendre-Hadamard condition ()2.7|) and discussing their relevance in the 
given context: 

Definition. An objects C l j[ B is called rank-two-positive if the following in- 
equality holds 

Cl B ni^ > eM 2 (2.48) 
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for a positive real number e. Here denotes an appropriate norm, f.e. 
\fi\ 2 := /if/if. We finally call an object Cabcd strong point-wise stable, if 

C ABC Blx AB lf D > e\Lfi\ 2 (2.49) 

for a positive real number e where Lfi AB = \([i AB + H- BA ) denotes the sym- 
metric part of ii. 

Note that by raising or lowering the indices we can write the rank-two pos- 
itivity condition as well as the Legendre-Hadamard condition in the form of 
the strong-point- wise-stability condition. Then we have the following picture: 
The Legendre-Hadamard condition is the weakest of these three, followed by 
strong point-wise stability. The strongest condition is rank-two-posit ivity. 

Note that this rank-two condition is strong enough to guarantee the valid- 
ity of the Garding inequality (|2.28j) . To show this we simply choose \x A = d^ A 
and integrate 1)2.48)1 over a region Q. This immediately gives the Garding 
inequality with e = K > 0. 

We will however show in a later chapter that in fact point-wise-stability 
suffices to obtain the Garding inequality and hence satisfies the requirements 
needed for the theorem. 



Chapter 3 

Geometry and Kinematics 



There have been several accounts to deal with relativistic elasticity, the best 
known maybe that of Carter and Quintana P|. We however will proceed along 
the route of Beig and Schmidt first published (for the space-time setup) in 



In this chapter we will introduce the basic concepts of relativistic elasticity 
due to the latter authors. The basic unknown will be the configuration of 
the elastic material, which is described as a mapping between manifolds. We 
will work in two different settings. On the one hand we will use the space- 
time-description (also known as Euler picture) where we assign to each point 
of space-time a particle. On the other hand we will work with the spatial 
description (or Lagrange picture), where for every instant of time the spatial 
position of the particle is determined. We will write down the Lagrangian 
for both cases. Later in chapter 14.11 the field equations will be derived from 
these Lagrangians. 

As we will see it is more natural to introduce the main objects of the 
theory in the space-time picture. This is due to the fact that the material 
description requires a foliation of space-time, therefore lacking the obvious 
covariance properties of the space-time description. 

There are two manifolds naturally involved in the formulation of this 
theory. One of them is space-time Ai, which is a four-dimensional manifold 
equipped with a Lorentzian metric g^ u . We assume that A4 is time-orientable. 
Let (/x = 0, . . . , 3) denote the coordinates on M.. 

The other one is called material manifold or simply body. We assign it 
the letter B. It is a three-dimensional manifold equipped with a volume-form 
Vabc- The coordinates on B are X A {A = 1,2,3). 

Other manifolds, which appear in this work are constructed out of these 
two. 

The material manifold can be viewed as an abstract collection of the 
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particles making up the continuous material. (It is the idea of the material, 
as Platon would say.) Usually people assume B to be a Riemannian manifold, 
but we find that the existence of a volume form (which is of course a weaker 
assumption) is sufficient. In standard non-relativistic elasticity the metric 
on B is used to define the natural (reference) deformation. In our approach 
we go the other way round: the natural deformation (possibly) defines the 
metric on B. This is more satisfactory, since one would not expect the natural 
deformation to have the same body metric for different space-time metrics, 
which would be a consequence of the standard approach. 

3.1 Space-time description 

In this section we treat elasticity in the space-time description. The un- 
knowns are a set of scalar fields on space-time taking values on the body. We 
will define the basic quantities and the energy-momentum tensor as well as 
the Lagrangian are derived. 

3.1.1 Basic setup 

The basic unknown in this setup is f A , a mapping from space-time onto the 
body: 

f:M - B (3.1) 
i-> f A (x) (3.2) 

Assumption. We place the following restrictions on f A : 

• d fl f A as a mapping from TM — > TB has maximal rank. 

• the null space of d tli f A is time-like w.r.t. the space-time metric. 

Definition. We call f A configuration. d ii f A will be called configuration gra- 
dient. 

Given a configuration f A . Then the conditions on the configuration gradi- 
ent uniquely determine a time-like direction in the tangent space TM.. Thus 
there exists a unique time-like future-pointing vector-field u 11 G TM with 
the properties 



u^f A = 
g^u^u" = -1 



(3.3) 
(3.4) 
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Definition. The u 11 as defined above is called four-velocity. 

Remark. It is important to note that in this context -u M is not a given vector- 
field, but rather a functional of the configuration. 

Next we look at the pullback of Vabc, the volume form on B, along the 
configuration, i.e. we consider 

v^x:= f A ^f B „f,xV A BC (3.5) 

which by definition is orthogonal to u^. Let e pv \ p be the volume-form coming 
from the space-time metric. Then 

e» uXp v vXp (3.6) 

is time-like and non-zero. Without loss of generality we assume that it is 
future-pointing. One easily finds that is orthogonal to d p f A and consequently 
proportional to u^. Thus there exists a positive quantity n, such that 

nu" := li^vuxp (3.7) 
Definition. We call n particle number density. 

Remark. As was the four-velocity, the particle number density is a func- 
tional of the configuration involving the configuration gradient. 

The definition is justified by the following observation: the vector-field 
nu^ satisfies a continuity equation: 

V M (W) = (3.8) 

This can be verified using ()3.7j) : 

V p (e^v uXp ) = {V^)v uXp + e^{V p v vXp ) (3.9) 

The first term vanishes, since the volume form on space-time is covariant 
constant. The second is zero, since V is closed and whence also v is closed. 

Next we define the strain. Similar to a metric (as the name suggests) 
it is a measure for the strain arising in the material as a consequence of 
configuration. 

Definition. The strain h AB corresponding to the configuration is defined as 

h AB :=d p f A g^d u f B . (3.10) 
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Note that h AB does not in general give rise to a rank two tensor on B, 
since it depends on the coordinates on space-time. 

Lemma 6. A necessary and sufficient condition for the strain to define a 
tensor on B is that the four-velocity is Born-rigid, i.e. that the Lie- 
derivative of metric projection orthogonal to -u M in the direction vanishes: 

£u{gnv + UnUu) = (3.11) 

The strain as defined in ()3.10j) is no tensor, since / is no diffeomorphism. 
But if the strain is constant along the flow of the four-velocity, then the strain 
is uniquely defined on the quotient (of space-time along the four- velocity). 
But the latter clearly is diffeomorphic to the body if the four-velocity is 
geodesic (one just has to think in terms of initial data). This is the case for 
solutions of the field equations as we will see. Thus the strain can be viewed 
as a tensor field. If on the other hand the strain really comes from a tensor 
field on B, then it must be constant along the flow of u. Both arguments 
come down to 

C u h AB = (3.12) 

Since the Lie-derivative of the configuration gradient (viewed as a collection 
of one-forms on space-time) vanishes according to 

C u d,f A = d^d u f A ) + u v d^d v] f A = (3.13) 

we find that this is equivalent to the statement that the Lie-derivative of the 
inverse metric lies in the kernel of the configuration gradient: 

= C u h AB = f A „f B , u C u g^ (3.14) 

Consequently there exists a scalar A and a vector v such that 

C u g^ = -Au»u v - 2v { »u u) (3.15) 

where without restriction of generality we can assume v to be orthogonal to 
u. (Otherwise just split v into a parallel and an orthogonal part and re-define 
A and v in the obvious way.) We can lower the indices: 

ZuQuu = + 2v( J1 u v ) (3-16) 

First note that A = according to 

A = u»u u £ u g^ = C u (g^uV) - 2u li C u u" = (3.17) 
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since u 2 = — 1 and C u u = for every vector field u. To compute v we 
contract ()3.16|) with u: 

-v,j, = u h '£u9\w = (3.18) 

Thus we find that 

£>v3iw = -2C u {u( M )u y ) (3.19) 
which is equivalent to the statement of the lemma. 

Remark. Note that the Born-rigidity condition constitutes a restriction on 
the space-time geometry, since the existence of a Born-rigid vector field is not 
trivially given. An example, where a Born rigid four-velocity may exist, is 
that of a space-time admitting a time-like Killing vector field. The existence 
of a homothetic vector field however turns out to be too weak. 

To see the last point of the remark, we insert 

£«fiV = c 9n* (3.20) 

into (|3.11j) . For nonzero C we obtain 

9nu = -2u(jj,u v ) (3.21) 

which contradicts the non- degeneracy of the metric. Note that Born-rigidity 
does in turn not imply homotheticity, since contracting ()3.11|) twice with 
f A ,^ leads to 

f AtM f Bv C u g^ = (3.22) 

independent of the values of A and B. Thus C u g ^ Cg for non-zero C, since 
this would again lead to a contradiction. 

Leaving behind the question of functional dependence and Born rigidity 
we note some further observations regarding the strain: 

• the strain is positive definite 

• therefore there exists an inverse, denoted by h^B- 

• since the four-velocity is annihilated by the configuration gradient, we 
have that h AB = f A ,^ f B , v (g^ u + su^u v ) for arbitrary real s. 

The first statement comes from the fact, that the configuration gradient in 
the definition of the strain projects out the (time like) -u^-direction, leaving 
behind the positive definite rest. The other points follow immediately. 
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We can use the strain to rewrite the particle number density as the pos- 
itive root of the determinant of the strain, taken w.r.t. the volume form on 
the body: 

6n 2 = det(h AB ) = V ABC V AIB , c >h AA 'h BB 'h cc ' (3.23) 

This can be seen by contracting ()3.7|) with itself and using the standard 
formulas (see f.e. jlH] p. 83) to evaluate the e-tensors. 



3.1.2 The space- time Lagrangian 

Our field equations will be the Euler-Lagrange equations associated with an 
action principle connected with the action 

&[/]:=/ p(f,df;g)y^tetg~d 4 x (3.24) 
Jm 

The Lagrangian density p is assumed to be a non-negative function of the 
configuration, of the configuration gradient and of the metric. An explicit 
x^ dependence only enters via the space-time metric. We subject p to the 
following assumptions: 

Assumption. The Lagrangian is covariant under space-time diffeomorphisms, 
i. e. it transforms as a scalar for given configurations. 

Remark. In non-relativistic elasticity this condition is called "material frame 
indifference" (see e.g. or 

An equivalent formulation (see f.e. [31]) of this assumption is 

= + ^AfV +^dT (3-25) 

for arbitrary vector fields £ M G TM.. If we write the Lie-derivatives in terms 
of partial derivatives we obtain 



By contraction with this implies 

V = (3.27) 



dg^ x 

This by the symmetry of the expression in turn implies 



<)P A AB f A ^f B w (3.28) 



dg 
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for some Aab- From this the following is clear: 

P(f \ f\ ; an = p(f A , fV ; <T + suW) (3.29) 

for arbitrary real s. Let s = 1. We claim that 

g t ,v + u ll u v = f A ^f B , u h AB (3.30) 

The objects on both sides are symmetric. We verify the claim by contracting 
with and f c , a g afJ- respectively. For both cases it is easy to verify the 
claim. Since those form an orthonormal basis of TAi, the claim is proved. 

The above formula then tells us that the Lagrangian can be rewritten as 
a function of the configuration, the configuration gradient and the strain: 

p(f,df;g) = a(f,df;h) (3.31) 

Plugging this into equation ()3.26j) we see that 

da 



df 



_ _ A 

A 



fV=0 (3.32) 



Since the configuration gradient is assumed to be of maximal rank, the first 
term has to vanish. Thus a is in fact explicitly independent of the configu- 
ration gradient (Of course it still enters implicitly via the strain). Keeping 
this in mind, we spit the Lagrangian p into 

p = ne (3.33) 

where by definition e is a positive function depending solely on the configu- 
ration and the strain 

e = e(f,h) (3.34) 

This factorization is justified, since the particle number density n can be 
written as the determinant of the strain (see ()3.23)l ). thus involves only the 
strain and (via the volume form on the body) the configuration. The action 
()3.24|) takes its final form 

S s [f}:= [ m(f,h)^fetg'd 4 x (3.35) 
Jm 

where the Lagrangian 

£ s = ne(f, h)y/-detg (3.36) 



depends solely on points of the base space (via the metric), on the strain and 
on the configuration. 



34 



CHAPTER 3. GEOMETRY AND KINEMATICS 



Definition. We call e stored energy function or simply stored energy. 

The stored energy function contains all the information about the mate- 
rial under consideration. Knowing it is equivalent to knowing the equation 
of state. 

We can derive some important quantities from the stored energy. First 
of all it provides us with an intrinsic definition of stress (depending only on 
the material): 

Definition. The second Piola-Kirchhoff stress tensor is defined as the first 
variation of the stored energy w.r.t. the strain: 



T ^ ■= ~^AB ( 3 - 37 ) 



As was the case for the strain, the Piola-Kirchhoff stress is in general no 
tensor on the body, since it depends on the space-time coordinates x^. We 
have the obvious symmetry property 

tab = t B a (3.38) 

We will mainly use this second Piola-Kirchhoff stress tensor and whenever we 
talk about the Piola stress we mean exactly this one unless said otherwise. 
For completeness we stress that one can define a first Piola-Kirchhoff tensor 
too. Let us denote it by a s . In the space-time description it is defined by 

de 

One can check that the first and second Piola-Kirchhoff stress tensors are 
connected via 

r)h CD 

<a = tcdZa- = ^cnSff^g^ (3.40) 



<a ■= 1T7^ (3-39) 



Note that c^u^ = 0. 



3.1.3 Elasticity operator 

The second variation will turn out to be the most important quantity for 
questions of local well-posedness in the neighborhood of stress-free deforma- 
tions: 

Definition. The second variation of the stored energy w.r.t. the strain is 
called elasticity operator. It is defined as 

d 2 e 

Uabcd ■= dhABdhCD (3-41) 
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Again this does not in general define a tensor field on the body. From 
the definition the following symmetries are evident: 

Uabcd = Ubacd = Uabdc — Ucdab (3.42) 

Remark. Note that the second Piola stress and the elasticity operator are 
introduced without referring directly to the configuration. This fact will allow 
us to use the very same definitions in the material setup. 

Before giving an important example for the elasticity tensor, we give a 
short definition concerning symmetries of the material. 

Definition. A material is called homogeneous, if the stored energy e depends 
solely on the strain and does no more include the configuration itself. This 
is equivalent to 

e = e{h) (3.43) 
We also require the volume form V to be independent of X. 

Assume there is given a metric Gab on B (f.e. coming from the strain of 
a special Born-rigid reference configuration). Further assume that the stored 
energy is covariant under diffeomorphisms of B. Then it depends solely on 
the invariants of the strain taken w.r.t. G. We cast this into the following 
definition: 

Definition. We call a material isotropic if the stored energy e depends on 
the configuration gradient only via its invariants, i.e. 

e = e(f,trh,trh 2 ,deth) (3.44) 

An important special case of such materials are the so called St.Vernard- 
Kirchhoff materials (see f.e. [Hj). They are homogenous and isotropic and 
their elasticity tensor is given by (remember that G denotes some metric on 
B) 

Uabcd = 2/xG a(cG 'd)b + ^GabGcd (3.45) 

where fx and A are constants. 

One can formally develop the stored energy for an isotropic (not neces- 
sarily homogeneous) material into a Taylor series around the metric G and 
find that the quadratic portion can always be written in the form (G AB is 
the inverse metric) 

(2fxG A{c G D)B + XGabGcd) {h AB - G AB ) {h CD - G CD ) (3.46) 

for fx and A being functions on B. If the material is homogeneous, then they 
are constant. 
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Definition. The constants fi and A appearing in the quadratic portion of the 
Taylor series of the stored energy for a homogenous isotropic material are 
called the Lame constants. 

This in particular means that the constants appearing in (|3.45J) are the 
Lame constants. For the special case of St.Vernant-Kirchhoff materials they 
completely determine the material. 

To get used to the notation we apply the three positivity notions intro- 
duced in section l2~T1 to the St.Vernard-Kirchhoff elasticity operator: What 
one obtains are conditions on the Lame constants: 

For the Legendre-Hadamard condition we have to contract ()3.45jl with 
a A b B a c b D . Then the condition reads 

(/i + A)(a • b) 2 + fia 2 b 2 > ea 2 b 2 (3.47) 

By splitting b in a part parallel and a part orthogonal to a, we find that the 
above condition holds for all choices of a and b, if the Lame constants satisfy 
the inequalities 

H>0 2fi + X>0 (3.48) 

These combinations, as we will see in section 14.1.11 are crucial entities for 
the hyperbolicity of the corresponding equations of motion. We will also see, 
that these conditions (|3.48J) are equivalent to the statement that plane waves 
in the linearized theory propagate at real (and not imaginary) speed. 

Next we do this computation for strong point- wise stability (12.49)) : The 
inequality, we have to satisfy for arbitrary m AB is (recall that all finite di- 
mensional norms are equivalent. See f.e. j2]) 

(2fiG A{c G D)B + XG AB G C D)m AB m CD > 2tG A{c G D)B m AB m CD (3.49) 

Evaluating both sides leads to 

(fi - e)n^n {AB) > (3.50) 
(2/x + 3A - e)m A > (3.51) 

where n denotes the trace-free portion of m. The above conditions are sat- 
isfied for all possible choices of m AB if 

(U>0 2/i + 3A>0 (3.52) 

We see that there is a tighter restriction on the possible values of A. 

Finally we come to rank-two-positivity. The corresponding inequality 
reads 

(2fiG A{c G D)B + XG AB G C D)m AB m CD > eG AC G DB m AB m CD (3.53) 
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At the first glance it may look similar to the case of strong-point-wise- 
stability, but this is not true, since on the right hand side antisymmetric 
terms do appear. But sine on the left hand side only symmetric terms are 
involved, the antisymmetric terms on the right hand side cannot be bounded. 
Thus there is no possible way to satisfy rank-two-positivity of the elasticity 
operator for St.Venant-Kirchhoff materials (This holds for general isotropic 
materials. For this check the following remarks.) Note however that it makes 
sense to apply the notion of rank-two positivity to certain combinations in- 
volving the elasticity operator. 

Remark. We stress that although we did the above example for the case of 
St. Vernant-Kirchhoff materials, the analogous treatment works for arbitrary 
isotropic materials. 

Remark. Note that G entered in the conditions. It thus may well happen 
that the conditions are met for a certain choice of G and cease to hold for a 
different choice. We have also seen, that some positivity conditions may just 
fail for certain types of materials. 



3.1.4 Energy momentum 

Next we derive the energy-momentum tensor for this model. We can derive 
it from the Lagrangian according to the formula (see f.e. [HI] p. 42 ) 



T^:=-^-f A ,u-pS^ (3.54) 



Recall that the Lagrangian p did only depend on the configuration and the 
strain. This allows us to write (recall the definition of the strain (|3.10jl ) 

dp A _ dp dh CD A _ 2 J^ fC f D a, (3 55) 



On the other hand 

dp dp dh CD dp c D 



Q gl xv Q h CD Qgt+v Q h CD 



r„r,u (3.56) 



Beside a factor 2 this is just the same term. Thus we can rewrite the energy- 
momentum tensor as 

T,u = 2^-pg, u (3.57) 

Note that this expression is manifestly symmetric. Contracting with the 
four-velocity and using ()3.27|) we find that 

T^u" = -pUf, (3.58) 
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Thus the energy-momentum tensor must be of the form 

Tyj, = pu^ + Sfu, (3.59) 

where S is symmetric and orthogonal to u, i.e. S^ u u u = 0. (S is sometimes 
called Cauchy-stress. ) This in turn implies that there exists an object Aab 
such that S/a, = / A , M f B , v A AB . 

Remark. Note that this representation of the energy-momentum tensor al- 
lows a good interpretation for the p entering the Lagrangian. It plays the role 
of the total energy density of the material in its rest system. 

To compute the actual form of A, we have to use p = ne to write out the 
right hand side in ()3.57|) : 



nr A B + QfoAB~ e ) f A '» f B iv ~P9^ (3-60) 

To compute the variation of the particle number density, we use a well known 
result that connects the variation of the determinant with the inverse tensor, 
namely (note that the strain is symmetric and that h A B is the inverse strain 
as defined in section 15.1.1)) 

On 1 

-nfiAB (3.61) 



dh AB 2 

The expression entering the energy momentum tensor is thus given by the 
symmetric tensor h^ u := hABf A ,^f B ,u- A closer look shows that 

V« M = (3.62) 

as well as 

KJ M " = h AB f\h BM = f M ,„ (3.63) 

While the four-velocity is annihilated, h^ u acts as identity on the configu- 
ration gradient and hence also on the orthogonal space. We conclude the 
formula 

?V = 9i*v + u v u v (3.64) 
Insertion of this into the expression for the energy-momentum tensor yields 

= pu^ + 2nr AB f A ^ f B ,u (3.65) 

This justifies to call the second Piola-Kirchhoff tensor tab stress, since it is 
equivalent to the Cauchy stress. Note that there are no currents involved in 
the energy-momentum tensor. 
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Remark. The energy-momentum tensor is divergence- free, if the configura- 
tion is stress-free and the corresponding four-velocity is geodesic. 

This can be shown by direct computation: 

V V T^ = V u (nu»)eu» + nv?V v {eu») + 2<T / V>WV f B » ) (3-66) 

The last term is zero, since the configuration is assumed to be stress-free. 
The first term vanishes due to the continuity equation ()3.8|) . The second 
term can be split into three parts: 

nu v (eV v u» + u^f A ,„ +u»r AB h AB , v\ (3.67) 

All these terms vanish. The last because the configuration is assumed stress- 
free, the second because the four-velocity is orthogonal to the configuration 
gradient and the first because the four-velocity is geodesic. 

Remark. From the proof of the last remark we can extract another interest- 
ing observation: The divergence of the energy-momentum tensor (this time 
for a general configuration) is annihilated by the material velocity. Thus 
V V T^ V = are practically only three equations, which are equivalent to 



This can be seen by contracting ()3.66jl with the four-velocity and then 
evaluating the resulting expression. We will however proceed along a shorter 
route. Since by the Leibnitz rule 



u 



t V v T^ = V„«0 - (3.68) 



we find that (using that the four-velocity annihilates the configuration gra- 
dient and is normalized) 

u^ v T^ v = -V v {pu v )2nr AB f A ^ f B , u W (3.69) 

We can split the total energy density according to p = ne. Then the first 
term becomes 

-V„(piO = -eV v {nu v ) - nu u V u e (3.70) 

The first term vanished due to the continuity equation which is identically 
satisfied in our theory. The second term gives 

Be 

-nu v V v e = -nu v ^-f A , v -nu»T AB V u h AB (3.71) 
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The first is zero because the four-velocity annihilates the configuration gra- 
dient, the second term prevails. We finally obtain 

u^ v T^ = -nu u r AB V v h AB - 2nr AB f A , fl f lV VV (3.72) 

This is altogether zero. To see this we first substitute 

u v V v h AB = 2g^u»f A , a VJ B \ P (3.73) 

Since the configuration is a collection of scalars, we can use that the covariant 
derivative is torsion-free. This leads to 

u v V v h AB = 2g a ^f( A , a V f3 f B \ u =-2f( A , a f B \ u V a u» (3.74) 

Thus the two terms in ()3.72|) cancel and the statement of the remark is shown 
to be true. 

3.1.5 Reference state 

In the definition of isotropy we needed to assume the existence of a metric G 
on the body. We now assume that G is generated by some special reference 
configuration: 

Definition. We call a configuration f A reference configuration whenever it 
is stress-free, i.e. 

T AB\ f= f = (3.75) 

and the corresponding four-velocity is geodesic. Quantities evaluated at this 
reference configuration will be denoted by an upper tilde (like e.g. f). Rela- 
tions, which hold only for this reference configuration will be denoted by an 
o above the relation (f.e. =). We sometimes refer to the reference configura- 
tion more generally as the reference state of the material (this notation will 
be used in the material setup too). 

Remark. There is another possible way of defining the reference configu- 
ration: one may assume it to be stress-free and in addition to generate a 
divergence-free energy-momentum tensor. Then the four-velocity is geodesic 
if the reference deformation solves the field equations. This can be seen from 
the following arguments. We shall however stick to our original definition. 

By the last remark of the foregoing section we know that the energy- 
momentum tensor for the reference state is divergence- free. We will see in 
section Bini that this is equivalent to the statement that the reference config- 
uration solves the field equations (which will be the Euler-Lagrange equations 
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corresponding to the Lagrangian ()3.36|) .). But we would like to show that 
also slight perturbations of this reference state give rise to solutions. This 
will be shown to be the case (at least for the unbounded setup), if the ref- 
erence configuration satisfies a certain minimality property, namely that the 
elasticity operator satisfies the Legendre-Hadamard condition: 

Definition. We call a reference configuration natural, whenever the elastic- 
ity tensor for this reference state satisfies the Legendre-Hadamard condition 
\2. 7| ). In formulas this reads 



for all a and b and for a positive e. 

For most of the following work we will assume the existence of such a 
natural configuration. 

Remark. These definitions did not claim Born-rigidity for the reference 
state's four-velocity. Thus the reference strain need not define a tensor field 
on B. This is a difference between the formulation given here and the stan- 
dard (non-relativistic) formulations of elasticity, where a (flat) metric on the 
body is assumed throughout. 

We have seen in section 13.1.11 that the existence of a Born-rigid vector 
field constitutes a restriction on the space-time geometry. We will cast this 
statement into a more concrete form in the following lemma: 

Lemma 7. // the reference strain defines a tensor field on the body, then 
space-time has to be stationary. 

We have shown in lemma |B1 that the strain defines a tensor if and only 
if the Born rigidity condition holds. Thus the problem is the following: we 
have to show that 



for a four- velocity satisfying « 2 = -1 as well as u u 'V u u fl = and the Born 
condition ()3.11|) . The latter implies the Killing equation whenever 



U ABCD a A a c b B b D > e\a\ 2 \b\ 2 



(3.76) 



(3.77) 








(3.78) 



This is equivalent to 




(3.79) 
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Both terms on the right hand side vanish due to the assumptions of the 
lemma. This in turn implies the Killing equation to be valid for the reference 
four- velocity. Consequently space-time has to admit a time-like Killing vector 
field since otherwise one would have a contradiction. Thus the metric has to 
be stationary at least. 

Remark. This lemma implies that a non-stationary space-time can not ad- 
mit a Born-rigid reference four-velocity. In such a setup the reference strain 
can never be a tensor on the body. 

We end this section with a slight adjustment of the original definition of 
isotropy: 

Definition. We call a material isotopic (w.r.t. the reference configuration) 
whenever it is isotropic w.r.t. the reference strain. 

For convenience we slightly change notation in the following way: if the 
material is isotropic, we redefine the Lame moduli introduced in (I3.45j) ac- 
cording to 



This definition of isotropy will be used in the following. 

3.2 Material description 

In this section we introduce the material description. The dependent vari- 
ables in this setup are time-dependent mappings from the body onto space- 
like hyper- surfaces of space-time. Roughly speaking they are the instanta- 
neous inverse to the configurations. In this section we will perform a 3 + 1 
decomposition of space-time, translate all main objects into this new descrip- 
tion and finally rewrite the Lagrangian for this setup. 

3.2.1 Basic setup 

Assume there exists a foliation of space-time by space-like hyper-surfaces. 
This means space-time is the product of a three-dimensional space-like man- 
ifold S with an interval /cR: 



Let the leafs of the foliation be the level sets of the function t. Introducing 
coordinates x % (i = 1,2,3) on S, we can write the space-time metric (using 
t =: x° as fourth coordinate) as 




(3.80) 



M = I x E 




Q ' fjii/ dx dx 



N 2 dt 2 + g^dx* + Y i dt)(dx j + Y j dt) 



(3.82) 
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The inverse metric reads 

gTd^ = -N-\d t - Y'dA 2 + g^didj (3.83) 

All functions depend on (t, x l ). Clearly we can write every metric in the given 
form. This is the ADM-representation of the metric (see ^H] for details). The 
function N is assumed to be positive and is called lapse while Y % is called shift. 
The geometrical meaning of the quantities entering the ADM-representation 
is the following: g^(t,x) is the induced metric on the leaf H t := {t} x S of 
the foliation, is its inverse, i.e. g^gjk = 5^. The norm of dt is given by 
—N~ 2 . Thus the unit co- normal is n^dx^ = ±Ndt. The (future-pointing) 
unit-normal is n^d^ = N~ 1 (d t — Y % di). 

Remark. One has to be careful about the notation: although g^ denote the 
spatial components of g^ u dx^dx u , the analogous statement fails to be true for 
the inverse metric. There the spatial components involve both lapse and shift 
as well as g lJ . 

Given a configuration f A . Since the corresponding material velocity 
is assumed to be time-like, we have that n^ L ^ 0. This means that 
is transversal to the foliation. This implies that for each instant of time t, 
the configuration is a diffeomorphism between the hyper-surface S t and the 
body, since the configuration gradient is assumed to be of maximal rank, 
thus generating a linear isomorphism between TS t and TB. 

It thus makes sense to write down the inverse of the configuration for 
every instant of time. We define 

Definition. The deformation F % corresponding to the configuration f A is 
defined by (remember that the coordinates on the body B were denoted by 
X A ) 

f A (t,F\t,X))=X A (3.84) 

We note that this definition involves the chosen foliation. We also note 
that by the above statements the deformation is well-defined, since we can 
solve ()3.84|) for it by means of the implicit function theorem. It is a collection 
of time dependent scalar maps from the body onto the leafs of the foliation. 

F:IxB -> {-JxEcM 
(t,X A ) » F%X) 

Remark. One could easily construct a diffeomorphism between space-time 
and I x B by extending the configuration and the deformation into four- 
dimensional mappings of the form (t, f) and (t, F). We will however (mostly) 
stick to our original formulation. 
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Definition. In accordance with the notion of configuration gradient, we call 
F 1 ,a deformation gradient. Note that (in contrary to the configuration gra- 
dient) the deformation gradient defined this way is a linear isomorphism. 

From the definition of the deformation f)3.84j) we can derive the formal 
connection between the configuration gradient and the deformation gradient: 
differentiating w.r.t. X B leads to 

f A ji (tj F)F i ,B — o~b (3-85) 
Contracting with F\a leads to 

F\ A f A n {t ) F)F\ B =F\ B (3.86) 

which in turn gives 

F^ A f A n (t,F) = S{ (3.87) 

We see that the configuration gradient can be obtained as the inverse of the 
deformation gradient and vice versa. 

Another consequence of the definition of F l is obtained by differentiating 
fl3.84|) with respect to t, namely (an upper dot denotes differentiation w.r.t. 

*) 

f A (t, F(t, X)) + f A n (t, F(t, X))F% X) = 0. (3.88) 

This equation connects the time derivatives of the configuration and the 
deformation. 

Remark. From now on we will drop the arguments of the configuration gra- 
dient. Whenever we write it down we understand it as the inverse of the 
deformation gradient by means of \3. Sty . \3.81\j and \3. Sty . 

For our purpose it turns out to be more convenient to introduce a different 
normalization for the four-velocity which is in accordance with the foliation. 
We introduce a new vector field v M , which is parallel to the four-velocity and 
satisfies v fl (dt) fl = 1. Thus 

vPdp = d t + v l di (3.89) 

Since the four- velocity is annihilated by the configuration gradient, the same 
is true for v: 

= f A + v l f A n (3.90) 
Subtracting this from equation ()3.88|) we obtain 

v% F(t, X))f A n (t, F{t, X)) = f A n (t, F(t, X))F% X) (3.91) 
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Since ff is non-degenerated, this implies 

v%F(t,X)) = F i (t,X) (3.92) 

Thus 

v"d li = d t + F i (t,f(t,x))d i (3.93) 

Definition. From now on we use the term four-velocity for t> M in the context 
of the material setup. If we refer to -u M we will speak of the normalized four- 
velocity. 

The norm g^v^v" of the four-velocity is given by 

v 2 = -N 2 + gij (Y l + F i )(Y j + F j ) (3.94) 

Introducing 

W*:= — + (3.95) 

we find that the four-velocity is time-like if and only if 

1 - W l gij W j =: 1 - W 2 > (3.96) 

Assumption. Throughout the rest of this thesis we will assume that this 
condition holds. 

This equation may be understood as a restriction on the possible choices 
for the time flow. Roughly speaking lapse and shift have to be chosen so that 
they balance the spatial evolution of the material. 

This quantity W l will appear throughout the following analysis. We will 
also frequently encounter it in the following context: 

7 :=(1-V^ 2 )- 1 (3.97) 

The condition ()3.9(i|l is then equivalent to 

1 < 7 < oo (3.98) 

We will find that this factor plays a similar role to the 7 factor of special 
relativity. (This can be seen clearly in the Newtonian limit. We refer to the 
appendix for details.) 

Next we deal with the strain. We can split the configuration gradient into 
its spatial part and the rest. Then the strain becomes 

h AB = f\ f B , 3 (g ij - ^Y l Y 3 ) - -^f A f B + -^f (A ,i f B) Y l (3.99) 
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Recall that the first term inside the bracket is the spatial part of the inverse 
metric (|3.83|) . Then we use ()3.88|) to replace the time derivatives of the 
configuration by those of the deformation. One obtains 

h AB = f A n f B ,j {g ij - WW) (3.100) 

If one understands the configuration gradient as being the inverse of the 
deformation gradient, then this is an expression for the strain only involving 
the deformation and the deformation gradient besides the metric. We have 
thus successfully translated the strain into the material description. 

Remark. Abusing notation we keep the letter h for the strain both in the 
space-time and in the material description, although strictly speaking those 
are different objects. The same remark applies to most of the following cases. 
But since it is clear from the context, which description is used, we stick to 
this simplified notation. 

The inverse strain is given by 

h AB = F a , A F\ B g ia9jb (g ij + -yUTW*) (3.101) 

To prove this first note that h A B as defined is clearly symmetric and positive 
definite, sice 7 is bounded from below. 

To show that it is indeed the inverse we contract with the strain: 

h AB h BC = f A , a f B , b (g ab - W a W b )F c , c F\ B g CJ gu{g l] + (3.102) 

Now we use ()3.85j) and ()3.87j) to evaluate the contraction of the configuration 
gradient with the deformation gradient. The result of this is 

h AB h BC = 5 A - f\ W a Wig l3 F\ c (1 - 7 + W 2 ) (3.103) 

The term in the bracket vanishes due to the definition of 7, which proves 
that h A B is indeed the expression for the inverse strain in the material setup. 

3.2.2 Material Lagrangian 

To translate the space-time action ()3.35|) into the material description, we 
have to find expressions for 

• the particle number density ()3.23|) 

• and the stored energy ()3.34|) . 



• We also have to transform the volume element in the proper way. 
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We start with the stored energy. It is a scalar build from the strain and 
the configuration, which themselves can be treated as a collection of scalars. 
Thus for the material description we have 

e = e(X,h) (3.104) 

It is an advantage of the material description that the dependence of the 
stored energy becomes more clear: the energy stored in the material depends 
on the strain acting at each point. Different points of the (non- homogeneous) 
material react differently to strain. 

Remark. From this expression for the internal energy it is clear, that we 
can simply use the definition given for the second Piola stress \3.31\j and the 
elasticity operator \3.41\) - We just have to replace the configuration f A by 
the image X A . 

Only the first Piola tensor causes trouble, since it involves a variation 
w.r.t. the configuration gradient. For the material setup we introduce a 
new definition via the variation of the internal energy w.r.t. the deformation 
gradient: 

<■ : = inh < 3 - 105 ' 

But again this definition is given just for reasons of completeness. 

Coming back to the original problem, namely that of translating the 
Lagrangian, we next show how to deal with the particle density. Recall that 
we could express it through the determinant of the strain. For the strain we 
now have formula (j3.100|) . All we need to do, is to compute the determinant 
of this expression w.r.t. Vabc, the volume form on the body. 

First we introduce k AB , the three-dimensional analogon of the strain: 

k AB = f A ,if B , jg ij (3.106) 

This expression is exactly the non-relativistic version of the strain (see f.e. 
|23j). In this sense we can understand the W l terms in ()3.100|) as relativistic 
corrections to the non-relativistic strain. The inverse of k AB is given by 

k A 1 B = F i , A F>, B g ij (3.107) 

which follows from basic linear algebra. 

According to (|3.23|) the determinant of the strain (apart from a factor 6) 
is given by h AA ' h BB ' h cc ' VabcVa'b'c Using ()3.100|) we write this as 

(2-' _ W a W a ') . . . ft ft: ■ ■ ■ VabcVa>b>c> (3.108) 
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Since Vabc is the volume form on the body, fafbfc^ABC is a volume form 
on the hyper- surface. Thus it must be proportional to the three-dimensional 
metric- volume-form e^. We define a scalar k to be this proportionality 
factor: 

fafbfcV ABC =■■ *e abc (3.109) 

We will find according to (j3.111|) that k is basically given by the particle 
density. From the above equation we find that (j3.108|) can be written as 

_ W a W a '){g bh ' - W b W b '){g cc ' - W c W c ')K 2 e abc e a > bld (3.110) 

This (besides the factor 6k 2 ) is the determinant of (g aa ' — W a W a '). Deter- 
minants of such objects can be easily computed using traces. The result 
is (1 — W 2 ). Thus the particle density n (which is the square root of the 
determinant of the strain) becomes 

n = k(1 - W 2 ) 1/2 = (3.111) 

This equation connects the particle density n with the factor k. 

Finally we turn to the volume form. First we give its expression in ADM 
variables, then we apply our transformation ()3.84j) on the spatial part. The 
result will be a volume form on the body. 

The space-time volume form can be written as 

A/det g^tfx = N^det g^dtd^x (3.112) 

Now ()3.109j) implies that the volume-form on the hyper-surfaces ^det gijd?x 
can be transformed into k' 1 V abc^X a dX B dX° ^ the volume form on B (mod- 
ulo k). 

Collecting all pieces, we can write the action (j3.35)) as an integral over 
the material manifold times (an interval I of) R: 

S M [F] = [ N^hVdtd 3 X (3.113) 

where V is a scalar on the body defined by Vd 3 X := VABcdX A dX B dX c . 
The corresponding Lagrangian is obviously given by 

C M = Nj'hv (3.114) 

A final word on the functional dependence: V is a scalar on the body, iV 
depends on time t and on the deformation. The stored energy is a function 
on the strain bundle over the body, i.e. it depends on the points of the 
body and on the strain. Finally we note that 7 depends on time, on the 
deformation and on the time-derivative of the deformation. This completes 
the listing of the basic objects of the theory in the material description. 
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3.2.3 Reference deformation 

This section gives a definition of reference deformation. We remain as closely 
as possible to the definition used for the reference state in the space-time 
description: 

Definition. A deformation F l is called reference deformation, if it is stress- 
free, i.e. 

tab\ f= f = U (3.115) 

and if the configuration corresponding to F l generates a geodesic normalized 
four-velocity. 

In analogy to the space-time description a reference deformation is called 
natural, if the corresponding elasticity operator \3.41\) satisfies the Legendre- 
Hadamard condition \2. 7| ). 

This is a rather pragmatic approach, since it involves configurations in 
the definition. We have chosen it to keep the definition simple and intuitive. 
The following lemma will give a characterization of a reference state solely 
using deformations. 

Lemma 8. A stress- free deformation F l is a reference deformation, if and 
only if it satisfies the equation 

r\? = (3.116) 

where are Christoffel symbols formed from the metric \S.120i) . The ex- 
plicit form of this condition is given by setting the right hand side of 
equal to zero. 

We prove the lemma by exploiting the relation between space-time and 
R x B: From what we have learned so far it is quite clear that the mapping 

M -> RxB (3.117) 
(t,x) h- (t,f(t,x)) (3.118) 

is a diffeomorphism with inverse 

(t,X)»{t,F(t,X)) (3.119) 

We will apply this diffeomorphism to the geodesic equation for the normalized 
four-velocity. The result will be the geodesic equation for the transformed 
objects. This equation will be of a particular simple form, so that the result 
can be readily read off. 
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For convenience we combine time t and points X A on the body into a 
four-dimensional object 

(X") := (X°, X A ) = (t, X A ) (3.120) 

this notation may seem misleading at the first sight, but it will prove to be 
very practical. It will be clear from the context whether Greek indices belong 
to space-time or R x B. 

Let the transformation of the normalized four- velocity be denoted by U^. 
Then by the standard transformation formula for vector fields we have 

dX^ 

U» = ^—u u (3.121) 

dx v v ' 

(Recall the different meaning of fi and v in this context). For the normalized 
four-velocity this leads to the particular simple expressions 

U°(t,X) =u°(t,F(t,X)) =:u (3.122) 
U A = f A ^u» = $ (3.123) 

Then the geodesic equation takes the form 

= ud t u + T° 00 u 2 (3.124) 
= T A u 2 (3.125) 

where the Christoffel symbols have to be taken w.r.t. the transformed metric 

dx^' dx v ' , 

G ^ = ^QxilQxi ( 3 - 126 ) 

We obtain 

(r ( -N 2 (l-W 2 ) NgijW i F^A \ m2? , 

~ { N gij W^, B gijF\ A F j , B J [6AZ{) 

for the transformed metric where W 2 = gijW l W j . The first thing we can do 
is to determine the factor u. We find that u^u v g^ v = — 1 becomes 

-1 = G^U»U U = -u 2 N 2 (l - W 2 ) (3.128) 

From this we conclude that u = N~ lr y^, so that we can write the geodesic 
equation as: 

= N^d t (N~ 1 ^) + T° 00 (3.129) 
= r^ (3.130) 
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Now we need the expressions for these Christoffel symbols. Note that we can 
rewrite the transformed metric ()3.127j) in ADM form using 

G A B = gi j F i , A F j , B (3.131) 

as three-metric. The inverse is given by 

G AB = g vf A , i f B , J (3.132) 

Then for the shift we take 

y a = g ab Yb (3.133) 

where 

Y B = NgijVTF^B (3.134) 

The lapse finally is given by 

X 2 = Y 2 -G 00 = N 2 (W 2 + 1-W 2 ) = N 2 (3.135) 
Note that the above calculation also implies 

r l = (1 - ^) (3.136) 

Now we can apply the formulas for the Christoffel symbols in the ADM- 
formalism from the appendix to rewrite the expressions entering the equa- 
tions as 

7^(iV a 7~ 1 ) - 2iV 2 r° = W^dtiN 2 ^ 1 ) - Y A ((N 2 ^ 1 )^ +2Y A ) (3.137) 
2N 2 T A = -Y A d t (N 2 1 ~ 1 ) + [N 2 G AB - Y A Y B ] ((N 2 ^ l ), B +21^^(3.138) 

From this we can show that it is sufficient to look at the second equation, 
since it implies the first. This is seen by contraction of (|3. 138(1 with Y A : 

27V 2 r^ o y A = -Y 2 d t (N*y~ 1 ) + (A 2 - Y 2 )Y B ((JVV 1 )^ +2Y B ^ (3.139) 

which is nothing but N 2 y~ l times the right hand side of (|3. 137(1 . Thus 
()3. 138(1 vanishes if and only if the normalized four-velocity is geodesic which 
completes the proof on the lemma. 
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Chapter 4 

Dynamics of Elastic Materials 



This chapter deals with the dynamics of elastic materials. We will derive the 
equations of motion governing the evolution of the materials, which are the 
Euler-Lagrange equations arising from variation of the actions given in the 
previous chapter. 

The linearized equations (in the space-time picture; the same is however 
true for the material picture, since the linearized equations are the same in 
both descriptions) are shown to allow for global solutions, which will be given 
in closed form. This will give some insight into the propagation of elastic 
distortions. Then we will show that the non-linear equations are hyperbolic 
in both descriptions, i.e. that the Cauchy problem is well-posed. 

After this we generalize this result in two directions: First the self- 
interaction due to gravitational forces is taken into account by analyzing the 
equations of elasticity coupled to the Einstein equations. This will be done 
in the space-time description. On the other hand we introduce boundaries, 
i.e. we specialize to materials of finite extension. This will be done in the 
material picture, since the boundary-conditions take a much simpler form in 
that setup. For both cases (the self-gravitating as well as the finite-volume) 
we show the Cauchy-problem to be well defined. 

For all the well-posedness results in this chapter we employ the theorems 
given in the chapter El 

4.1 Equations of motion 

This section is devoted to the derivation of the equations of motion, which 
are the Euler-Lagrange equations arising in the context of the Lagrangians 
(jH.Hfjjl and (j3.114j) . Special attention will be given to the principal part, since 
we know from chapter 2 that this part of the equations is the relevant one 
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for well-posedness. 

Let C denote a Lagrangian depending on the field ^(x^) and its first 
order derivative. Then (under certain regularity assumptions on the La- 
grangian) (p A is a critical point of the action if and only if it satisfies the 
Euler-Lagrange equations: 

d dC dC . - 
= f4 1) 

dx»d(j) A ^ d(f) A v ; 

This is a quasi-linear system of second order for the unknowns . In the 
following we will derive these equations for both the material and the space- 
time description. 

4.1.1 Space-time description 

We use the Lagrangian ()3.36|) to obtain the equations of motion for the 
space-time description. Recall the form of this Lagrangian: 



C s = ne(f, h)^-detg (4.2) 
Then the Euler-Lagrange equations can be written as 

^ABf >/«> : ~ ~QfA QfB f -,nv= Ra (4.3) 

where the lower order terms are given by 

_ d(ne) d(ne) B 1 dC s 

df A df A ^fz J y/=tiBbg a,t df\ l4 - 4j 

Note that the differentiation in the last term on the right hand side is a 
partial differentiation acting on those entering explicitly (via the metric). 
They do not lead to second order expressions. 

Remark. Note that we could as well use A^g instead of A^ in the principal 
part. But this would lead to problems once boundary conditions are involved. 
(This can be seen in section \J~b] there we require a certain matching be- 
tween principal part and boundary terms, which is violated for the considered 
symmetry assumption ) Therefore we stick to the full principal part. 

It is interesting to note that the fact that the principal part comes from 
a Lagrangian is reflected in the following symmetry of the principal part: 



^-AB — BA 



(4.5) 
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At this stage we pause to note a deep connection between the field equations 
as given here and the energy-momentum tensor: In general an alternative 
way to derive field equations is to claim divergence-less-ness of the energy 
momentum tensor. We have seen in section (jH.1.4j) that -u M Vt,T^ = is 
identically satisfied due to the continuity equation. In the following lemma 
we will see, that the vanishing of the remaining components of the divergence 
is equivalent to the field equations derived from the Lagrangian via variation. 

Lemma 9. A configuration, which solves the Euler- Lagrange equation has a 
divergence- free energy-momentum tensor and vice versa. 

As mentioned before we only need to care about the components of the 
divergence of the energy-momentum tensor orthogonal to the material veloc- 
ity. For this rest we first recall the definition of the energy-momentum tensor 
from (jnHU): 

W = jprf A >*-f*!i ( 4 - 6 ) 

Then 

i i 

V7 _ J>H _|_-TM fa _ rpn _ rpp , rpa IJ./3 rpfj.a (a n\ 

where the last equality hods due to the symmetry of the energy-momentum 
tensor. The first term on the right hand side can be written as 



= «7Z- U S\ -^71 f\ -duP (4-8) 



This is already very close to the desired result. As in section 13.1.41 we can 
write 

duP=^f A ,af B ,pg a ^, (4.9) 

since an explicit x-dependence only entered via the metric in the strain. The 
other two terms on the right hand side of (|4.7jl can be written as: 

1 dp fA n/3 n 1 ®P fA \ a (Am\ 



where the first term can be written as 

~'-f A ,» (4.11) 



dp (V-dety,.,, r y 



df A , a V- del g 
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This term together with the first two terms in ()4.8|) is the desired result. We 
have to show, that the second term in (j4.10j) cancels the last term in (|4.8|) . 
To see this recall from section 13. 1 .41 that 



dp 



f A = o ® P f A a afl f B (A 121 

qja J > v dh AB ' ,a ,u y^-^J 

Application of this to the two terms in question shows, that they cancel. The 
remaining parts can be cast into the form 



v M v 



\J— det g 



dp \ i — dp 



-detg —. ,„ -W-detg 



df A ,J^ v *df 



f A w (4-13) 



This together with the continuity equation implies the lemma. 

We need good information on the principal part in order to apply the 
existence theorems from section 12.11 Thus we will give a more explicit form 
of A^ B . In a first step we compute 

dne n dh MN 

Thus 

AW rr dh MN dh KL . v . . 

= + nUMNKL W7,W7, + (415) 

There are basically three parts: The first involves the stored energy, the 
second the elasticity tensor, and the last term is linear homogeneous in the 
stress. More precisely the following holds 

/ dn dh MN dn dh MN d 2 h MN \ 

For a stress-free state $^| T=0 = and the principal part consists only of 
the first two terms in (|4.15jl . 

To obtain a higher level of detail, we have to compute the involving partial 
derivatives. The easiest is the variation of the strain w.r.t. the configuration 
gradient. It is straightforward to see that 

^ = 25ff\ a g^ (4.17) 

Then the second variation yields 

f)2t,MN 
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The only missing ingredients are the variations of the particle number density. 
We have seen in section 13.1.41 that 

On I 

-nhAB (4.19) 



dh AB 2 

with Hab being the inverse strain. Thus 

dn 



nhANf N ,a!r (4-20) 



df A „ 

From this expression we can compute the second variation 
d 2 n , „ M „„, nM „„ dh 



nh BM f M ,p 9 ut3 h A Nf N ,a gT + n^j^f N , a <T + nh AB g^ 



df A „df* w ' p ° — "df* w 

Since Hmn is the inverse strain, we can use the formula 

dh AN , dh CD 



h A cwj7^h DN (4.21) 



This leads to 

-2h B[A h N)D f D , a g av (4.22) 



df B „ 

Inserting this into the formula for the second variation gives 
2nh BM h AN f M ,a /% g a[v gM + nh AB {gT - h MN f N , a / M )/3 g^g^) (4.23) 

This is where possible antisymmetric (in the upper index-pair) portions of 
the principal part enter. 

We can further simplify using two observations: first note that 

<P B :=h BM f M , a g a » (4.24) 

is the inverse of the configuration gradient on the orthogonal space of the 
four-velocity. In formulas this means 

r B f A ,u = S A (4.25) 

and 

<P u B f B „=S; + uX (4-26) 

While the first equation is immediate from the definition of the (jf B , the 
second equality is shown by contracting (t) u B f B ^ (which is symmetric w.r.t. 
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the metric) with the four- velocity and the configuration gradient respectively. 
The first contraction obviously vanishes, while the second gives 



'B 



f B i^ f A -,v — fu^BMf M ,a 9 aU f B ifi — f A ,fi (4.27) 



This proves the given formulas. 

On the other hand we can use this new information to simplify the 
bracket-terms in ([4.23)1 . We can write them as 

<T - h MN f N , a f M ,p gTsf" = -u"u v (4.28) 
Thus the final version for the second variation of the particle density is 



d n 



2nfy B ] - nh AB u^u u (4.29) 



df A ,adf B , t 



With this expression at hand, we have all what we need to cast the principal 
part into its final appearance: 

A% B = (2nfy v] B - nhABifu") e + AnU ANBL f N , a g a » f L , + % (4.30) 

This will be the form of the principal part relevant in the next sections. We 
have the following expression for ^ AB : 

% = In {r BN f N , a g^^A + W*,« 9^b + W) (4-31) 

Note that besides the antisymmetric very first term the stress-free part of the 
principal part is split into a term in it-direction and a term orthogonal to the 
four- velocity. One can continue in this direction by splitting the stress-terms 
in a similar manner. Recalling the definition of the <j) B , we find that we can 
split the metric according to 

<T = h MN f M , a g a »f N ,p of* - u»u v (4.32) 

We also find that the antisymmetric term in the principal part is orthogonal 
to the four-velocity: 

fyi = h M[ Ah B] Nf M ,a g a »f N ,p /" (4.33) 
Thus we can split the principal part into two orthogonal components: 

A% B = -n (eh AB + 2t ab ) vPu" + 2nM AMBN f M , a g a »f N , p g^ (4.34) 
where 

M AM BN = 2U A MBN + T AB h,MN + T BN h AM + T AM h B N + th M [ A h B ] N (4.35) 
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Note the following symmetry of M AMBN : 

M A mbn = M BNAM (4.36) 

The other symmetries of the elasticity operator (|3.1.3|) are lost by virtue of 
the other terms. 

We finally give the equations of motion in a convenient form 

(-nefiABWu" + AnU AMBN f M , a g^f N , p + d^d u f B = R A (4.37) 

where R A again denotes the lower order terms. 

4.1.2 Material picture 

From the Lagrangian (j3.114j) we derive the equations of motion for a elastic 
medium in the material description. This Lagrangian had the form 

C M = N-y-hv (4.38) 

In the above formula N is the lapse, V is the scalar volume element and 
7 _1 = (1 — W 2 ) where W was defined by NW l := F % + Y\ where Y l is the 
shift. 

As in section 13.2.31 we combine time t and points X A on the body into a 
four-dimensional object 

(X") := (X°, X A ) = (t, X A ) (4.39) 

Whenever we use the material picture (and whenever the contrary is not 
explicitly emphasized), we apply this notation. Greek indices thus range 
from zero to three and denote points on K. x B. 

In this notation the Euler-Lagrange equations can be written as 

Algd&Fi = Ri (4.40) 

where the lower order terms are given by 

! ^ Kt - dF* dF\,8Fi V y J [ } 

As before in the space-time description we will no longer care about this 
lower order terms. They are (as long as they stay smooth) of no importance 
for the question of local well-posedness and only come into play when the 
longtime-behavior is analyzed which is far beyond the scope of this thesis. 
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The principal part reads 



i 



Note that 7 only involves the time-derivative of the deformation and none of 
its other derivatives. Thus the F 1 ^ can only be found in the stored energy 
function (via the strain). With this property of the Lagrangian variation 
looks easier than in the space-time-description. But as we will see, this 
case is in fact the more complicated one. This is because the deformation 
gradients enter the Lagrangian in a less explicit way, which is a result of the 
transformations involved in the transition between the descriptions. 

In analogy to the space-time-description the principal part is given by 

where $^ is again linear homogeneous in the stress. It is of the form 

V? = r AB L-k dVB + + ^±^1) (4.44) 

" AB y dF* v BFi,„ BF* v dFi,„ BFi„BF* v ) y ' 

As for the space-time-description we will derive more explicit formulas by 
computing the variations. The expressions finally obtained will be the coun- 
terparts of ()4.30|) and (|4.34|) in the material description. 

In any case we need the variations of 7^2 = (1 — W 2 ) a and of the strain. 
For the former we first compute 

_L_ = -J^w a g ai (4.45) 
Q F i N y v ' 

From this the second variation can be obtained 

<9 7 5 _ / 7 5 TJl&jxrb 7^ ,,/, 



. . ±-W a W b - i^g ab ) g ia gib (4.46) 

OF* OF' V N 2 A - iJ Jj 

Now we turn to the variation of the strain: Note that the / , a entering 
the formula (j3.100j) for the strain are by means of ()3.85j) the inverse of the 
deformation gradient and thus only depend on F a ,A and not on F l . 
The first variations are given by 

f^ = -f/^/^ (4.47) 

pjUAB 

-2f, a f (A n f\ b (g ab - W a W b ) = -2f A n h B » (4.48) 



8F l 
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where we have used the variation of ()3.85|) w.r.t. the deformation gradient 
to compute the variation of f A , a - From the above expressions for the first 
variations the second variations follow: 

r) 2 h AB 9 

(4-49) 



<fIrW - 2f J n f*, h B » + 2f A n (f% h " + h B)Jf tj ) (4.50) 



dF^jdFi, 



<)2hW ^/ , .;/ /; '..,/', /'.,, 4r" (4.51) 



dF i , I dFi N- 

Using these formulas we can write down the principal part. Since the re- 
sulting expressions are rather lengthy we split into time-components, spatial 
components and mixed components: 



Time components 

We first compute the time-portion of the principal part (|4.43j) : 



Q2 -i XfrAB fjUMN 

4? = + 7 -**Wb • + *S? (4-52) 

13 dF^Fi OF* dFi tJ 



From the formulas we have obtained for the variation of 7 and of the strain, 
we find that this can be written as 

1 _i 

= - (WiW j + 9ij ) + ^u ABCD f A , a w a f B n f c JC w c f D , 3 +<&°° 

where the index at W was pulled with the induced three-metric: W{ = gijW 3 . 
For application it will be more convenient to replace the original principal 
part Afj by a new object defined by 

B^ B :=A^F\ A F\ B (4.53) 

Since the deformation gradient F l , A is an isomorphism, those two objects 
can be viewed as equivalent. For example positivity of Afj implies positivity 
of B^ B and vice versa. The reason for introducing B^ B will be clear, once 
we write it down. Using the formula (|3.101|) for the representation of the 
inverse strain we obtain 

1 _i 

BZ = -hAB^ + ^U AC BDf C ,cf D , d W c W d + ^ A ° B (4.54) 
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where ^ A V B := ^F\ A F j , B . To complete the time part we note that $°° is 
given by 

i 

= 2 ^ AB (-7-V A » f B ,j +2W {l f A , 3) f B , b W b ) (4.55) 

Thus we obtain 

i 

= Zj^tmn + 2W i F\ {A 5%f, b W b ) (4.56) 
Spatial components 

Next we come to the spatial part. This is easily obtained, since 7 does not 
depend on F\ A - It is given by 

aar _i t , dh KL dh MN ^ AR . . 

< = 7 . gl „,„ — — + «y (4.57) 

This takes the simple form 

A AB = ^U MNKL f M n h NA f K „ h LB + (4.58) 

which becomes 

Bff = ^- h UjMj N h MA h NB + Vjf (4.59) 
It remains to compute § AB and ^ AB ■ The former is given by 

*f = 2^r MN (f B n f M , 3 h NA + f M n f N , 3 h AB + f M n f A , 3 h NB ) (4.60) 
From this the latter is found to be 

*ff = {5 B 5fh NA + 8f8jh AB + 5f5 A h NB ) (4.61) 
This completes the computation of the spatial terms. 

Mixed components 

Now we turn to the mixed components. They are obtained from 

.OB *TT dh KL dh MN , A/%n . 

^ = T*U MNKL ^ WJ - + *°° (4.62) 

The other mixed components then follow from the symmetry of the principal 
part: 

< = A f ( 4 - 63 ) 
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Applying our formulas to compute the included variations we obtain 



A°f = ^U MNKL f M , i h NB f K , a f L , j W a + <!>°f < 1.0 J) 



This leads to 



Bfj = ^U INJK h NB f K , a W a + ¥°? (4.65) 



Finally we compute 



N 

and 



*?f = ^r MN [2 7 -V (B , f\a W\f M ^ +h BM f N , j Wi] 



X T,0B 2 7 2 



2 1 ~ 1 5ff N \ a W a Sf + h^S^WiF^A (4.66) 



With these formulas for the mixed components of the principal part we can 
proceed to the final task of this section: 



Full principal part 

We will add up the partial results from the proceeding subsections and will 
derive the full principal part including the stress-terms in analogy to ()4.34|) . 
Note that the key point there was to split the principal part into a portion 
proportional to the four-velocity and an orthogonal remainder. Here we will 
proceed along the same lines. This means that first of all we have to identify 
the analogous objects for the material setup. 

The form of the principal part ()4.43j) suggests, that one preferred direction 
will be dt- To give an explanation for this fact we recall our procedure in 
section 13.2.31 There we saw, that the transformation (t, x) <-> (t, X) given 
by the diffeomorphism 

(t, X) = (t, f(t, x)) (t, x) = (t, F(t, X)) (4.67) 

maps the normalized four-velocity to 

^ u 5% (4.68) 

Thus dt is nothing but (a multiple of) the normalized four-velocity mapped 
to ix B. From this one can guess that the transformation of f A ,^, (viewed 
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as a collection of vector fields on space-time) will determine the remainder. 
We find that 



This leads to the definitions 



MA 



» C " :=S^F a f c , a +5 M h MC 



(4.69) 



(4.70) 
(4.71) 



The coefficient entering the definition of U^ 1 guarantees the normalization. 
In order to keep the calculation as simple as possible, we assume Gaussian 
coordinates on space-time. This choice makes the computation a lot easier 
to overlook. 

Even for this simple choice of coordinates the actual calculation is rather 
lengthy, we will only catch the most important steps. We find it more con- 
venient to replace A^ by B^ B - This new object is given by 



AB — °0 °0 



2r CD ( 2 7 f F c f c , c f D , {l F j} F\ A F\ B - 7 -i<%'5 



a 7716 



-/Us7 2 e + 47 *U AaBb F a F 



(i r 3) J 



CIuDJ 



4~f-iU AC BDh ul h 



+2^tcb (S J A S C B h DI + 5 c A 5 D B h IJ + fi8 I B h EJ ) 



+ 



i^U ACBD F a f c , a h DJ + 



+ 



+ 2r CB ^F k g kn F n , A h CJ + 2y~h ( £ f J \ m F r ' 

[^U AC BDF a f D , a h CI + 
+ 2t ca ^F k g kn F n , B h CI + 2 7 -^g7 /) ,™ F r ' 
In analogy to ()4.34|) we try to decompose this in the following way 
B% B = -A AB UW + A^Uacbd^'^ + B ACBD <P C ^ D » 



(4.72) 
(4.73) 
(4.74) 
(4.75) 
(4.76) 
(4.77) 
(4.78) 
(4.79) 

(4.80) 



That this is indeed a valid ansatz will be seen in the course of the computa- 
tion: we will find that there are indeed no other terms than the ones given 
above. In particular there will be no U^(f) Cv '-terms. 

There are in principle two ways to proceed. One may guess A AB and 
then compute B A cbd, which is again a check of the validity of the ansatz for 
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Aab- On the other hand one may try to construct Bacbd by putting terms 
into the form <fi c ^. The remaining terms form then Aab- in both cases we 
have to check that the coefficients of the elasticity tensor are of the required 
form. This will be our first task. The parts involving Uacbd are 

5^ A 1 ^U A aB b F a F b + 6?6 v j4'y-*U A cBDh CI h DJ + (4.81) 
+6S6 v J 4r r -*U A cBDF a fZ: h DJ + 6ffi4ry-*U AC BDF a f?h CI (4.82) 

Now it is straightforward to find that this is of the desired form. To discuss 
the remaining terms we decide to make a proper ansatz for Aab since this 
turns out to be the simplest route. In reminiscence of the analogy for the 
part involving the elasticity operator we pick 



Aab = 7 2 (h A Bt + 2t A b) 



(4.83) 



Both these terms come from the variation twice w.r.t. F. The first term 
in the bracket is found immediately. The second requires a little more care 
since it appears with the (at the first glance) wrong factor. This is overcome 
by splitting 7 _1 = 1 — F 2 . The remaining terms then are 



2r CD [2l*F°£f$F j) F* A n + T^F 
[2^r CD (S J A 6 c B h DI + 5 c A 5 D B h IJ + 5 c A 8 B h DJ ) 
\2rcB (^F k g kn F n A h CJ + Zy-Wff$F» 



+ 
+ 



[2t ca (r*F k g kn F%h CI + 2^5^^) 



(4.84) 
(4.85) 
(4.86) 
(4.87) 



In the next step we will evaluate the Kronecker deltas and split according to 
the indices of the stress. This leads to 

2r AB (« 7 §F 2 + S^fy-W J + «7"* f J J™ + S^j^f^F m ) +(4.88) 



+ 2t CA 



6^F c ffF a g ab F B + 8$r f1 -H 1 B h VJ + (4.89) 



/ uCJ 



[^F k g kn F-h CI + 7 -^i/^ m )] +(4-90) 



5%5» r*F c tfF a g ab F b A + 5^r^ J A h CI + (4.91) 

+65% ^F k g kn F2h CJ + 7 -^i^F m )](4.92) 

All left to do, is to simplify the terms in the brackets. The first term is the 
easiest: it is proportional to tab^cd<^ Cii ^ iDv ■ The factor is obtained as 



TAB{...)=l-^TABh C D<f^ DV 



(4.93) 
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The other two terms are slightly more difficult but can also be handled with- 
out much trouble. We do the calculation for the second term, the third works 
just the same way. We want to show that 

t C a{- ■■) = l^r Ac h BD( p c y D » (4.94) 

To prove the validity of this claim we expand the right hand side: 

h BD f^ Dv = h BD f c f D 5^ + 5 N B h MC 5 M 5» N - <^/ c « - S^ M h BD f D h CA 

The second and the third term match the second and last term in the bracket 
in ([4.89)1 . The other terms need some handling. 

h BD f C f D = F b B (g bd + 1 F b F d )F d f^F c = 1 F b B g bd F d f^F c (4.95) 

which shows that the first term equals the first term in ()4.89j) . The remaining 
last term can be dealt with along the same lines. 
The very same arguments show that 

t cb {. ■ ■ ) = l^r C Bh AD( p c ^ Du (4.96) 

So altogether we find that the principal part can be written as 

Kb = -T l (h AB e + 2t ab ) U»U» + 2 7 "^ M ACBD <t> c "<t> Dv (4.97) 

where 

M AC bd = 2U AC bd + r AB h CD + r BD h AC + r Ac h BD (4.98) 

We find that the principal part takes the very same form as in (|4.34|) . apart 
from the missing h A [chr)]B term. Nevertheless we only have the symmetry 

M ACBD = Mbdac (4.99) 

Remark. The same result could be obtained including lapse and shift. One 
would just have to replace F l by W l and add the shift N at the right places. 

As was the case in the space-time-description, the "spatial" portion of the 
principal part lacks the full symmetry of the elasticity operator for general 
states. Only for special cases, f.e. stress-free states of the material, this full 
set of symmetries is given. 

This ends the derivation of the field equations. We can now turn to the 
discussion of local well-posedness for certain situations. Before doing so, we 
make some comments on the existence of conserved quantities. 
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4.2 Conserved quantities 

In this section we discuss conserved quantities. By the Noether theorem (see 
f.e. |51j ) we know that the existence of symmetries implies the existence of 
conserved quantities. Under the assumption of staticity we will derive such 
a conserved quantity, which will be the energy of the system. We will then 
linearize this energy at a stress-free reference solution. We show that the 
linearized energy is of the form "kinetic plus potential terms" . We will also 
show that the linearized material Lagrangian is of the form "kinetic minus 
potential terms" with the same expressions for the kinetic and potential 
terms, which reflects the connection between Lagrangian and Hamiltonian of 
classical mechanics. Finally we rewrite the linearized potential terms using 
a covariant derivative. This allows for a direct comparison with the non- 
relativistic case. Note that the derived energy is conserved for both finite 
and infinite matter distributions (this last statement of course depends on 
the boundary conditions; in particular we will see that it is true for natural 
boundary conditions as introduced in section 14.6.2 J) . 

4.2.1 Symmetries- conserved quantities 

Given a Killing vector we can contract it with the Energy-momentum 
tensor: T^ u ^ v . This new object is divergence free, if the Energy-momentum 
tensor itself is divergence free, which we assume in the following. (We have 
seen that the energy-momentum tensor is divergence-free for solutions of the 
Euler-Lagrange equations) We have 

V„ (T^Q = (4.100) 

Integrating this over a region Q in space-time, we obtain 

= J V M (T^£„)vol (4.101) 
h 

We can use Stoke's theorem to convert this into a surface integral 

= J V, (T^ v ) vol = J T^voL, (4.102) 
n an 



The volume element on the boundary can be written as vol M = det g^n^dPx. 

is the outer unit normal, while denotes the induced metric on dQ. We 
find that 

= ^ T^&n^l det g (i) \d 3 x (4.103) 
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Assume Q to be a cylinder. Then dQ consists of three parts: The bottom and 
the top and the mantle. If one can show that the integral over the mantle 
vanishes (due to boundary or fall-off conditions for the fields), one finds 
(taking into account the sign change in the normal) that the integral over 
the bottom coincides with the integral over the top. A conserved quantity is 
found. 



4.2.2 The conserved energy 

Assume the space-time to be static and that the Killing vector has norm 
— 1. Then we can introduce coordinates (t, x l ) such that t is the parameter 
along the flow of the Killing vector field and the x l are coordinates on the 
hyper-surfaces t = const. The metric can be written in the following form: 

d s 2 = -dt 2 + g ij (x)dx i dx j (4.104) 

The Killing vector is then given by <9 t . Integrating over the world-tube of 
the material (i.e. the inverse of the body under the configuration map), 
we assume that the integral over the mantle of the tube vanishes. This is 
guaranteed by fall-off or by the natural boundary conditions introduced in 
section FPU Thus the integral of T^^n^y/l det | dx over any space-like 
slice of the world-tube is conserved. By choosing the slices to have constant 
t we arrive at an energy, which by construction is conserved in time: 

E:= J T^^^\drtg~\d 3 x (4.105) 

f-H*m 

here f~ 1 (t)(B) is the inverse of the body B under the map / at the instant 
t. Modulo the determinant of the induced metric the integrand is given by 

T^Uv = ee(pu„u v + 2nr AB f A ^ f B , u ) (4.106) 

For each value of t this is a scalar on the hyper-surface t = const. We can 
compute its pullback along the deformation map F l . We had the following 
relations between quantities in the material and spatial picture: 

= ^=v M (4.107) 



V- 



■v 



From this it follows that 



V »Q^ = d t + piQ. (4.108) 
- V 2 = l-F 2 (4.109) 



= = (4.110) 

VI- F 2 
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and that we have 

T^U* = + 2nr AB E A E B (4.111) 

where 

E A :=f A , li e = f A = -f A nF i (4.H2) 

The pullback of the volume form on the t = const hyper-surface is given by 
(the volume-element on B is chosen such that V — 1) 

(8F k \ (I — F 2 )^ 
dX = K - ^-d^X (4.113) 

Now the energy can be rewritten as an integral over the body £>: (The integral 
coincides with the integral of the pull-back over the base space) 

E = / | + 2TAB r f A n f B ,i FFn d 3 X (4.114) 

Using the definition of 7 we can abbreviate this: 

E = y 7 1 (e + 2r AB f A n f B „ F i F j ^j d 3 X (4.115) 

B 

This energy is conserved in time. It is built solely from material quantities. 



4.2.3 Linearized energy 

In this section we will expand the energy in a neighborhood of a static and 
stress-free natural deformation. By this we mean that the reference con- 
figuration (the existence of which we take as given) is assumed to satisfy 

u^ = dt (4.116) 

(An equivalent way of putting the above equations is to write u^d^ — d t . It 
is however useful to stick to the notation used in ()4.116j) . For conventions 
regarding the reference state see section I3.1.5|) By the definition of the four- 
velocity ()4.116j) includes that the reference configuration and hence also the 
reference deformation is time-independent: 

^ = (4.117) 

We will find that the first order variation of the energy vanishes and that 
the second order terms are of the form kinetic energy plus potential energy. 
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Remark. As will be seen in later sections this particular reference state is a 
solution of the field equation satisfying the boundary conditions (both fall-off 
and natural). 

As usual we will write = for relations which only hold for this natural 
deformation. The perturbation will be denoted by 5F\ 

To see that the first variation of the energy vanishes note that we can 
easily compute the following partial derivatives: 

<h ' < 7 § ^ = (4.118) 



9 ^ = \l F ap b o 

dF i 2 





(4.119) 


de 8h AB a 
3FV ~ rAa dF\^ ~ 





(4.120) 


de _ dh AB a 
dF i ~ dF i 





(4.121) 



The first variation of the second term in the bracket in ()4.115|) vanishes, since 
it is quadratic in F % . There is also the Piola stress tensor in the third term, 
so it only contributes to third order. We may omit it in this analysis. We 
have found that the expanded energy is at least quadratic in the deformation 
field. 

We now compute the second variation of the energy: Therefore we need 
the second variations of 72 and the stored energy e: For the former we obtain 



<9 2 7 



2 



dF\,dFi, v 



« 3 7 ^ + 7 5 ^ = «fo (4-122) 



All other variations vanish for the obvious reason that they still involve time- 
derivatives of the field, which are zero for the natural reference deformation. 

The second variation of the stored energy with respect to the deformation 
gradient is 

d 2 e TT dh AB dh CD a dh AB 8h CD 

dW^dFT„ - UABCD dF^dFT- u +TAB ^ UABCD dF^dFr u (4123) 

We can write out the variations of the strain and get 
d 2 e 



dF^dFi, 



W ABCD h MB f A n h ND f c , 3 « (4.124) 
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The second variation with respect to the deformation reads 

d 2 e TT dh AB dh CD a TT dh AB dh CD 

Uabcd arii k^t + tab {■■■) = Uabcd a =j „ „, (4.125) 



dFdFi dF i dFi gpi dpj 

writing this in terms of the deformation we obtain 

<)1( = U A BCDf A ,af B , b g ab nf C ,cf D , d g Cd , 3 (4.126) 



dF*dFi 

Finally we need to take care of the mixed term: 

d 2 e TT dh AB dh CD , « q dh AB dh CD . 

Uabcd an ,. a „ t . + r AB (. . . ) = Uabcd an ,. arS (4.127) 



dF\^dFi L dFi ^v-y dFi dF\ 

which gives 

<9 2 e 



2U A BCDf A ,af B , b g ab v h MC f u n 5» M (4.128) 



These are all the terms entering the second variation of the energy. All other 
vanish. Summing up we find that the lowest non-trivial term in the expansion 
of the energy is given by the integral over the following expression: 

eg t3 5F l 5F3 + AUabcd (X ab + Y AB ) (X CD + Y CD ) (4.129) 

where 

X AB = - 1 -f A , a f B , b g ab n 5F* (4.130) 
Y AB = h MA f B n SF\ M (4.131) 
To write this in a more convenient form, we introduce a covariant derivative: 

S7 A V l = d A V i + T) k V k F A (4.132) 

The rk are the Christoffel symbols corresponding to the induced three-metric 
gij. We will show, that we can combine X AB and Y AB into this covariant 
derivative. 

We first define 

VT := 4U AB CDh MA h NC f B m f D , n (4.133) 
Then the second term in ()4.129|) becomes 

VZ N (SF m , M -\h MA F\ B f A , a f B , b g a \ 3 6F-*) ((. . . ) n N ) (4.134) 
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We will show that the second term in the brackets coincides with the con- 
nection terms. A short computation indeed gives 

h MA F\ B f A , a f B , b g a \ 3 5Fi = -F a , M g hm g ab , 3 5Fi (4.135) 

All we have to show is that g a b,j is the only non- vanishing term left from 
the Christoffel symbol T^. The rest, which we will prove to vanish is 
(—g a j,b+9jb,a)- This is antisymmetric in (a, b). We have to show that 

FiMS^V™ = 4:F a ,M g bm UABCDh MA h N( ~ ! f B , m f D , n (4.136) 

is symmetric in this pair of indices. But this is nearly obvious since the term 
on the right hand side becomes 

4h NC f D , n U ABCD f B n f A , k g ka g lb (4.137) 

which of course has the desired symmetry by the symmetries of the elasticity 
tensor. Thus the antisymmetric part of the Christoffel symbol is annihilated 
and we end up with 

\hMAF\ B f A ,af*g a \ 3 ?>Fi = -F a , M r™5Fi (4.138) 

Together with the partial derivative we can combine this to the covariant 
derivative. We finally have the following formula for the linearized energy 
(besides the ground state energy e, which is an additive constant): 

El = J (eg^F'SFJ + VZ N V M 5F m V N 5F n ^ d 3 X (4.139) 

B 

Note that El is of the form kinetic energy plus potential energy. Besides 
the covariant derivatives it coincides with the corresponding non-relativistic 
quantity. Using the strong point-wise stability condition (i.e. if the material's 
elasticity operator satisfies this condition) on the elasticity tensor of the 
natural deformation, we find that the linearized energy is non-negative. For 
non-static, non-isometric deformations it is positive. 

This very last remark needs some additional explanation. 

Definition. We call a linearized deformation isometric, whenever it is a 
Killing vector field w.r.t. the spatial metric g ab . 

Note that it makes sense to speak of the linearized deformation 5F l as 
of a vector field, since by its nature it is an element of the tangential bundle 
over the hyper- surfaces t = const. Going back to the original form of the 
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energy given in (|4.129jh we note that the part involving the elasticity operator 
vanishes if and only if A A< - AB - ) + Y^ AB ^ = 0. We will show that this is the case 
if and only if 5F l is a Killing field for the spatial metric. First note that 
X (AB) + y(AB) _ o if and only if 

—g^iSF 1 + f M , m g m{a SF b \ M =Q (4.140) 

Viewing 5F l as a field on physical space t = const rather than on the body, 
we find that SF\ m = f M , m SF 1 ,m- Lowering the indices in the above equation 
using the spatial metric leads to 

g ab , i 6F i + 2g m{ JF m , b) =0 (4.141) 

which is exactly the Killing equation £$Fg a b = 0. We have the following 
lemma 

Lemma 10. Assume that the strong point-wise stability condition holds. 
Then the linearized energy is positive if and only if the linearized deformation 
is non-static and non-isometric in the above sense. 

Taking into account the energy of the natural state, which is the first (i.e. 
zeroth order) term in the Taylor series, we find that the energy increases, 
whenever one has a non-trivial (i.e. non-isometric) perturbation of the natu- 
ral state. This implies that the assumed natural state (if it exists) minimizes 
the energy if the strong-point- wise stability condition holds. This sheds some 
light into the physical meaning of this condition. 

4.2.4 Linearized material Lagrangian 

Here we will derive the linearized material Lagrangian w.r.t. the given static 
natural state. We will see that it looks very similar to the linearized energy. 
Just the relative sign between the kinetic and the potential term changes. As 
was the case for the energy, we will find, that the first order perturbation of 
the Lagrangian vanishes and that the first non- vanishing order is the second. 

Recall that the fully non-linear Lagrangian for the material picture is 
given by 

% = 7^e (4.142) 

Comparison with the energy shows that the first term differs only in the sign 
of the exponent of 7. The first variation vanishes due to well known reason: 
The variation of the internal energy gives the Piola-stress tensor, while the 
variation of the 7-term holds time-derivatives of the natural deformation. 
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The second variations become 

d 2 £ M ,. 01v w 0h'" D 



Uabcd ktTi ktt, (4.143) 



dF*dFi dF i dFi 

d 2 £ M o „ dh AB dh CD 



Uabcd^-^- (4.144) 



dFt^dFi ^"dP^ dFi 
92 £m o sllsv TT dh AB 3h CD 

mjk = ~ 9i5<5 ° + UABOD d^mz (4145) 



Besides the minus in the last line in front of the spatial metric, this is just 
the same as we had for the linearization of the energy. We thus can directly 
write down the linearized action 

Sml = J (-egijSFSFi + V^ n N V M SF m V N 5F n ) d'X (4.146) 

B 

This linearized Lagrangian is of the form potential minus kinetic energy. 
From this linearized Lagrangian we can read off all the information needed 
for the local existence results proven in section l4~4l If the metric is in addition 
flat, we can introduce coordinates on space-time and on the body, such that 
the configuration gradient and the strain coincide with the Kronecker Delta. 
Then 

VZ N = W ABCD 5 MA 5 NC 5X (4.147) 

so that V^f n N coincides with the elasticity tensor in these particular coordinate 
systems. But then the linearized relativistic Lagrangian coincides with the 
non-relativistic linearized Lagrangian. 

Remark. We have seen that for aflat metric (i.e. elasticity on Minkowski 
space) there is no way to distinguish between the relativistic and the non- 
relativistic theory on the linearized level. Relativistic effects only enter to 
third or higher order. 



4.3 Linearized elasto-dynamics 

We linearize the space-time Lagrangian w.r.t. a stress-free static reference 
solution. From this linearized Lagrangian we will derive the corresponding 
linearized equations of motion. The actual computation will be similar to 
the one done in the previous section. There however we were working in 
the material description while here the space-time picture is used. We do 
not know a priori that these two descriptions lead to the same formulas in 
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the linearized case. That this is indeed true will be seen in the following 
computations. 

Then in the second part of this section an explicit solution for this lin- 
ear system will be derived for isotropic and homogeneous materials. For 
simplicity we assume space-time to be Minkowski space. 



4.3.1 Linearized equations of motion 

The original nonlinear Lagrangian on Minkowski space is according to ()3.36|) 
given by 

L = ne (4.148) 

where the stored energy e depends solely on the strain h AB due to the assump- 
tion of homogeneity. The particle number density n depends again solely on 
the strain, if the medium is assumed to be homogeneous. The strain itself 
is a function purely of the configuration gradient and the space-time metric 
which is known. Thus we can develop both the density and the internal 
energy into a formal Taylor series: 

n = " + a^W^ + df\Jf*, v lf=fSfA,fl +h " at - (4,149) 

' = 1 + 9pr\f=f s f A * +2 d f\ df* w l ^' 6fA ^ 6fB - +h ' at - (4 - 150) 

We recall from section E.l.ll that there always exist which form the inverse 
configuration gradient in the sense that 

f A ,M = Simd^f B , v =ht (4.151) 

where 

hT = rf v + (4.152) 

is the metric on the orthogonal space of the four-velocity. We assume that 
we can choose coordinates in such a way that u — d t in standard Minkowski 
coordinates. This means that the normalized four-velocity is assumed to be 
a hyper-surface-orthogonal time-like Killing field. Thus the corresponding 
configuration f A is a time-independent diffeomorphism between the body 
and the hyper-surfaces orthogonal to u^ 1 . In addition we assume f A to be 
stress-free, i.e. tab — 0. 

Using 0^ and the formulas derived in section E. 1.11 for the variation of the 
internal energy and the particle number, we find that 

n = n (l + A 5f A n -h AB 5f A 5f B + (4.153) 

e = e + 2U AM BNf M , m 5 mt f N , n S nj 5f A n 5f B „ + ... (4.154) 
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Thus the linearized Lagrangian (omitting the zeroth order term which does 
not contribute to the variation) is given by 

U ■= -~eh AB 6f A 5f B + AU AMBN f K \ m 5™f N , n 5 n >5f A n 5f B „ (4.155) 

This Lagrangian is quadratic in the linearized configuration gradient Sf A . As 
in the foregoing section the Lagrangian is of the form kinetic minus potential 
terms. 

Once the linearized Lagrangian is given it is easy to obtain the correspond- 
ing field equations which are the linearization of the original Euler-Lagrange 
equations. We obtain 

-fh AB d t d t 5f B + AU AMBN f M m 8 mi f N , n 5*1 Ms f = (4.156) 

Without restriction of generality we set e = 1 (this can be obtained by 
dividing by e then redefining the elasticity operator). We also use coordinates 
on the body, such that jf 1 = Sf 1 . Then the reference strain satisfies h AB = 
S AB . We also specialize to the isotropic case, where the internal energy 
depends on the strain only via its invariants. Then the elasticity tensor is 
determined by two constants, A and /i, the so called Lame constants (see 
section 13. 1.3j) . In the chosen coordinates we have 

4Uambn = ^am^bn + 2ii>8a(b8n)m (4.157) 

For convenience we write g B instead of 5f B . Then the equation reads 

-g A + ^iAg A + + X)V A (V ■ g) = (4.158) 

Note that this is hyperbolic if and only if the following two conditions on the 
Lame constants are satisfied: 

cj := fjt > and c\ := 2fi + A > (4.159) 

In the following section we will derive an explicit solution for the Cauchy 
problem. 

4.3.2 Explicit solution 

In this section we will derive a solution for equation f)4.158|) for given initial 
data 

g\t=o = J g\t=o = h (4.160) 

For explicit solutions in a more general setting (but only for the case j = 0) 
we refer to |16j . There the requirement of isotropy is dropped. The resulting 
system still has constant coefficients but is coupled in a more complicated 
way. 
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Zero initial configuration 

To keep the calculation clear and simple, we first assume that 

g\t=o = g\ t =o = h (4.161) 

The other case h = will be handled afterwards along similar lines. By 
linearity the solution to the full problem is then the superposition of these 
two solutions. 

Let the Fourier transformation of g be denoted by g. Then the Fourier 
transformation (in space only) of f)4.158|) is given by 

§ A + /ik 2 g A + (jjt+ \)k A {g- k) = (4.162) 

This equation is of the form 

V = -C A ~g B (4.163) 

where the (symmetric) matrix C is given by 

Cab = (fi + ^k A k B + l^k 2 5 AB (4.164) 

Indices are lowered and raised using the Kronecker delta (which coincides 
with the reference strain in the chosen coordinates). The solution of the 
above equation with the Fourier transform of (|4.161j) as initial data is of the 
form 

A /sinC 1 / 2 t\ A 7B 
9 =[—cTJ^) B h ( 4 - 165 ) 

This can be easily verified by direct computation. Note that s ™^ 7 / 1 / 2 is a 
matrix, and that in general it is not easy to actually compute this operator 
from the given C. But one readily checks that C is positive definite if the 
operator ()4.157|) is elliptic, i.e. as long as ()4.159|) holds. But this is assumed 
anyway, so the root C 1 / 2 is well defined. For the actual computation we split 
C into projection operators: 

C = (2fi + \)k 2 P + fik 2 Q (4.166) 

where the projectors 

Pab = k- 2 k A k B (4.167) 

Qab = S AB - k~ 2 k A k B (4.168) 
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satisfy the following relations 

P 2 = P Q 2 = Q (4.169) 
PQ = P + Q = 5 (4.170) 

Using this together with the definition of the sine via its Taylor series and 
the uniqueness of the inverse one can show that (recall the definition of the 
wave speeds q according to (|4.159jl ) 

C 1 ' 2 = Cl \k\P + c 2 \k\Q (4.171) 

C~ 1/2 = AnP + -ttiQ (4-172) 
ci\k\ c 2 \k\ 

sin C 1/2 t = sin Cl |A;|tP + sin c 2 1 (4.173) 

Note that in the last expression the operators P and Q are not affected by 
the sine. When taken together, the above statements imply 

sinC 1 /^ _ sinci|fc|t sinc 2 \k\t ( A -\7A\ 
C 1 / 2 ~ Cl \k\ c 2 \k\ * [ } 

Using that P + Q = 5 we find that g can be written as 

%A _ f^ Cl \k\t sinc 2 \k\t\ - k)kA + sinc 2 \k\ t}iA ^ ^ 



sin c\ 


k\t 


Cl 


k 


3 



c 2 \k\ 3 J c 2 \k\ 

The last term will provide a Dirac delta, the first terms become Heaviside 
functions (and later on by partial integration generate Kronecker deltas too). 
We first investigate the last integral, since it is the easier one. We have to 
compute the following integral: 



H A : = -4— / °^ ^° h A e ik - x d 3 k (4.176) 
(27r) 3 /2y C2 \k\ v ; 

By the standard Fourier formulas we can express this as the convolution of 
F, the inverse Fourier transformation of ( | A; | C2 ) — 1 sin C2I with the initial 
data: 

(2vr) 3 / 2 // = j F(x - y)h(y)d 3 y = J F(y)h(x - y)d 3 y (4.177) 

We thus have to compute the integral kernel F to obtain an explicit state- 
ment. The usual Fourier transformation formula reads 
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For convenience we introduce some notation: k = \k\, £ = |x| and z = cos# 
where 6 is the angle between x and k, so that x ■ k = £kcos#. By a suitable 
choice of coordinate axes we can introduce polar coordinates such that the 
above formula becomes 



oo 1 



p v = p^// *%r**** d * (4 - 179) 



-1 

We can carry out the z integration and obtain 

oo 

2 



F(x) = - — r—j- — - / sin c 2 nt sin k^cIk (4.180) 
{2nyi I c 2 t > J 
o 

The integrand has a primitive, which is not defined at infinity. Therefore we 
cannot hope at all to obtain a classical representation of the above integral. 
But assuming that the initial datum h is a smooth function of compact 
support we can interpret (|4.177j) as the action of the distribution F on the 
test function h. Thus we may as well search for distributional F. To do so 
we note that we can rewrite F using the Euler formula: 

oo 

1 ] ~ (27r)V2 2c2 £i [ + } 



Now we combine the first and the second term as well as the third and the 
last one into integrals over all of R: 



F(x) = "7^)172^ J ( etK(C2tH) ~ etK(C2t ~°) d * ( 4 - 181 ) 

— oo 

The integral is basically the well-known formula for the Fourier transform of 
the Dirac delta distribution. We thus finally obtain 

F(x) = - Q 1/2 JL [5( C2 t + - 5(c 2 t - 0] (4.182) 

This allows a great simplification in ([4.177)1 . Since we are interested in the 
solution for positive t, which corresponds to solving to the future, the first 
Dirac delta does not contribute (recall that £ = |x| > 0). The solution is 
then given by 

11=^-— f 5( C2t _ \ x _ y\)J^L d 3 y (4.183) 
4vrc 2 J \x-y\ 
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We can use the delta to transform the nominator and get it in front of the 
integral: 

One clearly sees that the initial disturbance travels along the sound-cone 
given by c%t — \x — y\ = 0. Next we introduce spherical coordinates. Then 
we can directly carry out the radial integration. This allows us to rewrite g 
as a surface integral over a ball of radius c^t with center x. 



47T C t JdB{x,c 2 t) 



which is t times the average integral of the initial datum over the sphere 
dB(x, c 2 t). Now we turn to the first two terms in (j4.175j) . They are a result 
of the coupling in the equation ()4.158|) and are a little more complicated to 
handle. We have to compute 

I A := * f ( "^'f - (h ■ k)k A e^k (4.186) 



(2?r) 3 / 2 J V 

Since this integral consists of two identical parts we will only compute one 
of them and then add the analogous expression for the other one. First we 
note that 

,A„ik-x j3 7 „ 1 & f sinci|A;|t ,~ ik-x*, 



sin ci 


k\t 


Cl 


A; 


13 



(h ■ k)k A e^k =~^xj ^T(h ■ k)e"«d>k (4.187) 

The remaining integral will again be written as convolution of a function 
derived from the initial data and an integral kernel obtained from the sine 
term. 

Let G denote the inverse Fourier transform of ^rpj^ and let H denote 
the inverse Fourier transform of h ■ k. Then 

sin Cl |fc|t ( ~ k yk-x d 3 k = ( G ( x _ y )H(y)d 3 y (4.188) 



ci\k\ 3 
and consequently 



/ ^r ( ^ ' k ) kAelk ' Xd " k = \j ^i G ( x - y) R Wy ( 4 - 189 ) 



We first determine H. We claim 



1 ~ \Ai 



H(y) = -d A h A (y) (4.190) 
i 
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To prove this we show that by partial integration (recall that h is assumed 
to be of compact support) 

J d A h A e- tk x d 3 x = i J h A k A e- tk x d 3 x (4.191) 

The k can be put in front of the integral. This proves that the Fourier 
transform of d A h A is given by ik A h A , which implies the claim. 

It thus only remains to compute the integral kernel G. We proceed similar 
to what we did with the last term of (|4.175jl . Along the same paths we find 
that 

oo 

(2^)V2 2ci £ J € )dK [4AyZ) 

— oo 

The integrals are nothing but one-dimensional Fourier transformations of 
k~ 2 . Thus we obtain 

G = ^1 ^ {(at + sgn( Cl t + - (cit - sgn( Cl t - 0] (4.193) 

we can give g a better known appearance by exchanging the sgn = 20 — 1 
by the Heaviside step function 6. This gives 

G = ~\fli " [{cit ~ ° e{cit ~°~ {cit + ° e(cit + 01 (4,194) 



Next we have to take the derivative with respect to x. 
dG Pk x r . Pk 1 



dx V 2 ci<^ 3 J V 2 



|e( Cl t-o + |e( Cl i + o 



(4.195) 



There are also Dirac deltas, but they only appear together with the argument 
(i.e. as x5(x)) and therefore do not contribute. Taking into account the terms 
in the first bracket we end up with 

^ = -/|^[ ( c ^ + - ^ - 0] (4-196) 

Note that the constant c\ only appears as an argument in the Heaviside 
functions. Having this expression for G, one half of / is determined by 
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The second half of I is just the negative of this expression only with c\ 
replaced by ci- At this stage we assume that c\ > 02- When these two 
expressions are added up, we find that the Heaviside functions add up in 
way, which only leaves B(x, cii) \ B(x, c^fi) as integration domain. (If C2 > C\ 
the balls would simply be interchanged.) Thus we have 

I A = ~ [ { 4^^(d ■ h){y)d?y (4.197) 



4?T J B(x, Cl t)\B(x,c 2 t) \ X ~ V\ 3 

This is possible, since (as already mentioned) the wave speeds q only appear 
as arguments of the Heaviside function. Now we can add // to obtain the 
full solution g: 

9 A = ^\il h A (y)dS y -±[ T^T^ • h)(y)d 3 y 



47T C%t JdB(x,c 2 t) J 47T J B(x,cit)\B(x,C2t) \ X ~ V, 

This formula still needs some cosmetics, since there appears a divergence of 
the initial data. This problem is easily overcome by partial integration: 

9 -h A <f> B = l h A ^r(f> B - I h A N A <p B (4.198) 



dy A Jv dy 



dv 



Note that in contrary to what was done before, we must not omit the bound- 
ary terms, since integration takes place over a finite domain. In our case V = 
B(x, C\t) \ B(x, Cqf). Thus the co-normal is given by \x — y\N = ±(y — x) 
where the sign depends whether we are at the outer or inner boundary. The 
derivative of 4> = (x — y)\x — y\~ 3 with respect to y is given by 

9 - B - 1 (-5 B A + Z {x - yr{x - y)A \ (4.199) 



dy A \x — y\ 3 \ \x — y\ \x — y 

For convenience we introduce 

h a = (y-x) A (4200) 

\x - y\ 

At the outer boundary this n coincides with the unit normal and on the inner 
boundary it is just the anti-normal. We write 

^ = |^ + 3-^) (4.201) 

This means that the volume term on the right hand side of 1)4.198)1 becomes 

I - — l —^{'in A n B -5 1 l)h A d i y (4.202) 

J B(x,dt)\B(x,c 2 t) \ x ~ y\ 
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The surface integral over dV splits into two integrals over spheres of radius 
Cit with center x. They are given by 

1 



eft 2 
1 



n B n A h A {y)dS y (4.203) 

dB(x,cit) 



n B n A h A (y)dS y (4.204) 

dB(x,C2t) 



Adding the last three integrals up according to (j4.198j) and combining the 
result with the original formula, we arrive at the final form for the unknown 
9- 



9 B = i- t -^-^n A n B - 5 B )h A (y)d 3 y + (4.205) 

4vr J B(x, Cl t)\B(x, C2 t) F -y\ 

+ T^4t / (5 B -n A n B )h A (y)dS y + 

47T C 2 t JdB(x,c 2 t) 

+ / n B n A h A (y)dS y 

47T C x l JdB(x,at) 

According to this formula, a given initial disturbance at a point y is noticed 
by an observer sitting at x in the following way: there will at first be a sharp 
signal followed by a continuous wave, which finally ceases with another sharp 
signal. There is also a grave difference in the nature of the sharp signals. One 
(ci) wave is transversal, while the other one (02) is longitudinal. Usually the 
longitudinal waves are called pressure waves and the transversal waves are 
known as shear waves (see f.e. [23 or another introductory textbook). 

Zero initial momentum 

In order to complete the picture we investigate the behavior of solutions to 
(j4.158|) with the following initial data (compare to (j4.161j0 : 

9\t=o = J 9\t=o = (4.206) 

Due to the linearity the solution for general initial data is then the sum of 
the two solutions obtained for the initial data (j4.161|) and ()4.206|) . For the 
above data, the solution of the Fourier transformed equations (I4.162J) is 

g = cos C 1/2 tj (4.207) 

The C can again be split into the projection operators P and Q which allows 
us to explicitly compute cos C 1//2 t. We finally obtain 

,,. ~ /cosci|A;|t cosc2|/dt\ ,~ 
g = cosc 2 \k\tj + y-^r TjjT 1 - ) ■ k)k (4.208) 
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The solution g is the inverse Fourier transform of this expression. We split 
into two parts, the first of them being 

1 := (2^72 / ^ cosc 2\k\te ik - x d 3 k (4.209) 



Thus 



(27r) 3 / 2 / = J E(y)j(x-y)d 3 y (4.210) 

where E is the inverse Fourier transform of cosc2|A;|t. Introducing spherical 
coordinates we again do the integration over the angles. We obtain 

E = ,^ 2 ^ 3/2 J k (e m(c2 ' + « + e lK(f - C2t) ) dn (4.211) 

Writing this as 

E = - c ^^3/2 ^ / (e tK(c2t+0 - e m «- C2 *)) dn (4.212) 



E{x) = - Q V2 j-d t [S(c 2 t + - - Pat)] (4.213) 



we obtain 

'7T\l/2 1 

This allows us to evaluate J: We can write the t-derivative in front of the 
convolution. We also omit the first Dirac delta, since (c 2 t + £) > 0. The 
remaining terms read 

= ±dt J ±8(\x -y\- c 2 t)j(y)d 3 y (4.214) 

We introduce spherical coordinates and use the Dirac delta to do the radial 
integration. This leads to 

J =^(i/ Kv)dSv) ( 4 - 215 ) 
4tt \c 2 t J dB ( x ,c 2 t) J 

To actually carry out the differentiation we transform the integral in one over 
&B(0,1): 

1 = IT 9 * ( * I j( x + C2ty)dSy] (4.216) 

47r V JdB(0,l) J 

Now we can interchange the integration with the differentiation, since the 
domain dB(0, 1) is independent of t. Then we transform the domain back. 
This finally gives 

4n(chf jdS " + 7if M dSy ) (4217) 

47T \C 2 t JdB(x,c 2 t) JdB(x,C2t) J 



4.3. LINEARIZED ELASTO-DYNAMICS 



85 



where denotes the outer normal derivative. 

Next we compute the inverse Fourier transform of the remaining terms 
in ()4.2()8|) .We start by computing 

U x := -1— [ C °^j fc|t (] ■ k)ke tk - x d 3 k (4.218) 



(2tt) 3 / 2 J |fcl 2 
Using the previous result we obtain 



(2nf' 2 II?(x) = - J -J^D^x - y)d ■ jd 3 y (4.219) 



where 



1 f cos Cl \k\t ik . . 



A(x) = j^j- 2 j —pjL-LWk (4.220) 



Along the meanwhile well known paths we convert this into a one-dimensional 
Fourier integral: 



-oo 

The integrals give (modulo constants) the sign distributions. Thus Di be- 
comes 

Dx(x) = J|l [9( Cl t + + 6(e - Cl t) - 1] (4.221) 
Consequently (since £ > we have 0(ci£ + £) = 1 and 5(£ + c x t) = 0) 

7T / a; _ ,. . a; 

5 (- ? ett-c l() + ? 



Sd(x) = ( -^e« - c,t) + ^<S(? - c,t) ) (4.222) 



This means that J/ := Jii — IJ2 is given by 

47r// = -/ (x-y)d-jds y + 

J B(x,at)\B(x,cat) \ x ~ y\ c l l JdB(x,dt) 

j (x-y)d-jdS y 

c 2 l JdB(x,c 2 t) 

We know from the previous case of initial data (|4.1fil|) how to deal with the 
volume integral. By partial integration we can split in up into a volume 
integral involving the undifferentiated initial data plus some surface terms. 
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Rearranging finally yields 



B 



L 



B(x,at)\B(x,c 2 t) \ x V 




-(3n A n B - 5 B )j A (y)d 3 y + (4.223) 



c l* JdB(x,cit) 




n B n A j A (y)dS y + 



+ — + [ [n A d A f-n B d Aj A ]dS y + 



C 2t JdB(x,c 2 t) 



+ — + [ n B d Aj A dS, 



c l* JdB(x,at) 



Then the full solution is the superposition of this solution with (|4.205J) . 

Note that properly speaking we have not proved the existence of a solution 
but have rather derived the form of the solution, if existence is assumed. In 
principle we would have to prove that the representation formulas we derived 
in this section really solve the system (j4.158j) and satisfy the initial conditions. 
But since this is a rather trivial task, once the formulas are given, we omit 
it at this stage. 



The minimum requirement a mathematical theory should fulfill to be viewed 
as physical reasonable is that the equations of motion describing the dynamics 
are well-posed in a certain sense. In our case the simplest non-linear setup is 
that of an infinitely extended elastic material on a given background space- 
time, since in that case we do not have to worry about effects arising from 
self-gravitation nor have to care about boundary terms. It is thus a crucial 
test for our matter model to provide this well-posedness for the field equations 
arising in the said setup. 

In this section we thus take a closer look at the equations of motion 
derived in the sect ions 14 . 1 . II and 14 . 1 . 21 We will compute the principal symbol, 
and check hyperbolicity of the equations of motion. Finally we apply the 
basic existence theorem from section 12.21 to obtain local well-posedness for 
the Cauchy problem for data close to those describing natural states of the 
material. 

Beside the results on well-posedness we will also derive some results on 
the underlying geometric structure of the equations. 



4.4 Unbounded setup 
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Throughout the following analysis we assume the internal energy to be a 
smooth function of its entries (i.e. of the configuration and the strain in the 
space-time description and of the body point and the strain in the material 
picture). 

We finally note that some of the results of this section can be found in 
0, or 0. 

Remark. Note that from now on we again deal with the fully non-linear 
case. 

4.4.1 Space-time description 

The equations of motion for the configuration took the form 

(-neh AB u»u v + AnU AMBN f M , a g a »f N , p g &v + %) d^dj 3 = R A (4.224) 

This is a quasi-linear second order system for the configuration f A . R A 
denoted the lower order terms and ^ AB was linear homogeneous in the stress. 

Second order well-posedness 

We will show that the principal part is regular hyperbolic in the neighborhood 
of a natural configuration. The existence of this natural configuration will 
be assumed in the following. 

To prove regular hyperbolicity for configurations close to the natural one 
(in a proper topology of course) it suffices to show that the natural configura- 
tion itself has the required property. By continuity (f.e. in the C 1 topology) 
the same is then true for all configurations sufficiently close. Remember that 
regular hyperbolicity implied weak hyperbolicity as was shown in section l2~2l 
The principal symbol for the natural configuration reads 

A AB (k) := -heh AB {u^f + AnU AMBN f M , a k a f N ,fs (4.225) 

This form of the principal symbol strongly suggest to try to prove regular 
hyperbolicity w.r.t the pair (u fl , —u v ). Indeed one finds that 

A AB {-u) = -nlh AB (4.226) 

which is negative definite. On the other hand for all / M satisfying = 
we have that 

A AB {1) = mU AMBN f M w r~f N , p f (4.227) 

which is positive definite by virtue of the Legendre-Hadamard condition im- 
plied in the definition of the natural configuration. Finally one checks the 
final requirement: w M (— it M ) = 1 > 0. Thus we have the following result: 
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Lemma 11. The system \4 . 224 ) ^ s regular hyperbolic w.r.t. the pair (ft M , 



-u. 



Then the corollary from section 12.21 immediately implies the following 
result. 



Theorem 3. The Cauchy problem for \4 . 224 ) ^ s well-posed in the following 



sense: Given a smooth natural configuration. Then there exists a foliation of 
space-time, such that the initial value problem for initial data (/o,/i) being 
close to the data induced by the natural configuration (Jo? A) in H s+1 x H s 
(for s > 3) give rise to a unique local (both in space and time ) solution, which 
depends continuously on the initial data. 

Remark. The problem with this result is that we only get existence and 
uniqueness locally in space. This can be improved if the we add an additional 
assumption. If the characteristic speeds of the material are smaller than the 
speed of light, one may obtain regular hyperbolicity w.r.t. (d t , di) all over the 
leafs of the foliation and not only in certain neighborhoods. For the case of 
an isotopic material this is explicitly carried out for the material setup in 
the following section in lemma XTb} but it would work just the same wy in the 
given context of the space-time description. 

Note that the connection between the Legendre-Hadamard condition and 
regular hyperbolicity is even deeper than we already know. The following 
holds 

Lemma 12. The system \4-224\l is regular hyperbolic w.r.t. (-u M , — u v ) for a 
reference configuration if and only if it is natural. If the material is isotropic 
w.r.t. the reference strain, then \4 . 224 ) ^ s regular hyperbolic if and only if 



the reference configuration is natural. 

The first part of the lemma is easy to prove: One part of the statement 
we have already checked. The other part is obvious from the structure of the 
principal part: suppose that the system is regular hyperbolic w.r.t — u v ). 
Then the Legendre-Hadamard condition for the elasticity operator follows 
directly from the definition. 

The second part of the lemma is a little bit harder. We will see that 
regular hyperbolicity w.r.t. any (X m ,clv) implies the Legendre-Hadamard 
condition. We have 

< A AB (k)m A m B < 4nU AMBN f M , a k a f N lf3 k^m A m B (4.228) 

for all k^X^ = 0. Since we assume isotropy, we can write this as 



< fim 2 K 2 + (/i + A)(m • k) 2 



(4.229) 
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where k a := f A , a k a and contractions are taken w.r.t. the strain. We can 
always find a k, which gives rise to a non-vanishing n. Then m = k includes 
2/j, + A > 0, while choosing m orthogonal to k gives rise to fi > 0. But those 
are just the conditions on the Lame moduli which establish the Legendre- 
Hadamard condition. The statement of the lemma is proven. 

Remark. This lemma does not generalize to non-isotropic materials. 

The reason for this lies in the following fact: suppose the system is regular 
hyperbolic w.r.t. a pair (X^, uo u ) where X M is not parallel to the four-velocity. 
Then the hyperbolicity conditions would imply the Legendre-Hadamard con- 
dition only on the subspace, which is the intersection between the orthogonal 
spaces of the four- velocity and the given X^. 

To recover the full Legendre-Hadamard condition one would in addition 
require < Uambn f^ia X a f N ,f3 X /3 m A m B to be true. But in general there is 
no way to obtain this inequality. One could of course try to derive a necessary 
condition on the elasticity operator which would result in the above inequality 
and thus would imply the statement of the lemma for such materials but this 
will not be done here. 

We make a further comment on possible choices for (X M ,av) for isotropic 
materials: 

Lemma 13. Let the material be isotropic. Then the system \4 ■ £%4\ ) ^ s regular 
hyperbolic w.r.t. the pair (X^,lo v ) if and only if the following inequalities 

(gT + (l ~ ^vTj < (4.230) 

as well as 

(g^ + (1 - c 2 ) UpuS) X»X V < (4.231) 

hold for c being the maximum of C\ = \J1\i + A and c 2 = y/JI while c is the 
minimum. For this lemma we also assume that 1 > c > c, i.e. that those 
constants are bounded from above by the speed of light ( the physical meaning 
is that the sound waves travel with physical velocities). uF" is the reference 
four-velocity. 

Remark. This lemma sheds some light on the geometric structure of the 
problem. It tells us that X^ as well as uj v have to lie inside certain sound- 
cones. Recall that we have found in the last section (at least for the linearized 
setup) that the elastic perturbations travel along these cones. 

To prove the lemma we first of all have to know the isotropic principal 
symbol. In principle (that means beside positive multiples) it is given by 

A(k) AB m A m B = -m 2 (u%) 2 + (fim 2 * 2 + (// + A)(m • k) 2 ) (4.232) 
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where k a = f A ,^ k^ and contractions are taken w.r.t. the reference strain. 
Without restriction of generality we assume m 2 = 1. 

Recall the definition: a system is regular hyperbolic w.r.t. (X 11 , uj v ) if and 
only if A{uj)ab is negative definite, while A(1)ab has to be positive definite for 
all l^X^ = 0. By continuity small changes in the argument k do not change 
the positivity properties of A{k)AB- A change of these positivity properties 
would reflect itselve in the condition det A{k)AB = 0. 

In general this determinant can be computed to be (see f.e. 

det A(k) AB = H 3 (sTW(^W ( 4 -233) 

where 

9? = 9» u + (l - ^) v?u v (4.234) 

(later on we will use g^ u and with the same meaning) The determinant 
vanishes if and only if g^k^ky = or g^k^ky = 0. Since the gf are non- 
degenerate these conditions define cones in the co-tangent space. 
The principal part taken w.r.t. — is negative definite. Since 

(gT + - u^u u = -i < (4.235) 

the principal part is negative definite w.r.t. all uj v satisfying the condition 
g fJ,v uj f j i uj u < 0. Thus all co- vectors lying inside cone g^k^ky = are sub- 
characteristic. 

The proof for is slightly more difficult: We need the following result 

Proposition 1. g^k^ky > for all k^X^ = if and only if g~yX^X u < 

This proposition however is trivially true since g^ v is a Lorentz-metric for 
1>£ 2 . 

The proposition gives the second inequality of the lemma. At first glance 
it might seem strange that one inequality involves c while the other one 
contains c. This comes from the fact, that in one case we approach the 
co-sound-cone from above and in the other case from below. 

Note that (strictly speaking) up to now we have only proven the "if '-part 
of the lemma. To complete the proof we simply note that positive (respec- 
tively negative) definite operators have positive (respectively negative) de- 
terminant. This observation gives the "only if '-part and thus completes the 
proof of the lemma. 

Remark. Note that the inner and outer cone interchange when viewed in 
tangent and co-tangent space respectively. Thus X^ lies in the inner cone, 
while u u lies in the inner co-cone. 
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Rank two vs rank one positivity 

Usually people in the relativity community try to prove well-posedness by 
casting the equations of motion into a symmetric hyperbolic first order sys- 
tem. We will not give a definition (which can e.g. be found in [2HD- Just 
note that there exists a complete local-in space-theory for such systems. We 
could of course try to apply this theory for the given setup as is done f.e. 
in But in general there may be a prize to pay for the transition to the 
first order system. In the following we shall briefly sketch, how this problem 
looks like and how it (at least in the special case of an isotropic material) 
can be overcome: 

If one prefers to cast the original equations into a first order system, one 
would need to contract the principal part not only with rank one objects but 
with general tensors with two indices to obtain hyperbolicity and hence well- 
posedness. Hence the original Legendre-Hadamard condition may in general 
fail to be sufficient, since it only covers rank one positivity. 

Note that the principal symbol and even the original second order system 
is unchanged if we exchange the elasticity tensor by 

Uabcd — ► Uabcd '■= Uabcd + J^abcd (4.236) 
if the additional term has the following symmetries: 

h-ABCD = h-ADCB = ~^CBAD = ^CDAB (4.237) 

We will see (at least in the case of an isotropic material) that a gauge transfor- 
mation of the kind introduced in ()4.236j) can be used to obtain hyperbolicity 
of the first order system by just using the Legendre-Hadamard condition. 

For an isotropic material we have seen that the elasticity operator can be 
written as (omitting the upper tilde for the reference strain) 

^Uabcd = e (2/jh A {ch D )B + A/iab^cd) (4.238) 

Remark. For convenience we will omit the factor 4 on the left hand side as 
well as the stored energy e from now on for the rest of this subsection. 

Since we do not assume the material to be homogeneous they need not 
be constant and may still be functions on the body. Hab is the inverse strain 
at the reference deformation. On flat background we may use coordinates 
such that Kab = $ab but we will stick to this more general notation. We 
know that the Legendre-Hadamard condition is equivalent to the following 
condition on the Lame moduli (compare section 13. 1.3j) : 



fi>0 2(i + A > 



(4.239) 
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Our goal is to establish a stronger positivity result for the transformed elas- 
ticity tensor Uabcd using these conditions on the Lame moduli. Namely we 
will show positivity of 

Uabcd^ AB ^ cd = U ABCD ^ AB <f D + A ABCD ^ AB <f D (4.240) 

Splitting into its antisymmetric {uj ab ) and symmetric but trace- free (x AB ) 
parts we obtain 

= U AB + X AB + 1 tr ^AB ( AM \) 

3 

where tr <fi = h AB <p AB is the trace of (j) w.r.t. the reference strain. The first 
term in (|4. 240)1 can then be written as 

y + A)(tr0) 2 + 2/i X CBX CB (4.242) 

The indices were lowered using the reference strain. This is nonnegative if 
and only if fx > and the bulk modulus k := 3 A + 2/x > 0. 

Remark. We note two things: 

• those conditions are stronger than \J^.239^ 

• they do not include the antisymmetric part of cj) AB 

We cannot use this expression to bound the contracted original elasticity 
tensor away from zero, since the left hand side vanishes for antisymmetric 
(f) AB . Thus the gauge term Aabcd must give rise to antisymmetric terms 
thus breaking the symmetry (isotropy) of the setup. 

The easiest and at the same time most natural choice for a gauge term is 

J^abcd = ah A[B h D]c (4.243) 

This term obeys the required symmetries as stated in ()4.237j) . The free 
function a will be specified in the following. Computing the second term in 
f!4.240p using this gauge term we obtain 

A ABCD <^ AB <f D = a 0(tr 0) 2 - X AB Xab + uJ AB u AB ^j (4.244) 

Adding this to the terms from the original elasticity tensor gives 

U{<f>, 0) = Q0u + a) + A J (tr 0) 2 + (2/i - a) X AB XAB + au AB u AB (4.245) 
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If we would chose a = 2/i then we could use the Legendre-Hadamard condi- 
tion to bound the first and the last term but the second term would vanish. 
This however would lead to problems as indicated in the remark above. We 
thus need a more refined argument. Assume the Legendre-Hadamard condi- 
tion holds. Then by continuity there exists a 5 > such that 

2/i + X - ^5 > (4.246) 
2/i - S > (4.247) 
Then we set a = 2/i — 5 and obtain 

U{<f>, 0) = f 2/i + A - (tr 0) 2 + 5 X AB Xab + (2/i - 5)cu AB u AB (4.248) 

Now all coefficients are positive by virtue of the Legendre-Hadamard condi- 
tion and all summands are nonnegative. 

We have shown how to construct a gauge term which leaves the second 
order system unchanged but allows for the Legendre-Hadamard condition to 
be sufficient for positivity and hence for hyperbolicity in the context of the 
first order system. 



4.4.2 Material description 

Now we turn to the material setup. In this section we will show that the 
equations of motion derived from the material Lagrangian are regular hyper- 
bolic close to a natural reference state. This will include well-posedness for 
the second order system. 

We could also try to cast the system into a symmetric hyperbolic system. 
Then in general stronger assumptions on the elasticity operator (and thus on 
the material) would be needed. For isotropic materials the last part of the 
previous section could be used to overcome this discrepancy. But (besides 
this remark) we will not go any deeper into the first order formulation, since 
the second order formulation works out well. 

As in the space-time setting we will find that there exists a preferred 
pair (X^,u u ), with respect to which the systems unfolds its hyperbolicity 
properties in the most natural way. Finally the connection between this very 
choice and the pair (« M , — u u ) used in the previous sections will be discussed. 

The equations of motion were derived in the form (see section 14.1.21 for 
the derivation) 

Agd&Fi = Ri (4.249) 

This is a second order quasi-linear system for the deformation F J . As before 
Ri denotes lower order terms. 
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Remark. Note that Greek indices are not space-time indices but rather de- 
note points ofM. x B. 

It turns out to be most convenient not to use the principal part directly 
to check the various positivity requirements but rather work with 

B%:=A^F\ A F\ B (4.250) 

Since F l , A is an isomorphism, this newly defined B^ B has the very same 
positivity properties as the original A^ had. The advantage is that the new 
object can be given in the more convenient form (see again section I4.1.2J) 

Kb = r k (-eh AB U»U» + AUacbd^^) + ^ B (4.251) 

where $^ is linear homogeneous in the stress. We have seen in section 
14.1.21 how to obtain and <\> Aii for Gaussian coordinates on space-time. It 
is straightforward to see that for general coordinates (i.e. for non-vanishing 
lapse and shift) we obtain very similar expressions 

U» = jV -1 7^o (4.252) 
<f» = N'H^W a f, a +8%h CM (4.253) 

For the rest of this section we will stick to this more general formulation. Note 
that 7 is thus again given by the more general expression 7 = (1 — W 2 )^ 1 
Then for a reference deformation the principal part reads 

^B% B = -eh AB UW + AU ACB D~f^ Du (4.254) 

We now try to find a pair (X^av) w.r.t. which the system is regular hy- 
perbolic. Note that the principal part consists of two portions. The first is 
non-positive, while the second is non-negative once the Legendre-Hadamard 
condition is assumed for the elasticity operator. Thus a good choice for each 
element of the pair (X M , uj v ) will annihilate one part of the principal symbol. 
For the reference deformation this is in general given by 

TB(k) AB = -eh AB (U^) 2 + AU ACBD ^ c ^ D "k u (4.255) 

From the form of the principal symbol it is clear that the choice for X M has 
to be 

X^ = (4.256) 
because then for all LX M = we have that 



7*B{1)ab = 4Uacbd4> C %0 Du I, 



(4.257) 
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which is positive by the Legendre-Hadamard condition unless c^fH^ = 0. 
This however can be excluded since 

4°% = h CM l M ^ (4.258) 

which is non-zero by the non-degeneracy of the strain. 

In order to obtain a negative principal part for uj v we have to assure that 
the part involving the elasticity operator does not enter the symbol. This is 
guaranteed only by the following condition: 

^ Ufl = o (4.259) 

We conclude 

com = -N- x u)Ji MC W a f G , a (4.260) 
Thus for fixing Uq through the (arbitrary) condition X^u^ = 1 we obtain 

u = NT h (dt - N- l h MC W a f c , a dX M ^j (4.261) 
Together this gives the following result 

Lemma 14. The system \4-24ty is regular hyperbolic w.r.t. the pair 

(N- 1 ^d t ,NTkdt-N- 1 h MC W a f c , a dX M )) (4.262) 
in a neighborhood of a natural deformation. 

Remark. Needless to say that the system is also regular hyperbolic w.r.t. 
all pairs sufficiently close to the one given above. Nevertheless the system 
\4-249jj satisfies the conditions for regular hyperbolicity in this optimal way 
only for the above choice. 

We can thus obtain the analogous result as in the previous section stating 
that the system ()4.249j) is regular hyperbolic near a natural deformation if 
and only if it is regular hyperbolic w.r.t. the above choices for X^ 1 and u v . 
The form of the corresponding lemmas can be taken word by word from the 
last section. Instead of spending too much time on this, we shall rather clarify 
the connection between the above choice for (X^,oj v ) and the four- velocity. 

For this purpose we make use of the transformation introduced in section 
I4.1.2[ namely the diffeomorphism 



(t, X) = (t, f{t, x)) (t, x) = (t, F(t, X)) (4.263) 
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We have seen in the section mentioned above, that the material velocity maps 
to which is a multiple of dt- We have also seen that 

f\aT^4^ (4.264) 

From this one may guess that — it M i— > a; M . To see this we use the following 
result from section IH. 2. HI 

(a ) - (G )-( ~ N2{1 ~ N ^ W ' F3 ^ \ U 265) 

From this it is easy to obtain that 

u v i- U»G, U = N-^Go, = -N-y-hl + ^g l3 W l F\ A ^ (4.266) 

To complete the argument we have to show that the spatial parts coincide 
(The time-part already fits). This can be seen by recalling the formula for 
the inverse strain derived in (jH.lOlj) : 

h MC f C , a W a = ( 9ij + jWiW^F^M W j = -yW^M (4.267) 

Using this calculation on the spatial art of u we immediately find that w M t— > 
— thus proving the claim. 

Remark. The above computation shows how the geometric importance of the 
four-velocity carries over to the material description. 

From our observations up to now the following lemma connecting the 
material and the space-time description is immediate: 

Lemma 15. The material system close to a natural deformation is regular 
hyperbolic w.r.t. (X^,u u ) as defined above if and only if the system in the 
space-time picture is regular hyperbolic w.r.t. (-u M , — u u ) where is the natu- 
ral normalized four-velocity . Moreover (X^jLUp) and (£t M , —u v ) are connected 
via the diffeomorphism introduced in section ^. 2. as indicated above. 

Remark. To prove this one just has to use that from regular hyperbolicity 
w.r.t. one of the above pairs the Legendre-Hadamard condition follows, which 
gives the hyperbolicity in the other description. 

We carry on with a local existence result. As in the space-time descrip- 
tion we can use the basic theorem given in section 12.21 to obtain local well- 
posedness: 
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Theorem 4. The Cauchy problem for \J^.2J^9^ is well-posed in the following 
sense: Given a natural deformation. Then there exists a (local) foliation of 
R x B, such that for initial data (Fq, F() close to the ones induced by the 
natural deformation on the leafs of the foliation in the H s+1 x H s topology 
(s > 3) there exists a unique regular (i.e. at least C 2 ) local (in space and 
time) solution which depends continuously on the initial data. 

Remark. The actual function spaces involved can be taken from the original 
theorem from section \2~2\ 

Remark. The reason why we only get local in space solutions in general lies 
in the fact that ( and hence u^) is not closed in general. 

There is another natural choice for the pair (X^, u u ), given by the foliation 
R x B. This pair is given by (d t , dt). In the following lemma we will see how 
to deal with this setup (at least for isotropic materials) 

Lemma 16. Let the material be isotropic w.r.t a natural deformation. Then 
the system \J^.2J^9^ evaluated in a neighborhood of the natural deformation is 
regular hyperbolic w.r.t (d t ,dt) if and only if the characteristic speeds c\ := 
2/x+A and c 2 , = t 1 are bounded from above in the following way: W~ 2 > c\ > 
as well as W~ 2 > c\ > 0. 

The following corollary follows immediately from the 1 > W 2 > 0: 

Corollary 3. The isotropic system is regular hyperbolic w.r.t. (d t ,dt) close 
to the natural deformation if the characteristic speeds Ci are bounded by the 
speed of light, i.e. 1 > Cj > 0. 

Remark. It is interesting that due to the lemma one may obtain regular 
hyperbolicity w.r.t. (d t ,dt) even if ci > 1. 

Remark. An analogous result can be obtained for the space-time setup for a 
suitable foliation. 

For the proof of the lemma we note that the Legendre-Hadamard con- 
dition guarantees that the principal symbol is positive definite w.r.t. all Z M 
satisfying Iq = 0. This has already been mentioned in the proof of lemma 
Full The new part consists in showing that dt is sub-characteristic. We have 
to show that the principal symbol (for convenience we omit the tilde) taken 
w.r.t. dt is negative definite, i.e. 

N 2 ^B(k) AB = -eh ABl + AU AC BDf C ,c W c f D , d W d < (4.268) 
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Since the coefficient on the left hand side is positive (recall that 1 > 7 > 0) 
we only have to care about the right hand side. Using the expression for the 
isotropic elasticity operator, the question comes down to 

-m 2 7 + ((/i + A)(m • wf + fim 2 w 2 ) < (4.269) 

for non-vanishing m A where products were taken w.r.t. the strain Hae, 
while w A := f A , a W a . The coefficient in front of (m ■ w) 2 lacks a clear 
sign (the Legendre-Hadamard condition only leads to the conditions given 
in the lemma). We thus first have to rewrite the above expression. This is 
done by splitting m A into a part parallel to w A and an orthogonal rest: 

m A = rn 1 W wA + m A (4270) 
w 2 

This leads to 

m 2 = (m 7 )2 +ml (4.271) 
w 2 

Thus besides positive multiples the principal symbol twice contracted with 
m A is given through 

ml + (m " W H 7 + (2fi + A) (m ■ w) 2 + ^m\w 2 (4.272) 
//- / 

now the Legendre-Hadamard condition can be used to guarantee positivity 
of the coefficients. Since m\ and (m • w) 2 can be chosen independently, the 
symbol is negative definite if and only if the following two inequalities hold: 

- m 2 ± -f + ixm\w 2 < (4.273) 
-{m ■ w) 2 7 + (2/i + A)(m • w) 2 w 2 < (4.274) 

Assuming the non-trivial case, those two inequalities are essentially of the 
type 

-7 + c 2 w 2 < (4.275) 

It remains to compute w 2 = W a W b f A , a f B ,bhAB and compare the result 
to 7. From the formula for the inverse strain given by equation (jH.lOlj) we 
obtain 

w 2 = 7H/ 2 (4.276) 
Thus the inequality we need becomes 

-1 + c 2 W 2 < (4.277) 
which is the claimed inequality. 

Remark. Note that we can use the lemma to establish a local existence result 
on all of B getting rid of the locality in space. The same argument works for 
the space-time description. There we can use it to obtain valid data on all of 
t = const. 
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4.5 Gravito-elastodynamics 

Once the initial value problem for the unbounded setup on a fixed background 
is shown to be well-posed, there are two natural ways to proceed. One 
can consider finite material bodies or one can investigate unbounded self- 
gravitating elastic materials. While the former task will be undertaken in 
the following section 14.61 we shall concentrate here on the latter challenge. 
We will show that the initial value problem for a self-gravitating unbounded 
elastic medium is well-posed under certain (mild) restrictions on the material, 
already well-known to us from the previous sections (namely the Legendre- 
Hadamard condition). 

We will first discuss the initial value problem for the Einstein field equa- 
tions. Following the exposition of Friedrich and Rendall {21J we will use 
the harmonic gauge to cast them into a quasi-linear wave equation, which is 
then shown to be regular hyperbolic. Then we will discuss issues like gauge 
propagation and proper initial data. 

Then in the second part we will show that the coupled system consisting 
of the Einstein (-matter field) equations (in harmonic gauge) and the elas- 
tic equations of motion in space-time description is regular hyperbolic and 
discuss how the arguments given for the un-coupled case carry over. 



4.5.1 Einstein field equations 

Let space-time be denoted by (M, g^)- The Einstein field equations de- 
termine the space-time geometry (encoded in the metric g^ v ) through the 
space-time entities (matter). They are obtained from the action principle (L 
denotes the matter Lagrangian) 

/ (Ry/-det g^ + 2nL)d 4 x (4.278) 

J M 

and read 

G>„ = kT^ (4.279) 

where G^ u is the Einstein tensor formed from the Ricci tensor R^ u , which 
depends on the metric and its first and second derivatives while T^ v is the 
energy-momentum tensor while the constant k = (G being the gravita- 
tional constant). An equivalent version reads 

= k [T^ - l -T a pg af3 g,A =: Mp, (4.280) 
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For vacuum the energy-momentum tensor vanishes and the Einstein field 
equations simplify to 

R^ = Q (4.281) 

Remark. For the moment we do not specify a particular matter model. In 
this spirit we do not give an explicit expression for the energy-momentum 
tensor. We rather assume it to be a functional of the matter fields to keep 
this discussion on a general level. Later on it will of course be the one derived 
for elastic matter \3. 6<ty . It is however important to assume that it does not 
contribute to the principal part. 

When written in terms of the metric equation f)4.280|) reads 
-\g Xp 9,u M + VfolV) + Tl^g^Yl + = M, v (4.282) 

where the contracted Christoffel symbols are defined by T M := ■ Since 

the Christoffel symbols involve first derivatives of the metric, the principal 
part is given by 

P^ s d,d s g a3 := \g xp g^x P + W) (4.283) 

Remark. We could exchange the covariant derivative by an ordinary one, 
because the difference is of lower order in the metric. 

Note that this is a quasi-linear expression, since the inverse metric de- 
pends on the metric itself but not on higher order terms. One can show (see 
e.g. [21]) that every hyper-surface is characteristic in the sense that for every 
choice of fc M the characteristic polynomial (see chapter F2.il for the definitions) 
vanishes: 

P(Jfe) = (4.284) 

This in particular implies that there exists neither a sub-characteristic co- 
vector nor can the equation in the form ()4.282|) be hyperbolic in the sense of 
any of our definitions. This defect can however be shown to be a result of the 
diffeomorphism invariance (i.e. of coordinate freedom) of the field equations. 

In order to obtain a reasonable system of differential equations we have to 
make a suitable choice of coordinates. The form of the principal part ()4.283|) 
suggest that a good coordinate system would allow us to get rid of the term 
dfjjT^, since we then would be left with a wave-equation. This is possible, 
since the Christoffel symbols are no tensorial objects. We can thus hope to 
find coordinates for which 

T M = (4.285) 
More precisely we have the following: 
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Lemma 17. Given a space-like hyper- surface. Let x a be coordinates in some 
open neighborhood U in this hyper- surface. Then there exists coordinates y^ 
on some neighborhood ofU in M with the property (y^)\u = (0, x 1 ) such that 
in these new coordinates 

P 4 = (4.286) 

Definition. The coordinates introduced in the fore-going lemma are called 
"harmonic coordinates" . 

We can prove the lemma by rewriting ()4.286j) (we construct harmonic 
coordinates): 

U x y*{x) = U y y» = -V = (4.287) 

where the second Christoffel symbol corresponds to the new coordinate sys- 
tem (y M ) . The first equality holds since coordinates are scalars on space-time. 
Using proper initial conditions one can obtain the solution claimed in the 
lemma (a way of seeing this is by application of the theorem in section l2~2"|) . 

In harmonic coordinates (which will from now on be denoted by x^) the 
principal part ()4.283|) becomes 

Af^d.dsg^ := \g xp g^x P (4.288) 

At the given stage it is easy to show that the principal part (after raising the 
indices {ijlv) in the above equation using the metric) in harmonic coordinates 
is regular hyperbolic: 

Lemma 18. The principal part \^ . 28 8\) is regular hyperbolic w.r.t. every pair 
(X^,u u ) where both X M and uj u are time-like and satisfy X^uj^ > 0. 

Given proper initial data we can thus apply the existence theorem from 
section 12.21 to obtain local well-posedness for the Einstein equations in har- 
monic gauge (assuming proper well-posedness of the matter equations). 

Remark. The procedure of reducing the Einstein field equations to a hyper- 
bolic system by a suitable choice of coordinates is called "hyperbolic reduc- 
tion". 

The story however does not end at this point. Two (three) questions 
remain open: 

• What are proper initial data? The related question: Are there any 
constraints or can the data be given freely? If there are constraints, do 
they propagate? 

• Does the gauge condition hold off the initial surface? (Is it preserved 
under the evolution?) 

These question will be answered in separate subsections. 
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Initial data and constraints 

One may at first guess that the metric and a suitable transversal derivative 
given on a suitable initial surface play the role of proper initial data. This 
however is not the case. We will see that an initial data set consist of a three- 
dimensional Riemannian manifold (S, gij), which in addition is equipped with 
a symmetric covariant two-tensor fcy plus a diffeomorphism identifying S 
with a space-like hyper-surface in M = R x S. Then plays the role of first 
fundamental form (induced metric) on this initial hyper-surface, while is 
the second fundamental form (external curvature). If matter is present, the 
initial data set has to be extended to include the matter data. 

While lapse and shift can be prescribed freely (each choice corresponds 
to a certain coordinate choice on M), and kij are subjected to certain 
constraint equations. From these "constrained" initial data one can construct 
initial data for (I4.282jl . 

To actually see this we have to carry out a 3+1 decomposition of the 
Einstein field equations ([4.280)1 . Assuming a foliation of space-time by space- 
like hyper-surfaces {t = const} we write the metric in its ADM representation 

g^dx»dx u = -N 2 dt 2 + g^dx 1 + Y i dt)(dx j + Y j dt) (4.289) 

Then the projector onto the tangent space of the hyper-surface {t = const} is 
= g^ u + n^n u where the unit-normal n^ is given by n M = N^ 1 (5q — <5f Y l ). 
Then g^ is the induced metric on {t = const}, while the second funda- 
mental form is given by 

h, = ^{d t -C Y ) gij (4.290) 

One finds that certain components of the Einstein equations (j4.279J) beside 
material quantities only depend on the first and second fundamental form 
and do not involve higher order transversal derivatives of the metric. First 
note the following trivial consequences of (|4.279|) : 

= n *n v (G^ - kT» v ) (4.291) 
= n^h\ (G M „ - kT^) (4.292) 

The interesting point about these equations is that they can be written as 

= r - k ab k ab + k 2 - 2kj (4.293) 

= 2D [a k a b] + K 3b (4.294) 

where 

T^n^ =: j T^n^h va = -j a k = k ij9 ij (4.295) 
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Here r denotes the Ricci scalar corresponding to the induced metric (one 
must use the Gauss and the Codazzi equation to arrive at the stated result. 
See [IB] for details) while D a is the covariant derivative associated with g^. 

Note that equations (j4.293|4.294j) form an under-determined system for 
first and second fundamental form, once the material quantities are known. 
In this sense (j4.293p4.294j) constitute necessary restrictions on possible initial 
data and are therefore called constraint equations. If the coupled system is 
considered, these equations also constrain the material initial data. 

It remains to show how to construct the initial data g^ and dtg^u for 
()4.282|) on {t = 0} once (g^, kij) (or rather the corresponding objects on the 
initial hyper-surface) are given. 

From the 3+1 decomposition of the Christoffel symbols one can derive 
the following equations for lapse and shift (see for the derivation) 

d t N — N,i Y i = AT 2 (fc-n M r^ (4.296) 
d t Y i -Y i , j Y j = N 2 {f - DMniV- W V Y V ) (4.297) 

Here 7* denote the contracted Christoffel symbols corresponding to the in- 
duced metric. 

Given lapse and shift and assuming harmonic gauge, the above equations 
give the time derivatives of lapse and shift in terms of (gij, k^, N, Y l ). What 
we still lack is an equation for d t g%j in terms of the data. For this we use the 
definition of the second fundamental form ()4.290|) . We can rewrite it in the 
form (C v denotes the Lie-derivative w.r.t. the vector field v l ) 

d t9ij = 2Nk ij + C Y9ij (4.298) 

This equation gives the sought time-derivative of the induced metric. 

We sum up: Given initial data (S, g^, k^) (we suppress the material data 
for the moment) which obey the constraints ()4. 29314. 294|) . These data induce 
corresponding objects on the initial hyper- surface, which is the isomorphic 
embedding of (S, g^) into space-time. Then for any sufficiently smooth choice 
of lapse and shift (this choice is pretty arbitrary, it just fixes the time flow off 
the initial surface) we can construct g^ u on the initial surface. The compo- 
nents of d t g^ v can be obtained using equations (j4.296l4.297l4.298l) under the 
restriction = (coming from the requirement for harmonic coordinates) 
as indicated above. 

We have seen that the relevant initial data are given by (S,gij,kij) to- 
gether with an isometric embedding into space-time. It has been shown how 
these data can be used to obtain suitable data for the problem as formulated 
in (JOBS). 
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Propagation of gauge 

Assume we have proper initial data on a hyper-surface, such that the theo- 
rem from section 12.21 provides us with local existence and uniqueness. (The 
accurate form of this statement will be presented in the last section.) Then 
we still have to show that this solution (obtained from the reduced system) 
actually is a solution to the full Einstein field equation ()4.28U|) . This is done 
by showing that if the contracted Christoffel symbols vanish initially, then 
they also vanish off the initial surface. 

For the solution of the reduced system g^ u we construct the un-contracted 
and contracted Christoffel symbols (which for the moment we do not assume 
to vanish). For these symbols we derive the Einstein tensor G^ u . Subtracting 
the energy-momentum tensor and the reduced field equations (which hold due 
to the assumptions) leads to (r M is assumed to be computed from g^„) 

G^ - kT^ v = G^ - kT^ v - (reduced equations) = V^T v) - ^g^VaT 01 

(The second term arises when we express g al3 g l _ iu ,ai3 in terms of covariant 
derivatives.) Now we employ the contracted Bianchi identity (see f.e. jSHJ 
or |Hj), which hold for purely geometrical reasons independently of the field 
equations: 

VG,, = (4.299) 

If the energy-momentum tensor T^ is divergence-free (which is the case for 
elastic materials according to lemma EJ) this leads to a homogeneous wave 
equation for the contracted Christoffel symbols: 

□r„ + R^Tu = (4.300) 

If we can show that and dtT^ both vanish on the initial surface, we may 
employ the basic existence and uniqueness result from section 12.21 obviously 
T M = is a solution for trivial initial data. The uniqueness part then assures 
that this is indeed the only solution. This means that g^ u found from this 
harmonic coordinate ansatz is a solution to the Einstein equations. 

To see that the initial data for ()4.3()()j) are indeed trivial, we recall equa- 
tions ()4.296|4.297|) . We used the harmonic gauge to obtain d t N and d t Y l 
in terms of the given data. Reverting this argument implies that the initial 
data have been chosen in a way such that V ^ vanishes initially. 

The argument for d t T ^ works in a similar way, but is more complicated, 
since the expression under consideration is second order in the metric. By 
this it is not sufficient anymore to argue in terms of initial data only. We 
have to make use of the field equations. These enter in the form introduced 
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in the beginning of this sub-section. Contracting with the unit-normal n M 
leads to 

n v (Gp, - kT^) = n v (v^T u) - ^g^jA (4.301) 

Recall that the left hand side only involves initial data when evaluated on 
the initial surface (this was shown in the last sub-section). But the initial 
data satisfy the constraint equations, so that we can conclude 

= n v V {M T u) - irvV a r Q (4.302) 

Since we already know that T M vanishes on the initial surface this can be 
written as (we stress that this equality only holds on the initial surface as do 
the following equations. We include no statements off the initial surface.) 

= n v d^V v) - ^n lt d a T a (4.303) 

Note that we can omit spatial derivatives of the contracted Christoffel sym- 
bols, since the latter vanish identically on all of t — 0. It is then a simple 
matter of inserting the definitions of the unit-normal together with the ADM- 
representation of the metric to obtain the desired result. We can conclude 
that dtY n has to vanish on the initial surface. 

Thus we have found that the data for ()4.300j) are trivial so that the only 
solution is trivial as well. We conclude that the metric obtained from the 
reduced equations is indeed a solution to the Einstein equations. 

Remark. Note that the only restriction on the matter model was the divergence- 
less-ness of the energy-momentum tensor. For elastic materials which satisfy 
the equations of motion this condition holds. This is true independently of 
the equations of motion for the metric. 

Propagation of constraints 

We have seen in the previous section that the constraints (j4.293l4.294j) nat- 
urally arise in the context of 3+1 decomposition. If we would have used the 
Einstein equations in their ADM-version (i.e. as evolution equation for the 
induced metric), we would have to assure that the constraints hold off the 
initial surface. 

In the given context however this requires no further work, since a so- 
lution to the initial value problem connected to the equation ()4.282|) gives 
rise to a solution g^ u of the full Einstein equations ()4.280j) . Thus at each 
instant of time the constraints are trivially satisfied and require no further 
argumentation. 
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This property still holds once matter fields (and the corresponding equa- 
tions of motion) are added to the system. 

Hyperbolicity for the ADM variables in harmonic gauge 

Although this is somehow "off- route" we will show how the harmonic gauge 
can be used to obtain a well-posed system for the ADM- variables. Therefore 
we first define 

g^^:Z^(N,Y k , gij ) (4.304) 



and introduce 



In this notation we have 
This gives for the derivatives 



> A ):=(N,Y\g lk ) (4.305) 
9^ = (4.306) 



<W = Z, vA <p A p (4.307) 



and 

W = V,W>? + Z, vA <p A pT (4.308) 
so the new principal part derived from ()4.288|) is given by 

Z pT Z, uA< p A pT . (4.309) 

Note that the mapping Z : (fi A i— > Z^^) is a diffeomorphism, so that Z^.a 
is a linear isomorphism. 

The symbol corresponding to ()4.309|) reads 

k 2 Z,, uA (4.310) 

In order to apply theoreml we have to modify this. The most natural way 
is by use of Z^x. For this purpose we first raise the free Greek indices in 
(j6.74j) by a non-degenerated contra- variant two tensor , which we will be 
specify later. It may depend on the ADM variables and their first but not 
on their second derivatives. The new symbol reads 

k 2 Z^, B G^'G^'Z^ A (4.311) 

In order to apply our definition of regular hyperbolicity we have to modify 
this. The most natural way is by use of Z^a- For this purpose we first 
raise the free Greek indices in (16.74)1 by a non-degenerated contra-variant 
two tensor G MM , which we will be specify later. It may depend on the ADM 
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variables and their first but not on their second derivatives. The new symbol 
reads 

k 2 Z^ u ,, B G^'G uu 'Z^ A (4.312) 

The proper choice of G^ v is crucial for the whole argument to work. A natural 
candidate is the following: 

G m> = z »v + 2n »n v (4.313) 

where = N~ 1 (—5q + Y i 8^) is the unit normal to the t = const hyper- 
surface. If we contract 1)4.312)1 with S(p A S(j) B , we obtain 

k 2 N~ 4 (4N 2 (6N) 2 + 2g ij 5Y i 6Yi + 5g lk 6g jig ' k g a ) (4.314) 

Since the terms inside the bracket is positive, the symbol is positive (negative) 
definite for space-like (time- like) choice of k^. Therefore the symbol is regular 
hyperbolic with respect to every pair (X m ,cj„) where both and u v are 
time-like and satisfy X^u^ > 0. 

4.5.2 Coupled 

Here we will collect the partial results from the fore-going sections to add 
them up into a local well-posedness result for the equations of gravito-elasto- 
dynamics. We know that the Einstein equations in harmonic coordinates 
form a regular hyperbolic second order system. We also know that the solu- 
tion of this system is indeed a solution to the Einstein field equations 1)4.280)1 
if the energy-momentum tensor is divergence-free (which in our case is trivial 
since the latter condition is equivalent to the matter equations of motion, see 
lemma EJ)- Finally we know which restrictions (namely the constraints) have 
to be imposed on the initial data and that these constraints propagate in the 
proper sense. 

We know from section 14.4.11 that equations governing the dynamics of 
an elastic material on a fixed background are regular hyperbolic close to a 
natural reference state (i.e. a stress-free state, which satisfies the Legendre- 
Hadamard condition). 

What we lack up to now is knowledge on the coupled system. For this 
purpose we take a closer look at the action describing the coupled scenario. 
For a self-gravitating elastic materials we rewrite f)4.278|) as (one could include 
a coupling constant in front of one of the two terms, but we set it equal to 
one, since we can f.e. absorb it into the stored energy) 



( (R + ne) a/— det g^d^x 

M 



(4.315) 
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At this stage we point out the following important detail: the elastic La- 
grangian does not depend on derivatives of the metric while the Ricci scalar 
is independent of the configuration at all. Thus the principal part decouples 
and in a harmonic coordinate system we find that we can write the system 
as 

9 T- a Z ) a » 9 - ( 7 ) - lot - (4316) 



where 

d 2 ne 

is the principal part for the material equations as given in ()4.30|) (We have 
divided the material principal part by J— det <7 Mi ,). As already indicated 
previous to lemma ITH1 it is more convenient to raise the indices in the portion 
of the principal part involving the second order derivatives of the metric. 
(There appear extra terms, but they are omitted here because they are at 
most of first order.) The new principal part reads 



9 °\'° p i u \ a a I 9 



W \j (4-318) 



Afj J " v \ / 

Next we recall the material principal part for a reference configuration (re- 
call that all quantities corresponding to the reference configuration f A were 
labeled by an upper tilde; see ()4.37|) for the principal part) 

AZ = -nehAB^u" + 4nU AMBN f M , a g a »f N ,p (4.319) 

We know from the section 14.4.11 that this is regular hyperbolic w.r.t the pair 
(m m , —u v ). Since the four- velocity m m is time-like by construction, we find 
that the principal part 1)4.318)1 evaluated at a natural configuration is regular 
hyperbolic w.r.t. (m^, — u v ). This is due to the diagonal form of the principal 
part which itself is a result of the uncoupling of the equations in the highest 
order. By continuity (the stored energy is assumed to be sufficiently smooth) 
the same is then true for all configurations sufficiently close. We thus have 
the following statement: 



Lemma 19. The principal part \^ . 31 is regular hyperbolic w.r.t. the pair 
figuration. 



[u^, —u^ for all pairs (g^, f A ) where f A is close enough to a natural con- 



Remark. Note that the definition of natural configuration did involve the 
Legendre-Hadamard condition. 
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We can thus apply the existence theorem (or the corollary, depending on 
whether is closed) from section l2~2l to obtain local well-posedness. 

There remains just one small obstacle regarding proper initial data: we 
have to identify the material quantities used in the constraint equations 
()4.293I4.294|) in terms of the configuration. For this we have to chose a 
proper initial surface. Let us for convenience assume that is closed. 

Remark. Carefully note that this very last assumption implies the existence 
of a reference metric g^ with respect to which the above statement is true. 
This particular tilde here does not stand for the reference configuration but 
for the reference metric. 

Then redefine := g^vF 1 . This one- form then by definition constitutes 
a foliation of space-time by space-like hyper-surfaces S t (which we assume 
to be isomorphic to M 3 ). For those hyper-surfaces the four- velocity vF is the 
unit-normal w.r.t. g^ v . Then the constraints (j4.293p4.294j) can be written as 
(note that in this particular foliation u b = 0) 

= r - k ab k ab + k 2 - 2KT IMV u tl u v (4.320) 
= 2D [a kt } + KT^utfyb (4.321) 

Now we can state the main result of this section: 

Theorem 5. Let the stored energy depend smoothly on its arguments. Given 
a reference state (g^ u , f A ), where f A is a natural configuration such that u v : = 
g^vF is closed. Then this reference state defines a foliation of space-time 
by space-like hyper-surfaces S t . Given an initial data set (S, g^, kij, f A , f A ) 
satisfying the constraints \4 ■ 32C\4 ■ 321)) together with a diffeomorphism map- 
ping S onto the initial hyper-surface So (defined using the reference state 
(g^v, f A ) )■ Then the embedded data give rise to data (g^, d t g^, f A , d t f A )\t=o 
for the system M.SIGQ . If the latter are close to the data induced by the refer- 
ence state (g^, f A ) in the topology H S+1 (R 3 ) x # S (M 3 ) x H S+1 (M. 3 ) x H S (R 3 ) 
for s > 3 they generate a unique local solution (g^v, f A ) to the coupled system 
\4-28C\ \4 ■ 224\ ) which depends continuously on the initial data. 

For the actual function space for the solution we refer to the existence 
theorem from section 12.21 

Remark. Once the theorem is given is becomes clear what we meant by the 
phrase "close to": namely that the corresponding data on a suitable hyper- 
surface would be sufficiently close in the topology given in the theorem. An ex- 
ample would be states which are close to the reference state in the supremum- 
norm. 
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We could have tried to obtain a similar result using the corollary from 
section 12.21 If is not assumed to be closed, then one obtains local well- 
posedness. We shall formulate this the following way: 

Lemma 20. Given a reference state (g^ u , f A ), where f A is a natural config- 
uration. Then there exists a (local in space and time) foliation of space-time 
such that all data close to the the ones induces by the reference state give rise 
to a unique local (in space and time) solution for the Cauchy problem. 

4.6 Bounded domains 

In the context of matter models most naturally the question arises whether 
the equations of motion admit solutions subjected to certain boundary condi- 
tions. In this context it is again the most natural choice to consider so called 
natural boundary condition. This means that the normal stress is assumed 
to vanish. An elastic body like a planet (neglecting the atmosphere) would 
be described by such a model. Nevertheless one should stress that a planet 
like earth is still well described by the Newtonian version of the equations 
given here. On the other hand the next bigger objects, where relativistic cor- 
rections come into play, are stars, which are described well enough by fluid 
dynamics. 

Fluid dynamics can by kinematically treated as a simple subcase of an 
elastic material (the stored energy depends only on the particle density, see 
e.g. Nevertheless the techniques used in this work are not suitable to 

discuss the dynamics of fluids. This leaves us with the only physical example, 
which seems to be realized in our universe: neutron stars. Those extremely 
dense objects are known to admit solid crusts ([HE] or [23]) and it may well 
be the case that the core is solid as well. This provides a physical motivation 
for the given problem. Of course we do not treat the self-gravitating case but 
the setup discussed here has to be rather seen as a natural first step towards 
this challenge. 

There is another reason to talk about finite elastic materials, which is a 
rather pragmatic one: there is a small number of physical reasonable theories 
where a well-posed initial-boundary- value problem can be considered without 
too much trouble. From this point of view one can argue that it is reasonable 
to consider the given problem for the sole reason that we can do so. 

But instead of giving further arguments to justify the following analysis 
we shall go in medias res. We finally note that the results of this section can 
be too found in [E]. 

We will use the equations of motion in the material picture. The reader 
will find that in comparison to the foregoing sections we will place stronger 
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restrictions on the space-time geometry. The reason for this is technical and 
will be discussed "on the way". 

This section is organized as follows: First we recall the basic setup in- 
cluding the relevant equations of motion. Then we discuss suitable boundary 
conditions, which in our case happen to be natural ones. This discussion is 
then followed by some remarks concerning the connection between the Gard- 
ing inequality and strong point-wise stability as introduced in section 12.11 
Finally we will show that the corollary to the basic existence theorem due 
to Koch given in section 12.31 applies to the system under consideration, 
which gives local well-posedness of the initial-boundary value problem. 

4.6.1 Basics and equations of motion 

It is convenient for the problem at hand to use the equations of motion cast 
in the material description. The reason lies in the boundary conditions. We 
will find in the following section that while the boundary in space-time is 
free, thus the problem in that setup is that of solving a free-boundary value 
problem, it is fixed in the material picture and given through the boundary 
of the material manifold. 

In the given scenario the basic unknowns are the deformation mappings 
F, which are time-dependent mappings from a compact manifold B = BUdB 
with smooth boundary dB onto a bounded region S of the standard leaf of 
the foliation, i.e. 

F : B x R -> S (4.322) 
(X A ,t) i-> F\X,t) (4.323) 

As in the unbounded case the body is equipped with a volume-form Vabc, 
which is smooth on B. 

For this setup we can define the basic kinematic objects just the same 
way as for an unbounded medium. We thus use the definitions as given in 
section 13.21 

For technical reasons we restrict ourselves to so called "ultra-static" space- 
times. This means we assume the existence of a hyper-surface-orthogonal 
normalized time-like Killing vector field. This restriction may seem harsh 
but we stress that the sole existence of a reference deformation, which ad- 
mits a Born-rigid normalized four-velocity, already provides a quite strong 
restriction on the space-time geometry. If this reference deformation ought 
to be a solution of the field equations, we found in section EJ that this nec- 
essarily implies the existence of a normalized time-like Killing vector field. 
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For technical reasons, which will be explained at a later stage it is very 
convenient to assume this Killing field additionally to be hyper-surface- 
orthogonal. 

Under this assumption a suitable choice of coordinates allows us to write 
the metric as 

ds 2 = -dt 2 + g ij (x)dx i dx j (4.324) 

Note that this form of the metric already leads to a notable simplification 
of the formalism: Since both the lapse and the shift take their most trivial 
form N = 1, Y i = 0, we find that 



which includes 
as well as 



W* = F l (4.325) 
h AB = f A , a f B , b (g ab -F a F b ) (4.326) 



7 =(1-F 2 )- 1 (4.327) 

Remark. Remember that f A , a are to be understood as the inverse of F a ,A- 

We finally recall the Lagrangian 

L = j'hv (4.328) 

where the stored energy e is a smooth function of X A and of the strain. 
Then the field equations are 

d f dL ^ dL (4.329) 



dX" \dF\J dF i 
Thus the principal part is determined by 

^ - wB^ (433o) 

4.6.2 Boundary conditions 

We now turn to boundary conditions. In non-relativistic elasticity theory 
there are essentially two kinds of boundary conditions: 

• boundary conditions of place: the boundary of the material in space- 
time (i.e. the set f^ 1 (dB) ) is prescribed. 

• traction boundary conditions: the normal traction, i.e. the normal 
components of the stress are prescribed. 
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When the normal traction is zero, we speak of "natural" boundary condi- 
tions. For obvious reasons they are the ones appropriate for modeling a freely 
floating finite object like a star or a planet. This reason together with a sec- 
ond one strongly suggests this choice: all other boundary conditions but the 
natural ones become inconsistent once one couples to the Einstein equations. 

For the space-time description these natural boundary conditions are 
given by 

T;n v \ f -± m = (4.331) 

where T^ v is the energy-momentum tensor as given in ()3.65j) while denotes 
the co-normal of f' 1 (dB). To see this note that the material velocity is 
tangential to the inverse images of points of B under the configuration map. 
Thus it is annihilated by the co-normal. Hence the above equation becomes 
(recall that the configuration gradient is assumed to be of maximal rank) 

TABf B , v grn v \ f -i m = U (4.332) 

which is just the statement that the normal stress vanishes. 

There is however a huge disadvantage about these boundary conditions: 
the boundary f~ l (dB) is free, i.e. it is not determined a priori, but it is 
rather a part of the problem to find the actual boundary. 

This problem is overcome by switching to the material description. Then 
the boundary conditions are evaluated on dB x M, which is fixed. To actually 
obtain these new boundary conditions we have to translate ()4.332|) into the 
material description. In section 13.2.31 we have encountered the diffeomor- 
phism (t,X) = (t,f(t,x)). Applying it to the space-time co-normal n M we 
find that it is mapped to the one- form N^, which by definition is connected 
to in the following way: 

N^dX" = (N + N A f A )dt + N A f A , a dx a = (4.333) 

Recalling now that the four- velocity was mapped to a multiple of dt, we find 
that the equation n^vF 1 = now reads Nq = 0, which at once implies the 
relation 

n, = f A ,»N A (4.334) 

Thus f!4.332j) becomes 

TABh BC N c \ dB = (4.335) 

These are the new boundary conditions. They state that the normal compo- 
nents of the stress measured along the body have to vanish for all time. 

There is a deep connection between the equations of motion and the 
boundary condition. To reveal this connection we note that 

dL i dh MN 

T*Vt mn —— (4.336) 



dF\ A 1 '""dF* 
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Using formula ()4.47|) for the variation of the strain this becomes 

dL 



dF 1 



2^Vr MN f M n h NA (4.337) 



This means the following: if we replace ()4.335j) by the equivalent equation 
(recall that / was an isomorphism) 

-%y-tyf A ti r AB h BC N c \ 9B = (4.338) 

we find that we can rewrite the boundary condition in a slightly different 

way as 

<>L -N A \ dB = -^-N^\ aB = (4.339) 



dF\ A ^ dF\, 

since Nq — 0. Due to this correspondence ()4.338|) (sometimes the version 
()4.339j) ) will be the version of the boundary condition we are going to use in 
the following. 

So the problem is to find a solution to the initial-boundary value problem 
given by the equations of motion ()4.329|) subjected to the boundary condi- 
tionsdOB. 

Remark. Note that the statements of this section are independent of the 
restricted space-time geometry, but apply to the general setup. 



4.6.3 Reference deformation 

It is natural to assume the existence of a "natural state" of the material, 
which may be viewed as a trivial realization of the (abstract) material man- 
ifold in space-time. 

In the given context we assume that there exists a static reference de- 
formation, such that the elasticity operator satisfies the strong point-wise 
stability condition ()2.49|) . 

We formulate these requirements first in terms of the space-time descrip- 
tion: There exists a configuration f A , such that the four- velocity coincides 
with the Killing field: 

u = d t (4.340) 
This implies that the reference configuration is time-independent: 

dJ A = (4.341) 

From section 13.1.51 we know that the reference configuration by definition 
has to be stress-free with the corresponding normalized four-velocity be- 
ing geodesic. But the latter point together with the property of being a 
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static Killing field implies the Born condition, which in turn means that the 
reference strain can be viewed as a metric field on the body. (A simple way 
of seeing this is by recognizing that h B is time-independent as well as the 
reference configuration, which implies that it becomes a field on B, once 
translated into the material description.) 

Remark. Note that this nice property of the reference configuration is an 
artefact of this special setup. For general space-times as investigated in the 
sections dealing with the unbounded setup, there will neither exist a Born- 
rigid motion nor will the strain define a tensor field on the body. 

Remark. Note that the requirements on the elasticity operator employed here 
are stronger than those needed for the unbounded setup. We will see that it 
is a proper way of satisfying the Garding inequality, which is a necessary 
component for the Koch theorem. 

Once the properties of the reference state are known in the space-time 
description it is quite easy to determine the analogous statements for the 
material description: 

There exists a reference deformation F\ which is time- independent. The 
reference strain h AB is a metric field on the body, and satisfies the following 
conditions: 

e(X, h) =: e{X) > (4.342) 

as well as 

T AB {X,~h) = (4.343) 
Finally strong point-wise stability for 

U(X, h) A BCD =: U{X) ABC d (4.344) 

reads (suppressing the X-dependence; also compare ()2.49j) ) 

U AB cd^ AB ^ CD > Ch c{ Ah B) Dti AB ^ CD (4.345) 

for C > 0. 

Remark. We stress that strong point-wise stability implies the Legendre- 
Hadamard condition but not vice versa. 

It is worth noting that our definition of the reference deformation did not 
imply any statement on the boundary conditions. However one finds that 
the following is true: 
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Lemma 21. The reference deformation also solves the boundary value prob- 
lem l[4-329\4-33£fy . In particular this means that initial data given on the 
hyper- surface {t = 0} of the form (F A ,0) generate the reference deforma- 
tion. 

The first part of the lemma is a trivial consequence of (|4. 343(1 . The second 
part however is not that trivial. It is obvious that the reference deformation 
is a solution to (|4.329|4.338(1 with the claimed initial data. But in order to 
exclude that there are other solutions with the same initial data, which do not 
coincide with the reference deformation we have to use the uniqueness part 
of the Koch theorem. We will show in the following that the requirements of 
this theorem are met for the reference deformation (and hence by continuity 
for deformation which are close enough). This will result in the main theorem 
of this section, which implies the uniqueness statement needed here. 

The principal part given in (|4.97|4.98(1 takes a relative simple form for 
the reference deformation. (Remember that in these formulas not the prin- 
cipal part A^j was given but rather B^ B := A^F 1 ,aF\b- Another slight 
difference comes from the fact that in the definition given there the factor 
V coming from the volume form was eliminated, which is not convenient 
anymore in this context due to the presence of the boundary conditions.) 

Using the convention that every object evaluated in the reference defor- 
mation is marked by an upper tilde, f.e. 

*S\f-s = AS ( 4 - 346 ) 

we find that 

I°° = -Vega (4-347) 

as well as 

Af = WUcEDF~h Ac h DB f E n f F ,i (4-348) 

The remaining terms of the principal part vanish when taken for the reference 
deformation: 

A°f = (4.349) 

Apart from other reasons the last equations holds as a result of the special 
form of the metric. If a lapse would have been present, the elasticity operator 
would have entered and A\ A would not necessarily have vanished. 

This would have resulted in the problems we know from section 14.4.21 
concerning a proper coordinate choice. This question of coordinate choice 
gains extra-finesse since in these different coordinates the boundary would 
most probably loose its convenient coordinate form. All these problems were 
bypassed by the restrictions placed on the metric in this section. 
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We finally observe that by ()4.347j) the coefficient of the time-derivative 
Afj is negative definite for the reference deformation and thus by continuity 
for all deformation sufficiently close. 

4.6.4 Coerciveness 

A crucial entity for the Koch theorem is the coerciveness of the spatial prin- 
cipal part. At this stage the equation of state (i.e. the elasticity operator) 
comes into play. We will show that in the given context materials those stored 
energy satisfies strong point-wise stability ()4. 345(1 at the reference state give 
rise to this coerciveness condition. 

More concrete we have to show that the Garding inequality ()2.28|) 

/ Afd A fd B ^d 3 X > e\\(f>\\ 2 Hl - «| 10111a (4.350) 

JB 

holds for all smooth <p l where the constants e > and k > 0. (H 1 = W 1,2 (B) 
and L 2 = W°' 2 (B); see e.g. jUJ) 

Remark. Note that the inequality involves no derivatives of the deformation 
higher than first order. By this property it becomes a restriction on possible 
initial data. 

We will first derive the Garding inequality for the reference deformation: 
/ WUcEDFh AC ~h DB f E n f F ^d A ^d B 4P > e\\<f>\\ 2 m - k\\4\\*> (4.351) 

J B 

By continuity (|4. 350(1 is then true for all deformation sufficiently close (the 
precise form of this statement will be given later). This includes that the 
Garding inequality holds for all initial data sufficiently close to those gener- 
ated by the reference deformation. 

To prove (|4.351() we have to invoke Korn's inequality (see f.e. |28j ) or 
(to be more accurate) a slight generalization to Riemannian manifolds due 
to PH]: Given a vector field if) A on the body, which can be viewed as the 
Riemannian manifold (£>, /z-ab)- Then define L(ip) by 

L^) AB = 2hc( A ^ C ,B) (4.352) 

Then the Korn inequality tells us that there exists a positive constant C, 
such, tjlicttj 

\m)\\h + > CM\ 2 H1 (4.353) 



118 



CHAPTER 4. DYNAMICS OF ELASTIC MATERIALS 



Going back to (|4.351|) we note that by (|4. 345(1 the integral on the left hand 
side is taken over the product of two positive functions (the scalar V and the 
rest). Since V is also a continuous function on a compact set, we can estimate 
it from below by its minimum, which is positive. Using (|4. 345(1 again on the 
remaining integral we still have to show that 

/ h C (Ah B )Df B ,i (f>*, E h AE f D ,j 0V h CF (4.354) 
Jb 

is bounded from below in a suitable way. For this purpose we assume without 
restriction of generality that coordinates on B were chosen in a way such that 
F l becomes the identity map. 

Remark. It is notable that this specialization directly leads to the classical 
setup of non-relativistic elasticity theory in its classical formulation: that the 
body is identified with a certain part of space. 

For this choice of coordinates the deformation gradient coincides with the 
Kronecker delta: 

F l , A =S\ (4.355) 

which clearly includes 

f A ,i=5f (4.356) 
Applying this to the integrand in (|4. 354(1 leads us to 

/ h c[A h B)D f B n <P\ E h AE f D ^, F h CF > C\\LCf A n<p l )\\\* (4.357) 
Jb 

Note that we would obtain equality, if the norm || • \ \ L 2 was defined using 
the reference strain. If one decides to use a flat metric 5ab to construct 
these norms (as is the usual way) equality does not hold anymore. In any 
case however the equivalence of finite-dimensional norms tells us that the 
above inequality holds for some positive constant C (which depends on the 
metric used to construct the norm on the right hand side). Note that the 
proper choice of coordinates was crucial for the argument, since it allowed us 
to interchange differentiation with contraction with the inverse deformation 
gradient. 

Now we can apply the Korn inequality 1)4.353(1 to the right hand side of 
1)4.357(1 . Since f A ,i = 5f, we find that 1)4.351(1 is true, which in turn implies 
validity of 1)4.350(1 for deformations close enough. 
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4.6.5 Main result 

In this section we will show that all the requirements of the Koch theorem 
as given in section 12.31 are met in the current setup. This will ultimately 
result in a local well-posedness for the relativistic elastic initial-boundary 
value problem. 

For convenience we recall these requirements (already stated in the terms 
used here). First note that the system ()4.329l 14.339)) has exactly the re- 
quired connection between boundary conditions and equations of motion. In 
addition we have to check 

• the symmetry = A v -^ 

• hyperbolicity (and coerciveness) of the principal part 

The first point is trivially met, since our equations of motion come from 
an action principle which reflects itself in ()4.33()|) from which the desired 
symmetry can readily be read off. 

The second point consists of two parts: First we have to show that A™ is 
negative definite. For the reference deformation this is evident from (|4.347jl . 
From this we concluded the same property for all deformations "sufficiently 
close" . This statement will be sharpened in the following: 

Assume we are given initial data, i.e. a pair (F ,Fi) 6 H s+l {B) x H S {B) 
for s > 3. Then by Sobolev embedding we find that (F , F x ) e C\B) x C°(B) 
(See f.e. JHj). But this means that the deformation gradient is a continuous 
function on the initial surface {0} x B. 

Denoting the data induced by the reference deformation (which itself 
clearly is assumed to be sufficiently smooth) by (F , 0), we find that all data 
close to the reference data in the topology H S+1 {B) x H S {B) are close in 
C 1 (i3) x C°(B) (which is a point- wise statement). We conclude that there 
exists a neighborhood of the reference data in H S+1 (B) x H S (B) for which 
A®® is negative. 

On the other hand we have to show coerciveness for the spatial component 
of the principal part. By the very same argument as we used for time- 
components, we conclude that coerciveness holds for all data close to the 
reference data in H S+1 (B) x H S (B). 

We have found that all the requirements of the theorem are met, so that 
we obtain the following result: 

Theorem 6. Given a space-time and a material manifold as described in 
section \4-6'.l\ Let the stored energy be a smooth function of all its entries. In 
addition assume that there exists a smooth static reference solution as defined 
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in section \4 . 6'. c\ such that the elasticity operator satisfies the strong-point- 
wise stability condition \4-34ty - Then for all data close to the data induced 
by the reference deformation in the topology H S+1 {B) x H S {B) that satisfy 
the compatibility conditions up to order s there exists a unique T > and 
a unique solution F l e C 2 ([0,T) x B) to the initial-boundary value problem 
\4-338j) which depends continuously on the initial data. 



We end with two remarks concerning the compatibility conditions: First 
note that by the lemma 0] from section |2~^1 we know that there exists a great 
number of possible initial data obeying the compatibility conditions. 

Secondly we stress that in the proof of lemma 0] the compatibility condi- 
tions were solved by prescribing certain normal derivatives of the initial data 
along the boundary (we refer to section [2~T>1 for details). In this sense the 
boundary conditions ()4.338|) are of Neumann type. 

For details concerning the compatibility conditions we refer to section l2~3l 



Chapter 5 

Summary, Conclusions and 
Outlook 

What have we learned in the course of the thesis and in what ways can 
the results given here hope to be generalized? This will be answered in the 
following. 

5.1 Summary and conclusions 

In this thesis we treated relativistic elasticity as a Lagrangian field theory. 
The key objects were configurations (deformations respectively), which are 
maps between space-time and an abstract material manifold. During all 
of our applications we assumed the existence of a certain reference state, 
which played the role of a natural realization of the material manifold in 
space(-time). This assumption constitutes the connection to the classical 
non-relativistic setup, where the material in rest is identified with a cer- 
tain region in space. Note however that in principle the (kinetic) theory as 
presented here can be formulated without the assumption of any reference 
state. 

There remains another grave difference to the classical non-relativistic 
setup: There the material manifold is (spoken in our language) viewed as 
a Riemannian manifold equipped with a flat metric. The reference state is 
defined to be the state where the strain coincides with this given flat metric. 

For the relativistic formulation we only assume the existence of a volume 
form on the material manifold. This assumption turns out to be sufficient to 
treat kinematics. If the strain is to define a metric on the material manifold, 
there arise certain restrictions on the space-time geometry (there has to exist 
a Born-rigid vector- field) . 
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As in the non-relativistic theory the elastic properties of the material are 
basically encoded in the stored energy function, which is a measure for the 
elastic energy in dependence of the strain. The difference between relativistic 
and non-relativistic theory lies in the definition of the strain (which can be 
viewed as a metric for the deformed state) itself. While in the non-relativistic 
theory it is a purely spatial object (compare &ab as defined in ()3.106|) ) in 
the given context it is constructed out of the full space-time metric (compare 
(JHUnj)) including also time-derivatives of the configuration. 

In the second part of this thesis we derived the equations of motion gov- 
erning the time-evolution of the elastic material from an action principle 
arising from the variation of a Lagrangian which was basically given by the 
total energy density. The resulting equations were then analyzed both on the 
linearized and non-linear level. For the linearized equations on Minkowski 
space-time we could give an explicit formula for the solution of the initial 
value problem corresponding to an un-bounded homogeneous and isotropic 
material. On the non-linear level we discussed conserved quantities and the 
connection between the thus defined energy and Killing vector fields (gener- 
ating trivial spatial configurations). 

Finally we proved well-posedness results for a variety of scenarios includ- 
ing both bounded and unbounded elastic materials both on a fixed gravita- 
tional background and self-gravitating. 

We found that the initial boundary value problem required stronger re- 
strictions on the material's elastic properties than the other setups, for which 
a rather mild assumption (real propagation speeds for linearized elastic wave 
solutions) turned out to be sufficient. 

In this sense we can say that the theory as presented in this thesis is 
well-suited for the description of relativistic elastic materials. 

5.2 Outlook and open problems 

Given the results of this thesis what else remains there to do? In principle 
there are basically two ways to proceed. One may either choose to go into 
detailed analysis for certain kinds of materials (equations of state) in hope 
of obtaining some results which could be compared with observational data 
(f.e. concerning neutron stars). This can be seen as a complementary ap- 
proach to the "ab initio" attempts of obtaining suitable equations of state 
(see [23] ). It might lead to an understanding of certain effects like the well- 
known "glitches" (see 38j and references therein). This task however would 
require numerical effort rather than analytical. But even for the numerical 
treatment well-posedness of the equations of motion is a necessary ingredient 
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On the other hand one can of course carry on without numerics. Such 
analytical work could go in two directions. 

First one might try to generalize the result from section 14.61 to include 
self-gravitation. This task by far exceeds the one related to the finite body 
on fixed background, since one would have to care about the propagation 
of the constraints along the boundary. We have learned in section 14.61 that 
it is favorable to formulate the boundary value problem in the material de- 
scription. Section l4"31 however teaches us that the self-gravitating problem is 
conveniently addressed in the space-time setup. One might even ask how to 
formulate the self-gravitating problem in the material description in terms of 
partial differential equations. The diffeomorphism between space-time and 
{material manifold times time} introduced in section 13.2.31 may be of some 
help in this conquest. 

Secondly one might try to derive global results. At least for small initial 
data (i.e. data close to a natural reference state) it should be possible to 
prove global existence on Minkowski space-time. For the non-relativistic 
theory such results already exist for a growing number of scenarios. (See jUj 
for an overview.) 

There are at least two important ingredients in the proofs of such re- 
sult. One has to require a certain "null-condition", which prevents the 
self-interaction of the eigen-modes of the elastic material (see [21] for the 
definition for wave equations and jlHj for isotropic elastic equations). In 
non-relativistic theory this condition is known to be necessary for the long- 
time existence of solutions ([SB] (it should be possible to obtain analogous 
non-existence results for the present setup). 

Under the assumption of this null-condition (a generalization of this con- 
dition also plays an important role in general relativity, see [H], jlOj, |39j ) 
one then uses generalized energy estimates (first introduced by Klainerman 
in the study of non-linear wave equations |3S1) which are connected to the 
symmetry group of isomorphisms of space-time. In the non-relativistic theory 
this symmetry group is the Galilei group and hence only 6-dimensional. One 
has to employ special projection techniques (see [IB]) to obtain additional 
bounds. In the present relativistic version we have the full 10-dimensional 
Poincare group at our disposal. This might lead to a simplification of the 
task of obtaining the required estimates. One can thus hope that it will be 
rather straightforward to obtain a global existence result for the relativistic 
setup by following the non-relativistic proof. 

The ultimate challenge would be the derivation of statements on the 
asymptotic structure for elastic materials. For unbounded materials one 
might expect that the global solution for small initial data approaches the 
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unstressed reference solution (on any open neighborhood) as time passes on 
(perturbation due to non-trivial initial data travel to infinity). The analo- 
gous problem for finite materials implies further problems. Will the material 
approach a natural state? What about translations of the center of mass? 
How to describe such a phenomenon? What is to be expected in the self- 
gravitating setup? 

Finally (far beyond the scope of any stern approach): what about global 
asymptotic behavior for a self-gravitating finite elastic material? What is 
there to be expected? 

This is a whole load of interesting and challenging questions, both from 
mathematical and physical point of view. In this sense the results obtained 
up to now can be seen as a (necessary) first step on the way to answering 
the questions of the type stated above. 



Chapter 6 
Appendix 



This appendix contains some minor results which do not directly fit into the 
main body of this work. It contains an example which is incorporated to 
become familiar with the objects defined in chapter El as well as remarks on 
energy conditions, Newtonian limit of the field equation and finally some 
formulas arising in the context of the 3+1 decomposition. 

6.1 An easy example: 1 + 1 dimensions 

This section contains a simple 1+1 dimensional example to get used to the 
definitions introduced in chapter El For convenience we will specialize to 
static space-times, in particular we concentrate on the elastic equations on 
Schwarzschild background. Note however that we do not employ any " small- 
stress" assumptions. 

6.1.1 Basics 

The unknown is a mapping / from two dimensional space-time M onto a one 
dimensional manifold B called body. The body is equipped with a volume 
form VdX, which is a one-form in the given context. 



We assume the tangent mapping df to be of maximal rank with time-like 
null-space. In the given context this simply means df 7^ 0. Therefore (as 
long as M is time-orientable) there exists a unique future-pointing time-like 
vector-field satisfying the two conditions 



/: M 



x a X = f(x a ) 



B 



(6.1) 
(6.2) 



u a f ta = g(u,u) = -l. 



(6.3) 
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The strain hdXdX is determined by the function h which (in accordance 
with (|3.10)l ) is given by 

h = UUg> lu (6.4) 
Then the particle density becomes 



n 



vVh (6.5) 



We also assume the existence of a reference configuration /. Then we can 
choose coordinates on the body in such a way that h = 1. For convenience 
we assume that V = 1 in these coordinates, such that the reference strain is 
a flat metric on the body and V is the corresponding volume form. 

Then a material which is homogeneous and isotropic w.r.t. this reference 
configuration is given by a stored energy which depends only on the invariants 
of the strain (taken w.r.t. the reference strain, which is a flat metric). Since 
in one dimension there is only one independent invariant we have: 

e(n) = e{Vh) (6.6) 

Remark. In this sense there is no difference between fluids and other elastic 
materials in this low dimension. The stored energy does (in the isotropic 
case) only depend on the particle number density. 

Consequently we have that the stress is given by 

T dh 2ndn 1 ' ' 

For the energy-momentum tensor this implies 

de 

= pu^u u + -T^fw f,u (6-8) 

From this point on we will specialize to static setups. What we exactly mean 
by this we explain in the following: 

6.1.2 Static setup 

We assume the existence of a Killing field which in addition is assumed 
to be proportional to the four-velocity: 

u = ^ (6.9) 

where 

<^ 2 :=(-<?(£, (6.10) 
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Under these assumptions we can easily derive the equations of motion. We 
have seen in section 14. 1 . II that they are equivalent to divergence- less-ness of 
the energy-momentum tensor. Recalling the continuity equation (J3.8)) this 
condition is equivalent to 

de_ 
dn' 



= nu»(eu u y, li +[—ffU ;„ (6.11) 



Since £^ is a Killing field, we have =0. As a consequence u^\ v = 

4>ivi^ + <Pi^]u yields -u^ = 0. We also find that u^u v ]^ = \4> 2 ((j)~ 2 ) w . From 
this we find that ()6.11|) contracted with the four-velocity leads to 

= -™%,+^//*/;^ ( 6 - 12 ) 



uV ^=~— fffW (6.13) 



But since 

1 de 
n dn' 

we find that T^;^u u = is satisfied identically by virtue of the continuity 
equation as was to be expected from section 13.1.41 

The remaining component is the only non-trivial one. It is obtained by 
contraction of (|6.11J) with f, u . 

= -neuWf^ + (^f» f, v ) ;„ f," (6.14) 

The bracket term includes both first and second variations of the stored 
energy. Besides these non-vanishing stress-terms the above equation looks a 
lot like the equations of motion as given in section 14.1.11 They constitute a 
second order equation for the configuration /. In fact it isn't even a PDE 
but can rather be cast in the form of an ODE, since we assumed the material 
velocity to be proportional to a Killing field. To see this more clearly we again 
specialize our example by fixing a particular background metric, namely the 
Schwarzschild metric. 

6.1.3 Elasticity on Schwarzschild background 

We now use the Schwarzschild metric 

2m\ , 9 / 2m \ 1 , 9 
1 i / 1 l 



ds = ~{ 1 ~ ~v) dt vj ( 6 ' 15 ^ 

Therefore we use the following Killing field 

£ = d t (6.16) 



2m x 



(6.17) 
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The first consequence is that the configuration is time- independent, i.e. 

/,*=e7,M=<rvv, M =o (6.i8) 

From this we find that 

u"«7;^=-^" 2 ),r/,r (6.19) 

Thus ()6.14j) becomes a single ordinary differential equation in r. 
We find that the second term in (|6.14j) becomes 

^- 2 f, 2 r ) ,r f,r ^ + ~f^. ,r (6.20) 

Altogether this means that ()6.14|) is given by (using the identity n 2 = h = 

= \n{r%rfrr (e + n^j + {^n^ , r nf 1 (6.21) 
Multiplying by <\>n~ x leads to 

O = (ln0- 1 ), r n^e + n^^ + (^n 2 ) , r (6.22) 

Assuming that (e + nj^) ^ 0, we conclude that ()6.14|) is equivalent to 

= (ln^" 1 ),^ MnU + n^JJ , r (6.23) 

Thus we can solve by integration. We obtain 

de 

e + n— = C0 (6.24) 
on 

where the constant C has to be determined via the boundary conditions 
entering the integration. This is an algebraic equation for n and hence for 
f, r . If we can solve algebraically for the configuration gradient, we can obtain 
the configuration itself via integration. 
For convenience we define 

A(n) := e + ne m = p m (6.25) 

Using the definition for the stress 
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we can translate A into 



A{n) 



P + P 



n 



(6.27) 



Recall that we used the assumption A ^ 0. In fact we will also assume that 
A is positive: 

A > (6.28) 



and monotone 



dA - 

dn 



(6.29) 



This implies that the energy density p has to be convex. 

We now turn to the analysis of ()6.24|) for certain special setups. 

First we concentrate on the case of a hanging body (supported at the 
upper end). 



[nun) 



It is natural to assume r u and rii to be given. The first is the suspension 
point. The second is given through p(ni) = 0, since we require the pressure 
to vanish at the lower end. Note that we can solve p{ni) = for ni, since by 



dp 
dn 



n 



dA 

dn 



(6.30) 



the pressure p is monotone in n. For the described scenario (j6.24j) becomes 



A(n u ) = A^^r 1 ^) =: d 1 



2m\ 2 

n 



(6.31) 



where C\ > 0. This equation only gives us information on the endpoints of 
the material. 

To obtain a formula which gives the particle density for general r (away 
from the endpoints), we assume for the moment that r\ is given (this assump- 
tion has to be removed later on). Then ()6.24|) becomes 



A(n) = A(n z )0" 1 (r i )0(r) 



(6.32) 
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Assume for a moment that m = (i.e. vanishing gravity). Then n = rii 
which means that the particle density is constant along the stick and thus 
that there is no pressure present when gravity is turned off. 

Remark. Recall that m was uniquely determined by the condition p{n{) = 0. 

Since we are interested in a finite body we assume that its intrinsic length 
L defined by 

L = [ V = [ dX (6.33) 

J B JB 

is finite. This implies (the configuration is viewed as mapping between body 
and the space orthogonal to the Killing flow) 

L= f*V= f,dr> = / —mdr> (6.34) 

Jf-HB) J n Jrt (1-^) V 

We can use this equation to get rid of the ri entering formula ()6.32|) . We will 
see that we can solve ()6.34j) for ri in terms of r u , ni, L and m. 
For this we have to show that 

dL r Tu dn 
^ — = / 4>{r)—dr - n{n)(j)(ri) (6.35) 

on J ri vn 

This condition is necessary and will give rise to a restriction on A and hence 
on the stored energy e. 

To see this note that the second term is negative while for the first term 
fEHU) gives 

^ = A{nM{r fJ^ = A{nmrd ^ (6.36) 

I I I 

which is positive. By the monotonicity of A (recall |^ > 0) this includes the 
analogous statement on the variation of n. Thus from our original (mild) 
assumptions on A we cannot conclude a sign for the right hand side of ()6.351 
We thus assume the validity of ()6.35|) in the following. It is reasonable to 
assume that ^ to be negative, since this implies that a longer material is 
also longer in space. In a more explicit form this condition reads 

dL C Tu A m 
0> 7T= ^^{nf-dr - m<f>(n) (6.37) 
Ofi J n A, n r l 

Without analyzing this condition any further at the moment we shall assume 
it as given, since it obviously plays the role of a minimal requirement for a 
reasonable model. 
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To assign a coordinate- independent meaning to a certain length in space(- 
time), we define an intrinsic length using the induced volume element of 
t = const surfaces of space-time: 



dr 



2m W2' 
r 



(6.38) 



Since both r u and r; can be obtained from the data, we can assign a spatial 
length to the material. 

One would expect that the elastic material is stretched when gravity be- 
comes stronger. To actually see this we compute l m , m . For this computation 
recall that n depends on m via ()6.34|) . Consequently 



dr 



n r 



2m \ 3 / 2 



2m 

n 



1/2 



(6.39) 



The first term is positive, which is in accordance with what one would expect. 
In the second term we have to get rid of the m derivative of r/, since we know 
a priori nothing about its sign. This is done by differentiating ()6.34j) with 
respect to m and using the resulting equation to determine the sign of ri, m . 
Since L is by construction independent from m, we have 



dL 

= -K-r hm + 

on 



n 



+ 



A 

A,r, 



4? <p(nf 



n 



dr 



(6.40) 



To obtain the desired result, namely that r;, m < we have to show that the 
integral has is negative. 

Since doing these estimates for general m would go beyond the scope of a 
simple example we shall be satisfied with a simpler case, namely by showing 
the above claims for small m. For this first note that ()6.37j) is indeed satisfied: 



dL, 
dri 



m=0 



-n, < 



(6.41) 



By continuity the above inequality and the following results also hold for 
small m. 

But we can do better. In fact we can derive an explicit formula for r\. 
We can use that the solution for the elastic equation ()6.32j) for m = implies 
that the particle density has to be constant to evaluate ()6.34|) at m = 0: 



H m=0 — r v 



L 

ni 



(6.42) 
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Thus in order to obtain an explicit expression for ()6.39|) for m = 0, we only 
have to compute ti,m \ m =o This is done by evaluating ()6.40j) at m — 0. 



^Urn |m=0 



1 

- + 



.4 



nA m 



\n=ni 



1 



dr 



(6.43) 



Thus according to ()6.39p we have 

A 



I 



mim |m=0 



nA, 



\n=ni 



) dr 

ri r 



(6.44) 



Since the integrand is positive we conclude that there an actual growth in 
the length of the material in space at least in first order. 

Using the above results we can expand the length of the elastic material 
into a formal Taylor series in the Schwarzschild mass m. Using only the given 
data, which are (n ; , r u ) as well as A and L we obtain 



dr A 
h m — — 

l r nA,, 



\n=ni 



1 



r,. - 



- - dr + h.o.t. (6.45) 



■in 



The next example is simpler in nature, since we need not use l m . Assume 
now to hang the body into the gravitational field. Then place the same body 
such that its lower end stands, where the hanging body ended. 



y^ui Tu ) 



(n ,ri) 



{n , R u ) 



{n h Ri = n) 



We may now compare the two solutions corresponding to the picture 
(Compare (|6.31j) ). where n/j(r) denotes the particle density along the hanging 
body (compare ()6.32j) ). while n s {r) corresponds to the hanging body (note 
that in this setup now no defines the pressure-free state): 



A{n h )(j){Ti) = A{n )(f) 
A(n s )(f)(R u ) = A(n )<p 



(6.46) 
(6.47) 



These equations give rih and n s as functions of r. To prove that the standing 
material does not reach the original suspension point of the hanging material 
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given by r u , we proceed by contradiction. Since both scenarios n s and 
correspond to the same material, they both lead to the same intrinsic length 
L. Assume now that Rjj > r u . Since the integrand is positive this includes 

pr u pRu pr u 

/ nh4>dr = I n s (j)dr > / n s (j)dr (6.48) 

Jri Jr t Jr t 

This means that there have to exist at least some points f where nh(f) > 
n s (f). By the monotonicity of A this includes A(rih{r)) > A(n s (r)). 

Now we make use of the equations ()6.46|) and ()6.47|) . Taken together they 
imply 

A{n h )4>{n) = A(n s )<P(R u ) (6.49) 
Combining this with the last observation we obtain 

<P(R U ) > (6.50) 

This is equivalent to 

R u < n (6.51) 

which is a contradiction to the setup. Thus we have proved the claim that 
r u > R u and have thus shown that the hanging material is the longer one. 

Remark. Note that one can formulate the above problem in a coordinate 
invariant way by comparing the spatial length l m . But it is easy to see that 
this problem can be reduced to the one handled here. 

Remark. Note that in these examples we had no grave restrictions on the 
equation of state nor on the state itself. In contrary to the main body of this 
work the results derived here hold for materials subjected to high stress. 



6.2 A short remark on energy conditions 

This section will be concerned with the energy conditions and whether they 
hold for elastic materials. This will be a rather short discussion since the 
energy- momentum tensor given by (see equation 1)3.65}) ): 

T^ v = pu^Uy + 2nT AB f A , fl f B ,u (6.52) 

includes the stress. Our conditions did only involve the elasticity operator. 
But the following is trivially true: for states sufficiently close to a reference 
state, the energy-momentum tensor is arbitrarily close to that of a pressure- 
free perfect fluid. 

In this sense all the energy conditions hold. For a stress-free state we 
have 
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• the weak energy condition: 

= p{u^f > o (6.53) 

for all time- like . 

• the dominant energy condition, since in addition to the weak energy 
condition 

-T;V» = -p(u^)u v (6.54) 
is future-pointing and time-like for time-like since V^u^ < 0. 

• the strong energy condition: for all with V^V 1 = —1 we have 

Wir = P {u ■ vf > p(-u 2 K-v 2 ) = P = -t^sT > -\t£ 

Remark. By the above observations at least some of the states appearing in 
our analysis in the last chapter obey the energy conditions. 

It may be interesting to see to which degree the states for which local exis- 
tence can be obtained satisfy energy conditions (and which of them in partic- 
ular). Both statements use smallness assumptions on the stress. Formulating 
those in a proper way may (or may not) lead to interesting agreement. 



6.3 Non-relativistic limits 

Here we show that the non-relativistic limit of the elastic field equations 
coincides with the known non-relativistic field equations. Since the non- 
relativistic limit for the space-time description can be found in [5], we will 
concentrate on the material description. The non-relativistic field equations 
are given by (see f.e. 



dX A \dF\ 

which can be obtained from an action principle of the form (see again 

/ ( ~\f 2 + Kk) J dtd 3 X (6.56) 

We will do the non-relativistic limit only for the Lagrangian, since the cor- 
responding limit for the equations of motion is then a trivial consequence if 
the involved quantities are regular. 



6.3. NON-RELATIVISTIC LIMITS 



135 



The action for the material description according to ()3.113|) is given by 
S[F] = - I Nj-hVdtd 3 X (6.57) 

Note the modification: we inserted an extra - to keep the limit regular. 
Otherwise the total Lagrangian would diverge with order c. 

In a first step we set V — 1. The metric is chosen to be Minkowskian, i.e. 

r) fiU dx' J 'dx l/ = —c 2 dt 2 + 5ijdx l dx^ (6.58) 

Note that this choice already implies that the foliation, w.r.t. to which 
the material description was introduced (recall that such a foliation was 
necessary) is flat. 

By the above restrictions N = c, while the shift vanishes. We find that 

j- 1 = 1 - W 2 = 1 - ^F i 6 ij F j (6.59) 

Thus developing into a formal Taylor series in terms of d := - we find that 

7 -2 = 1 - d 2 ]-F 2 + h.o.t. (6.60) 

It remains to evaluate the stored energy. For this purpose we have to say a 
word on the strain, which is given by 

h AB = k AB _ W A W B (6 01) 

where k AB = f A , a S ab f B ,b is the non-relativistic version of the strain, while 
w A = df A , a F a (the d comes from the fact that we use t now instead of ct as 
time variable). 

Note that by the expression for w B the strain is of order 

h AB = k AB + Q ^ (a62) 

The stored energy consists of a rest-energy portion and a variable part e 

e(h) = c 2 + 1(h) (6.63) 
Using the above observation on the strain, we find that 

e(h) = c 2 + e(k) + 0(d 2 ) (6.64) 
Thus the Lagrangian becomes 

C = (c 2 + e(fc))(l - d 2 ^F 2 ) + 0(d 2 ) = (c 2 - F 2 + e{k)) + 0(d 2 ) (6.65) 
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There are essentially three parts: the first part is constant and thus does 
not contribute to the field equations, since it vanishes under variations. The 
last term vanishes in the non-relativistic limit d —>■ (which is of course 
equivalent to c — > oo). Thus the action one obtains from the non-relativistic 
limit is 

Snonrel[F}= [ ( - F 2 + e(k)) dtd? X (6.66) 

which obviously coincides with the action given above. 

Remark. Comparing this non-relativistic Lagrangian with the relativistic 
one, we find the following: if one linearizes the above Lagrangian, the re- 
sulting field equations would exactly coincide with the ones derived from the 
linearized material Lagrangian fl^.i^6| ). There is thus no way to distinguish 
between relativistic and non-relativistic elasticity on the linearized level. 



6.4 3+1 decompositions 

At some stage of this work (section f3.2.3|) we needed the form of the Christof- 
fel symbols in a 3+1 decomposition. For matters of completeness these shall 
hence be given at this point. 

Let the metric be given in the ADM form 

ds 2 = -N 2 dt 2 + g ij (dx i + Y i dt)(dx j + Y j dt) (6.67) 

Then the inverse metric reads 

g^d^ = -N- 2 (d t - Y l d t ) 2 + jididj (6.68) 

Then one can compute the Christoffel symbols in terms of lapse, shift and 
induced metric. The result is 

2iV 2 r° = —Q t {-N 2 + Y 2 ) + r {-{-N 2 + Y 2 ), { +2Y^j (6.69) 

2iV 2 r°, = -(-N 2 + F% +Y l (2Y [hj] + 9lj ) (6.70) 

2iV 2 r°. = 9ij - 2Y (i , 3) -Y l g ivl +2Y l g l(hj) (6.71) 

d 

2N 2 r 00 = Y l — (-N 2 + Y 2 ) + (N 2 g d - Y l Y l ) ((N 2 - Y 2 ),i +2Y t ) (6.72) 

2iV 2 r 0j = -Y\N 2 - Y 2 ) n + (N 2 g a - Y { Y l ) {2Y [hj] + 9lj ) (6.73) 
2N 2 T) k = -Y% k + 2Y%, k) + (N 2 g u - Y*Y l ) {-g jM +2g l{k , 5) ) (6.74) 
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